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Preface 

These being lecture notes for a summer school, one should not seek original material in them. Rather, 
the most one could hope to find would be the insight arising from incorporating a unified approach 
(based on reduction by symmetry of Hamilton's principle) with some novel applications. I hope the 
reader will find insight in the lecture notes, which are meant to be informal, more like stepping stones 
than a proper path. 

Many excellent encyclopedic texts have already been published on the foundations of this subject 
and its links to symplectic and Poisson geometry. See, for example, [AbMal978], [Arl979], [GuStl984], 

[JoSa98], [LiMal987], [MaRal994], [McSal995] and many more. In fact, the scope encompassed by the 
modern literature on this subject is a bit overwhelming. In following the symmetry-reduction theme 
in geometric mechanics from the Euler-Poincare viewpoint, I have tried to select only the material 
the student will find absolutely necessary for solving the problems and exercises, at the level of a be- 
ginning postgraduate student. The primary references are [Mal992], [MaRal994], [Le2003], [B12004], 
[RaTuSbSoTe2005] . Other very useful references are [ArKhl998] and [012000] . The reader may see 
the strong influences of all these references in these lecture notes, but expressed at a considerably 
lower level of mathematical sophistication than the originals. 

The scope of these lectures is quite limited: a list of the topics in geometric mechanics not included 
in these lectures would fill volumes! The necessary elements of calculus on smooth manifolds and the 
basics of Lie group theory are only briefly described here, because these topics were discussed in more 
depth by other lecturers at the summer school. Occasional handouts are included that add a bit more 
depth in certain key topics. The main subject of these lecture notes is the use of Lie symmetries in 
Hamilton's principle to derive symmetry-reduced equations of motion and to analyze their solutions. 
The Legendre transformation provides the Hamiltonian formulation of these equations in terms of 
Lie-Poisson brackets. 

For example, we consider Lagrangians in Hamilton's principle deflned on the tangent space TG of 
a Lie group G. Invariance of such a Lagrangian under the action of G leads to the symmetry-reduced 
Euler-Lagrange equations called the Euler-Poincare equations. In this case, the invariant Lagrangian 
is defined on the Lie algebra of the group and its Euler-Poincare equations are defined on the dual Lie 
algebra, where dual is defined by the operation of taking variational derivative. On the Hamiltonian 
side, the Euler-Poincare equations are Lic-Poisson and they possess accompanying momentum maps, 
which encode both their conservation laws and the geometry of their solution space. 

The standard Euler-Poincare examples are treated, including particle dynamics, the rigid body, 
the heavy top and geodesic motion on Lie groups. Additional topics deal with Fermat's principle, 
the Poisson bracket, polarized optical traveling waves, deformable bodies (Riemann ellipsoids) and 
shallow water waves, including the integrable shallow water wave system known as the Camassa-Holm 
equation. The lectures end with the semidirect-product Euler-Poincare reduction theorem for ideal 
fiuid dynamics. This theorem introduces the Euler-Poincare variational principle for incompressible 
and compressible motions of ideal fluids, with applications to geophysical fluids. It also leads to their 
Lie-Poisson Hamiltonian formulation. 

Some of these lectures were first given at the MASIE (Mechanics and Symmetry in Europe) summer 
school in 2000 [Ho2005]. I am grateful to the MASIE participants for their helpful remarks and 
suggestions which led to many improvements in those lectures. For their feedback and comments, I 
am also grateful to my colleagues at Imperial College London, especially Colin Cotter, Matthew Dixon, 
J. D. Gibbon, J. Gibbons, G. Gottwald, J.T. Stuart, J.-L. Thiffeault, Cesare Tronci and the students 
who attended these lectures in my classes at Imperial College. After each class, the students were 
requested to turn in a response sheet on which they answered two questions. These questions were, 
"What was this class about?" and "What question would you like to see pursued in the class?" The 
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answers to these questions helped keep the lectures on track with the interests and understanding of 
the students and it enfranchised the students because they themselves selected the material in several 
of the lectures. 

I am enormously grateful to many friends and colleagues whose encouragement, advice and support 
have helped sustained my interest in this field over the years. I am particularly grateful to J. E. 
Marsden, T. S. Ratiu and A. Weinstein for their faithful comraderie in many research endeavors. 
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1 Introduction 

1.1 Road map for the course 

• Spaces - Smooth Manifolds 

• Motion - Flows (pt^ (l>s = 4>t+s of Lie groups acting on smooth manifolds 

• Laws of Motion and discussion of solutions 

• Newton's Laws 

— Newton: dp/dt = F, for momentum p and prescribed force F (on historically) 

— Optimal motion 

* Euler-Lagrange equations - optimal "action" (Hamilton's principle) 

* Geodesic motion - optimal with respect to kinetic energy metric 

• Lagrangian and Hamiltonian Formalism 

— Newton's Law of motion 

— Euler-Lagrange theorem 

— Noether theorem 

— Euler-Poincare theorem 

— Kelvin-Noether theorem 

• Applications and examples 

— Geodesic motion on a Riemannian manifold 

— Rigid body - geodesic motion on 50(3) 

— Other geodesic motion, e.g., Riemann ellipsoids on GL{3,R) 

— Heavy top 

• Lagrangian mechanics on Lie groups &; Euler-Poincare (EP) equations 

— EP(G), EP equations for geodesies on a Lie group G 

— EPDiff(M) for geodesies on Diff(M) 

— Pulsons, the singular solutions of EPDiff(M)) wrt any norm 

— Peakons, the singular solitons for EPDiff(M, i?^), wrt the norm 

— EPDiff(M") and singular geodesies 

— Diffeons and momentum maps for EPDiff(M") 

• Euler-Poincare (EP) equations for continua 

— EP semidirect-product reduction theorem 

— Kelvin-Noether circulation theorem 

— EP equations with advected parameters for geophysical fluid dynamics 
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Hamilton's principle of stationary action: 



Lagrangians on TM'^^; 
Euler-Lagrange equations 
Noether's theorem 
Symmetry =^ cons, laws 
Legendre transformation 
Hamilton's canonical equations 
Poisson brackets 
Symplcctic manifold 
Momentum map 
Reduction by symmetry 



G -invariant Lagrangians on TG: 

Euler-Poincare equations 
Kelvin-Noether theorem 
Cons, laws are built-in 
Legendre transformation 
Lie-Hamilton equations 
Lie-Poisson brackets 
Poisson manifold 
Momentum map 
Reduction to coadjoint orbits 



1.2 Motivation for the geometric approach 

We begin with a series of outline sketches to motivate the geometric approach taken in the course and 
explain more about its content. 

Why is the geometric approach useful? 

• Defines problems on manifolds 

— coordinatc-frcc 

* don't have to re-do calculations when changing coordinates 

* more compact 

* unified framework for expressing ideas and using symmetry 

• "First principles" approach 

— variational principles 

— systematic - unified approach 

e.g. similarity between tops and fluid dynamics (semi-direct product), and MHD, and . . . 

— POWER 

Geometric constructions can give useful answers without us having to find and work with 
complicated explicit solutions, e.g. stability of rigid body equilibria. 

Course Outhne 

• Geometrical Structure of Classical Mechanics 

— Smooth manifolds 

* calculus 

* tangent vectors 

* action principles 

— Lie groups 

* flow property 4>t+s = (l)t°(t>s 

* symmetries encode conservation laws into geometry 

* richer than vector spaces 

— Variational principles with symmetries 

* Euler-Lagrange equations — > Euler-Poincare equations (more compact) 

* Two main formulations: 
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Lagrangian side: 
Hamilton's principle 
Noether's theorem 
symmetry =^ cons, laws 
momentum maps 

(These two views are mutual beneficial!) 
Applications and Modelling 

- oscillators Sz resonance (e.g., LASER) 

- tops - integrable case 

— fluids 

shallow water waves 

— waves < optical pulses 
solitons 



Hamiltonian side: 
Lie-Poisson brackets 
cons, laws symmetries 
momentum maps 
Jacobi identity 



Range of topics 

Rigid body 

• Euler-Lagrange and Euler-Poincare equations 

• Kelvin-Noether theorem 

• Lie-Poisson bracket, Casimirs &; coadjoint orbits 

• Reconstruction and momentum maps 

• The symmetric form of the rigid body equations (Q = QCl, P = PCI) 

• bracket and intersecting level surfaces 

± = VC xVH = V{aC + /3H) x V{jC + eH) , for ae - /?7 = 1 

Examples: 

(1) Conversion: rigid body <^=^ pendulum, 
(3) Fermat's principle and ray optics, 

(2) Self-induced transparency. 

• Nonlinear oscillators: the n : m resonance 

• SU (2) rigid body, Cayley-Klein parameters and Hopf fibration 

• The Poincare sphere for polarization dynamics 

• 3-wave resonance, Maxwell-Bloch equations, cavity resonators, symmetry reduction and the Hopf 
fibration 

• 4-wave resonance, coupled Hopf fibrations, coupled Poincare spheres and coupled rigid bodies 

• Higher dimensional rigid bodies 

- Manakov integrable top on 0(n) and its spectral problem 
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• Semi-rigid bodies - geodesic motion on GL(3) and Riemann ellipsoids 

• Reduction with respect to subgroups of GL{3) and Calogero equations 

Heavy top 

• Euler-Poincare variational principle for the heavy top 

• Kaluza-Klein formulation of the heavy top 

Utility 

Kirchhoff elastica, underwater vehicles, liquid crystals, stratified fiows, polarization dynamics of telcom 
optical pulses 

General theory 

• Euler-Poincare semidirect-product reduction theorem 

• Semidirect-product Lie-Poisson formulation 

Shallow water waves 

• CH equation - peakons (geodesies) 

• EPDiff equation - (also geodesies) 

Fluid dynamics 

• Euler-Poincare variational principle for incompressible ideal fluids 

• Euler-Poincare variational principle for compressible ideal fluids 

Outlook: The variational principles and the Poisson brackets for the rigid body and the heavy top 
provide models of a general construction associated to Euler-Poincare reduction with respect to any 
Lie group. The Hamiltonian counterpart will be the semidirect-product Lie-Poisson formulation. We 
will often refer to the rigid body and the heavy top for interpretation and enhanced understanding of 
the general results. 

2 Review Newton, Lagrange &; Hamilton 

• Newton's Law: m,q = F{q,q), inertial frames, uniform motion, etc. 

• Lagrange's equations: = f| for Lagrangian L{q,q,t). 
Defined on the tangent bundle^ 

TQ of the configuration space Q with coordinates {q, q) G TQ, the solution is a curve (or 
trajectory) in Q parameterized by time t. The tangent vector of the curve q{t) through each 

^The terms tangent bundle and cotangent bundle are defined in Section 5. For now, we may think of the tangent 
bundle as the space of positions and velocities. Likewise, the cotangent bundle is the space of positions and momenta. 
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point q G Q is the velocity q along the trajectory that passes though the point q at time t. This 
vector is written q eTqQ. 

Lagrange's equations may be expressed compactly in terms of vector fields and one-forms (differ- 
entials). Namely, the Lagrangian vector field Xl = + F{q,q)^ acts on the one-form dq) 
just as a time-derivative does, to yield 

d/dL^\ / d 9L\ _^ d dL dL 

dt\dq ) \dt dq) \ dq J dt dq dq 



Hamiltonian H{p ■ q) = pq — L and Hamilton's canonical equations: 

. _ dH . _ dH 

^ dp ^ ^ dq 

The configuration space Q has coordinates q E Q. Its phase space, or cotangent bundle T*Q 

has coordinates {q,p) G T*Q. 

Hamilton's canonical equations are associated to the canonical Poisson bracket for functions 
on phase space, by 

r ... r ... N ... dF dH dF dH 

The canonical Poisson bracket has the following familiar properties, which may be readily veri- 
fied: 

1. It is bilinear, 

2. skew symmetric, {F , H} = — {H , F}, 

3. satisfies the Leibnitz rule (chain rule), 

{FG, H} = {F, H}G + F{G, H} 

for the product of any two phase space functions F and G, 

4. and satisfies the Jacobi identity 

{F,{G, H}} + {G,{H, F}} + {H,{F,G}} = 

for any three phase space functions F, G and H. 

Its Leibnitz property (chain rule) property means the canonical Poisson bracket is a type of 
derivative. This derivation property of the Poisson bracket allows its use in defining the Hamil- 
tonian vector field Xh, by 

' dp dq dq dp ' 

for any phase space function H. The action of Xjj on phase space functions is given by 

P = XhP. q = XHq, and Xh{FG) = {XhF)G + FXhG = FG + FG . 

Thus, solutions of Hamilton's canonical equations arc the characteristic paths of the first order 
linear partial differential operator Xh- That is, Xh corresponds to the time derivative along 
these characteristic paths, given by 

dt - - (2 1) 

" dH/dp ~ -dH/dq ^ ^ 

The union of these paths in phase space is called the flow of the Hamiltonian vector field Xh- 
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Proposition 2.1 (Poisson bracket as commutator of Hamiltonian vector fields) The Pois- 
son bracket {F , H} is associated to the commutator of the corresponding Hamiltonian vector 
fields Xp and Xh by 

^{F,H} = XhXf — XpXn ='■ — [Xp , Xh] 
Proof. Verified by direct computation. ■ 

Corollary 2.2 Thus, the Jacohi identity for the canonical Poisson bracket {• , •} is associated 
to the Jacobi identity for the commutator [■ , •] of the corresponding Hamiltonian vector fields, 

[Xp , [Xg , Xh] + [Xq , [Xh , Xp] + [Xh , [Xp , Xg] = . 

Proof. This is the Lie algebra property of Hamiltonian vector fields, as verified by direct 
computation. ■ 

2.1 Differential forms 

The differential, or exterior derivative of a function F on phase space is written 

dF = Fqdq + Fpdp , 

in which subscripts denote partial derivatives. For the Hamiltonian itself, the exterior derivative 
and the canonical equations yield 

dH = Hqdq + Hpdp = — pdq + qdp . 

The action of a Hamiltonian vector field Xh on a phase space function F commutes with its 
differential, or exterior derivative. Thus, 

diXnF) = XnidF) . 

This means Xh may also act as a time derivative on differential forms defined on phase space. For 
example, it acts on the time-dependent one-form pdg(t) along solutions of Hamilton's equations 
as, 

Xnipdq) = ^{pdq) = pdq + pdq 

= pdq — qdp + d{pq) 

= — Hqdq - Hpdp + d{pq) 

= d{-H + pq)=:dL{q,p) 

upon substituting Hamilton's canonical equations. 

The exterior derivative of the one-form pdq yields the canonical, or symplectic two-form^ 

d{pdq) = dp hdq 

^The properties of differential forms are summarized in the handouts in sections 3 and 15. 
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Here we have used the chain rule for the exterior derivative and its property that (P = 0. (The 
latter amounts to equality of cross derivatives for continuous functions.) The result is written 

in terms of the wedge product A, which combines two one- forms (the line elements dq and dp) 
into a two- form (the oriented surface element dp A dq = — dq A dp). As a result, the two- form 
Lo = dq A dp representing area in phase space is conserved along the Hamiltonian flows: 



This proves 

Theorem 2.3 (PoincEire's theorem) Hamiltonian flows preserve area in phase space. 

Definition 2.4 (Symplectic two- form) The phase space area u = dq Adp is called the sym- 
plectic two-form. 

Definition 2.5 (Symplectic flows) Flows that preserve area in phase space are said to be 
symplectic. 



Xh {dq A dp) 



d^ 
dt 



(dq Adp) =0 



Remark 2.6 (Poincare's theorem) Hamiltonian flows are symplectic. 
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3 Handout on exterior calculus, symplectic forms and Poincare's 
theorem in higher dimensions 

Exterior calculus on symplectic manifolds is the geometric language of Hamiltonian mechanics. As an 
introduction and motivation for more detailed study, we begin with a preliminary discussion. 

In differential geometry, the operation of contraction denoted as _l introduces a pairing between 
vector fields and differential forms. Contraction is also called substitution of a vector field into a 
differential form. For example, there are the dual relations, 

dq -1 dq = 1 = dp -1 dp , and dq A dp = Q = dp A dq 

A Hamiltonian vector field: 

XH = q^^+P^^ = Hpdq-Hqdp = {-,H} 

satisfies 



Xh -Jdq = Hp and Xh -Jdp= — Hq 

The rule for contraction or substitution of a vector field into a differential form is to sum the substitu- 
tions of Xh over the permutations of the factors in the differential form that bring the corresponding 
dual basis element into its leftmost position. For example, substitution of the Hamiltonian vector field 
Xh into the symplectic form uj = dq h dp yields 

Xh -^LJ = Xh -^{dq ^ dp) = {Xh -I dq) dp - {Xh -I dp) dq 

In this example, Xh -idq = Hp and Xh -i dp = — Hg, so 

Xh -iio = Hpdp + Hqdq = dH 

which follows because dq -\ dq = 1 = dp -\ dp and dq -\ dp = = dp -\ dq. This calculation proves 

Theorem 3.1 (Hamiltonian vector field) The Hamiltonian vector field Xh = {■ , H} satisfies 

Xh — I w = dH with u = dq h dp (3.1) 

Relation (3.1) may be taken as the definition of a Hamiltonian vector field. 



As a consequence of this formula, the flow of Xh preserves the closed exact two form ui for 
any Hamiltonian H. This preservation may be verified by a formal calculation using (3.1). Along 
{dq/dt,dp/dt) = {q,p) = {Hp,—Hq), we have 

^ = dq A dp + dq A dp = dHp A dp — dq A dHq 
= d{Hp dp + Hq dq) = d{XH -\u) = d{dH) = 

The first step uses the chain rule for differential forms and the third and last steps use the property 

of the exterior derivative d that = for continuous forms. The latter is due to equality of cross 
derivatives Hpq = Hqp and antisymmetry of the wedge product: dq A dp = —dp A dq. 

Consequently, the relation d{XH -I a;) = (PH = for Hamiltonian vector fields shows 
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Theorem 3.2 (Poincare's theorem for one degree of freedom) 

The flow of a Hamiltonian vector field is symplectic, which means it preserves the phase-space area, 
or two -form, uj = dq f\ dp. 

Definition 3.3 (Cartan's formula for the Lie derivative) The operation of Lie derivative of 
a differential form uj by a vector field Xh is defined by 

£xh^ ■= d{XH -loo) + Xh -idoj (3.2) 

Corollary 3.4 Because duj = 0, the symplectic property duj/dt = d{XH -Ilo) = in Poincare's 
Theorem 3.2 may be rewritten using Lie derivative notation as 

= ^ = £xhI^ := d{XH -Juj) + Xh -iduj =: {(HvXh) oo . (3.3) 
The last equality defines the divergence of the vector field Xh in terms of the Lie derivative. 
Remark 3.5 

• Relation (3.3) associates Hamiltonian dynamics with the symplectic flow in phase space of the 
Hamiltonian vector field Xh, which is divergenceless with respect to the symplectic form uj. 

• The Lie derivative operation defined in (3.3) is equivalent to the time derivative along the charac- 
teristic paths (flow) of the first order linear partial differential operator Xh, which are obtained 
from its characteristic equations in (2.1). This is the dynamical meaning of the Lie derivative 
£xh (3.2) for which invariance £xh^ — gives the geometric definition of symplectic flows 
in phase space. 

Theorem 3.6 (Poincare's theorem for N degrees of freedom) 

For a system of N particles, or N degrees of freedom, the flow of a Hamiltonian vector field preserves 
each subvolume in the phase space r*M^. That is, let uJn = dqn A dpn be the symplectic form expressed 
in terms of the position and momentum of the n~th particle. Then 



dt 



0, for UM = ^n=l^n, < N . 



The proof of the preservation of these Poincare invariants ujm with M = 1, 2, . . . , iV follows the 
same pattern as the verification above for a single degree of freedom. Basically, this is because each 
factor UJn = dqn A dpn in the wedge product of symplectic forms is preserved by its corresponding 
Hamiltonian flow in the sum 

M „ „ M M 

^H = Y. i^nTT- +PnW-) = ^ i^Pn^n " =^Xh„ = {-,H} 

n=l ^P"" n=l n=l 

That is, £xh^ujm vanishes for each term in the sum £xh'^m = 'I2n=i ■^Xh^^m since dq^ -idqn = 
Smn = dp^ -I dpn and dq^ _Jdpn = = dp^ _l dqn. 
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4 Fermat's theorem in geometrical ray optics 

4.1 Fermat's principle: Rays take paths of least optical length 

In geometrical optics, the ray path is determined by Fermat's principle of least optical length, 

6jni.,y,z)ds = 0. 

Here n{x, y, z) is the index of refraction at the spatial point {x, y, z) and ds is the element of arc length 
along the ray path through that point. Choosing coordinates so that the z— axis coincides with the 
optical axis (the general direction of propagation) , gives 

ds = [{dxf + {dyf + (dz) = [1 + + yV dz , 

with X = dx/dz and y = dy/dz. Thus, Fermat's principle can be written in Lagrangian form, with z 
playing the role of time, 

sjLix,y,x,y,z)dz = 0. 

Here, the optical Lagrangian is, 

L{x, y, i, y, z) = n{x, y, z)[l + + yY''^ =■ "/t > 
or, equivalently, in two-dimensional vector notation with q = (x, y), 

L(q,q,z) = n(q,z)[l + |q|Y''^ =: "/7 with 7 = [1 + |q| Y^^^ < 1 • 

Consequently, the vector Euler-Lagrange equation of the light rays is 

d / (iq\ _ d / (iq\ _ dn 
ds \ ds) ^ dz \^ dz ) 5q 



The momentum p canonically conjugate to the ray path position q in an "image plane", or on an 
"image screen" , at a fixed value of z is given by 

dL 

p = — r = ra7q which satisfies IpP = n^(l — 7^) . 

aq 

This implies the velocity q = p/(n^ — |p|^)^^^ • 

Hence, the momentum is real-valued and the Lagrangian is hyperregular, provided — |p|^ > 0. 
When = |pp, the ray trajectory is vertical and has grazing incidence with the image screen. 

Defining sin^ = dz/ds = 7 leads to |p| = ncosO, and gives the following geometrical picture of 
the ray path. Along the optical axis (the z— axis) each image plane normal to the axis is pierced at 
a point q = (x, y) by a vector of magnitude n(q, z) tangent to the ray path. This vector makes an 
angle 9 to the plane. The projection of this vector onto the image plane is the canonical momentum 
p. This picture of the ray paths captures all but the rays of grazing incidence to the image planes. 
Such grazing rays are ignored in what follows. 

Passing now via the usual Legendre transformation from the Lagrangian to the Hamiltonian de- 
scription gives 

H = p ■ q — L = n7|q|^ — n/7 = —nj = — [n(q, 2;)^ — |pp] 
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Thus, in the geometrical picture, the component of the tangent vector of the ray-path along the optical 
axis is (minus) the Hamiltonian, i.e. n{q,z)sin9 = —H. 

The phase space description of the ray path now follows from Hamilton's equations, 



dH -1 . dH -1 dn 

q = ^ = p = - 



2 



dp H^' ^ 5q 2H 

Remark 4.1 (Translation invariant media) If n = n(q), so that the medium is translation in- 
variant along the optical axis, z, then H = —n sin 9 is conserved. ( Conservation of H at an interface 
is Snell's law.) For translation-invariant media, the vector ray-path equation simplifies to 

1 duP' 

q = — r — — , Newtonian dynamics for q G . 

2ij^ aq 

Thus, in this case geometrical ray tracing reduces to "Newtonian dynamics" in z, with potential — n^(q) 
and with "time" reseated along each path by the value of \/2H determined from the initial conditions 
for each ray. 

4.2 Axisymmetric, translation invariant materials 

In axisymmetric, translation invariant media, the index of refraction is a function of the radius alone. 
Axisymmetry implies an additional constant of motion and, hence, reduction of the Hamiltonian 
system for the light rays to phase plane analysis. For such media, the index of refraction satisfies 

n(q, z) = n(r) , = |q| • 

Passing to polar coordinates (r, 0) with q = {x,y) = r(cos^, sin0) leads in the usual way to 



I |2 2 I 

IpI =Pt + 



Consequently, the optical Hamiltonian, 



H = -[n{rf-pl-pl/ 



21 1/2 



r 

is independent of the azimuthal angle (p; so its canonically conjugate "angular momentum" p^j, is 
conserved. 

Using the relation q • p = rpr leads to an interpretation of in terms of the image-screen phase 
space variables p and q. Namely, 

|pxq|2 = |p|2|q|2_(p.q)2=p2 

The conserved quantity p,/, = p x q = ypx — xpy is called the skewness function, or the Petzval invariant 

for axisymmetric media. Vanishing of p^ occurs for meridional rays, for which p and q are coUinear in 
the image plane. On the other hand, p^f, takes its maximum value for sagittal rays, for which p • q = 0, 
so that p and q are orthogonal in the image plane. 

Exercise 4.2 (Axisymmetric, translation invariant materials) Write Hamilton's canonical equa- 
tions for axisymmetric, translation invariant media. Solve these equations for the case of an optical 
fiber with radially graded index of refraction in the following form: 

n^(r) = -\- {n — i/r^)^ , X, u = constants, 

by reducing the problem to phase plane analysis. How does the phase space portrait differ between 
Pcf, = and ^ 0? Show that for p^p ^ the problem reduces to a Duffing oscillator in a rotating 
frame, up to a resettling of time by the value of the Hamiltonian on each ray "orbit. " 



D. D. Holm 



Imperial College London 



Applications of Poisson Geometry 



18 



4.3 The Petzval invariant and its Poisson bracket relations 

The skewness function 

S = P4, = p X q = yp.^ - xpy 

generates rotations of phase space, of q and p jointly, each in its plane, around the optical axis. Its 
square, (called the Petzval invariant) is conserved for ray optics in axisymmetric media. That is, 
{S^, H} = for optical Hamiltonians of the form, 

H=-[n{\ci\r-\p\'Y^'. 
We define the axisymmetric invariant coordinates by the map r*M^ i— > (q, p) <—>■ {X, Y, Z) , 

X = \qf>0, Y=\pf>0, Z = p-q. 
The following Poisson bracket relations hold 

{S^X} = 0, {S^Y} = 0, {S^Z} = 0, 

since rotations preserve dot products. In terms of these invariant coordinates, the Petzval invariant 
and optical Hamiltonian satisfy 

S"^ = XY - Z"^ >0, and H'^ = n'^{X) - Y > . 

The level sets of S*^ arc hyperboloids of revolution around the X = Y axis, extending up through the 
interior of the 5 = cone, and lying between the X— and Y— axes. The level sets of depend on 
the functional form of the index of refraction, but they are ^—independent. 

4.4 Poisson bracket for ray optics 

The Poisson brackets among the axisymmetric variables X, Y and Z close among themselves, 

{X,Y} = AZ, {Y,Z} = -2Y, {Z,X} = -2X. 

These Poisson brackets derive from a single Poisson bracket for X = {X, Y, Z) given by 

{F, H} = -VS^ VFxVH 

Consequently, we may re-express the equations of Hamiltonian ray optics in axisymmetric media with 
H = H{X, Y) as 

X = V^^ X VH . 

with Casimir S'^, for which {S"^,!!} = 0, for every H. Thus, the flow preserves volume (divX = 0) 
and the evolution takes place on intersections of level surfaces of the axisymmetric media invariants 
52 and H{X,Y). 

4.5 Recognition of the Lie-Poisson bracket for ray optics 

The Casimir invariant 5^ = XY — Z"^ is quadratic. In such cases, one may write the Poisson 
bracket in the suggestive form, 

, dF dH 
{F,H}--C,^X, — — 
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In this particular case, = 4, = 2 and = 2 and the rest either vanish, or are obtained from 
antisymmetry of C^j under exchange of any pair of its indices. These values are the structure constants 
of any of the Lie algebras sp(2,M), so(2,l), 5^(1, 1), or sZ(2,M). Thus, the reduced description of 
Hamiltonian ray optics in terms of axisymmetric variables is said to be "Lie-Poisson" on the dual 
space of any of these Lie algebras, say, sp{2, M)* for definiteness. We will have more to say about 
Lie-Poisson brackets later, when we reach the Euler-Poincare reduction theorem. 

Exercise 4.3 Consider the Poisson bracket 

{f,h} = - Vc-Vf xVh (4.1) 

Let c = x^ • Cx be a quadratic form on R^, and let C be the associated symmetric 3x3 matrix. Show 
that this is the Lie-Poisson bracket for the Lie algebra structure 

[u, v]c = C(u X v) 

What is the underlying matrix Lie algebra? What are the coadjoint orbits of this Lie algebra? 

Remark 4.4 (Coadjoint orbits) As one might expect, the coadjoint orbits of the group SP{2,M) 

are the hyperboloids corresponding to the level sets of S^. 

Remark 4.5 As we shall see later, the map r*M^ 1— > sp{2, M)* taking (q, p) 1— > (X, Y, Z) is an example 
of a momentum map. 

5 Geometrical Structure of Classical Mechanics 
5.1 Manifolds 

Configuration space: coordinates q G M, where M is a smooth manifold. 

o (/>~^ is a smooth change of variables. 
For later, smooth coordinate transformations: q ^ Q with dQ = "^dq 

Definition 5.1 A smooth manifold M is a set of points together with a finite (or perhaps countable) 
set of subsets Ua C M and 1-to-l mappings (pa '■ Ua ^ such that 

L []^Ua = M 

2. For every nonempty intersection Ua n Up, the set (pa {Ua n Up) is an open subset ofW^ and the 
1-to-l mapping (pp o (^~^ is a smooth function on cpa {Ua H Up) . 

Remark 5.2 The sets Ua in the definition are called coordinate charts. The mappings (pa o-i"^ called 
coordinate functions or local coordinates. A collection of charts satisfying 1 and 2 is called an 
atlas. Condition 3 allows the definition of manifold to be made independently of a choice of atlas. A 
set of charts satisfying 1 and 2 can always be extended to a maximal set; so, in practice, conditions 1 
and 2 define the manifold. 

Example 5.3 Manifolds often arise as intersections of zero level sets 

M={x|/i(x) = 0, i = l,...,k}, 
for a given set of functions fi : M" — M, i = 1, . . . , A;. 

If the gradients V fi are linearly independent, or more generally if the rank of{'Vf{x)} is a constant 
r for all x, then M is a smooth manifold of dimension n — r. The proof uses the Implicit Function 
Theorem to show that an {n — r) — dimensional coordinate chart may be defined in a neighborhood of 
each point on M. In this situation, the set M is called a submanifold o/M" (see [Le2003]). 
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Definition 5.4 If r = k, then the map {fi} is called a submersion. 
Exercise 5.5 Prove that all submersions are submanifolds (see [Le2003]). 
Definition 5.6 (Tangent space to level sets) Let 

M={x|/i(x) = 0, i = l,...,fc} 
he a manifold in M". The tangent space at each x G M, is defined by 

Note: we use the summation convention, i.e. repeated indices are summed over their range. 
Remark 5.7 The tangent space is a linear vector space. 
Example 5.8 (Tangent space to the sphere in M^) 

Example 5.9 (Tangent space to the sphere in M^) The sphere 5^ is the set of points {x,y,z) G 
solving x'^ + y"^ + = 1. The tangent space to the sphere at such a point {x,y,z) is the plane 
containing vectors {u, v, w) satisfying xu + yv + zw = 0. 

Definition 5.10 [Tangent bundle] The tangent bundle of a manifold M, denoted by TM, is the 
smooth manifold whose underlying set is the disjoint union of the tangent spaces to M at the points 
X G M; that is, 

TM = (J T^M 

Thus, a point of TM is a vector v which is tangent to M at some point x G M. 

Example 5.11 (Tangent bundle TS'^ of S^) The tangent bundle TS^ of S'^ G is the union of 

the tangent spaces of S'^: 

TS^ = { (x, y, z] u, v,w) ^M.^ | + + = 1 and xu + yv + zw = 0^ . 

Remark 5.12 (Dimension of tangent bundle TS"^) Defining TS'^ requires two independent con- 
ditions in R^; so dimTS'^ = 4. 

Exercise 5.13 Define the sphere 5""^ in M". What is the dimension of its tangent space TS^~^ ? 

Example 5.14 (The two stereographic projections of S"^ M^) 

The unit sphere 

S^ = {{x,y,z) : x^ + y^ + z^ = l} 
is a smooth two-dimensional manifold realized as a submersion in M^. Let 

Un = S^\{0, 0,1}, and Us = S^\{0, 0,-1} 

be the subsets obtained by deleting the North and South poles of S'^, respectively. Let 

Xn-Un^ (Cat, m) e , and xs ■ Us ^ {^s, Vs) e 
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be stereographic projections from the North and South poles onto the equatorial plane, z = 0. 

Thus, one may place two different coordinate patches in S"^ intersecting everywhere except at the points 

along the z—axis at z = 1 (North pole) and z = —1 (South pole). 

In the equatorial plane z = 0, one may define two sets of (right-handed) coordinates, 

K-.Ua^ m2\{0} , a = N,S, 

obtained by the following two stereographic projections from the North and South poles: 

(1) ( valid everywhere except z = 1) 

4>N{x,y,z) = {^N, vn) = 

(2) (valid everywhere except z = —1) 

Mx,y,^) = i^s,vs)=(^j^^, ^) . 

(The two complex planes are identified differently with the plane z = 0. An orientation-reversal is 
necessary to maintain consistent coordinates on the sphere.) 

One may check directly that on the overlap C/jv fl J/g the map, 

4>N o : M2\|o} ^ m2^{0} 

is a smooth diffeomorphism, given by the inversion 

Exercise 5.15 Construct the mapping from {Cn,i1n) — *■ {^s, Vs) o-'nd verify that it is a diffeomorphism 
in M2\{0}. Hint: {I + z){l - z) = 1 - z^ = + y"^ . 

Answer 5.16 

I - z 1 
{is, - VS) = T— Hn, vn) = ^2,2 (^^' • 

^ 'T z s jv "I" Vn 

The map {^NjVn) — ^ i^S^Vs) is smooth and invertible except at {^NjVn) = (0,0). 

Example 5.17 // we start with two identical circles in the xz-plane, of radius r and centered at 
x = ±2r, then rotate them round the z axis in M^, we get a torus, written T^. It's a manifold. 

Exercise 5.18 // we begin with a figure eight in the xz-plane, along the x axis and centered at the 

origin, and spin it round the z axis in M^, we get a "pinched surface" that looks like a sphere that has 
been "pinched" so that the north and south poles touch. Is this a manifold? Prove it. 

Answer 5.19 The origin has a neighbourhood diffeomorphic to a double cone. This is not diffeomor- 
phic to M^. A proof of this is that, if the origin of the cone is removed, two components remain; while 
if the origin of is removed, only one component remains. 

Remarlt 5.20 The sphere will appear in several examples as a reduced space in which motion takes 
place after applying a symmetry. Reduction by symmetry is associated with a classical topic in celestial 
mechanics known as normal form theory. Reduction may be "singular, " in which case it leads to 
"pointed" spaces that are smooth manifolds except at one or more points. For example different 
resonances of coupled oscillators correspond to the following reduced spaces: 1:1 resonance - sphere; 
1:2 resonance - pinched sphere with one cone point; 1:3 resonance - pinched sphere with one cusp 
point; 2:3 resonance - pinched sphere with one cone point and one cusp point. 
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5.2 Motion: Tangent Vectors and Flows 

Envisioning our later considerations of dynamical systems, we shall consider motion along curves c(t) 
parameterised by time t on a smooth manifold M. Suppose these curves are trajectories of a flow cpt of 
a vector field. We anticipate this means (l)t (c(0)) = c{t) and (f)to (f)g = <^f^g (flow property). The flow 
will be tangent to M along the curve. To deal with such flows, we will need the concept of tangent 
vectors. 

Recall from Definition 5.10 that the tangent bundle of M is 

TM = [j T^M. 

xeM 

We will now add a bit more to that definition. The tangent bundle is an example of a more general 
structure than a manifold. 

Definition 5.21 (Bundle) A bundle consists of a manifold B, another manifold M called the "base 
space" and a projection between them H : B ^ M. Locally, in small enough regions of x the inverse 
images of the projection 11 exist. These are called the fibers of the bundle. Thus, subsets of the 
bundle B locally have the structure of a Cartesian product. An example is {B, M, IT) consisting of 
(R2,R\n : R2 -> Ri). In this case, n : (x, y) G — > X G M^. Likewise, the tangent bundle consists 
of M, TM and a map tm ■ TM ^ M. 

Let X = (x-*^ , . . . , x") be local coordinates on M, and let v = (y^,...,v'^) be components of a 
tangent vector. 

T^M = G M" I ^ - v = 0,1 = 1,..., 

for 

M = (x e I fi{x) = 0,i = l,...,m} 

These 2n numbers {x,v) give local coordinates on TM, whose dimension is dim TM = 2dimM. The 
tangent bundle projection is a map tm '■ TAI M which takes a tangent vector v to a point 
x € M where the tangent vector v is attached (that is, v G TxM). The inverse of this projection 
T^{x) is called the fiber over x in the tangent bundle. 

Vector fields, integral curves and fiows 

Definition 5.22 A vector field on a manifold M is a map X : M ^ TM that assigns a vector X{x) 
at each point x G M. This implies that tm o X = Id. 

Definition 5.23 An integral curve of X with initial conditions xq at t = Q is a differentiable map 
c :]a, b[-^ M, where ]a, b[ is an open interval containing 0, such that c(0) = and c'{t) = X {c{t)) for 
all t g]o, b[. 

Remark 5.24 A standard result from the theory of ordinary differential equations states that X being 
Lipschitz implies its integral curves are unique and C^ [CoLel984j. The integral curves c{t) are 
differentiable for smooth X. 
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5.3 Summary 

Definition 5.25 The flow of X is the collection of maps (pt '■ M ^ M , where t — > 4)t{x) is the 
integral curve of X with initial condition x. 

Remark 5.26 

1. Existence and uniqueness results for solutions of c'{t) = X{c{t)) guarantee that flow (j) of X is 
smooth in {x, t), for smooth X. 

2. Uniqueness implies the flow property 

4>t+s = ct>to<l>s, (FP) 

for initial condition (f)o = Id. 

3. The flow property (FP) generalizes to the nonlinear case the familiar linear situation where 
M is a vector space, X(x) = Ax is a linear vector field for a hounded linear operator A, and 
4>tix) = e^*x. 

Differentials of functions and the cotangent bundle 

We are now ready to define differentials of smootli functions and tlie cotangent bundle. 

Let / : M ^ R be a smootli function. Wc differentiate / at x G M to obtain T^f : T^M 7)(^.)IR. 
As is standard, we identify Tj^^:)!^ with M itself, thereby obtaining a linear map df{x) : TxM M. 
The result df{x) is an element of the cotangent space T*M, the dual space of the tangent space T^M. 
The natural pairing between elements of the tangent space and the cotangent space is denoted as 
(•,•): T*M X T^M ^ R. 

In coordinates, the linear map df{x) : T^M — > M may be written as the directional derivative, 

{df{x),v) = df{x)-v=^^-v\ 

for all V G TxM. (Reminder: the summation convention is intended over repeated indices.) Hence, 
elements df{x) G T*M are dual to vectors v G TxM with respect to the pairing (• , •). 

Definition 5.27 df is the differential of the function f. 

Definition 5.28 The dual space of the tangent bundle TM is the cotangent bundle T*M. That is, 

{TxM)* = T*M and T*M = [J T*M . 

X 

Thus, replacing v G TxM with df G T*M, for all a; G M and for all smooth functions f : M —>■ W, 
yields the cotangent bundle T*M. 
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Differential bases When the basis of vector fields is denoted as ^ for i = 1, . . . , n, its dual basis 
is often denoted as dx*. In this notation, the differential of a function at a point x G M is expressed 
as 

df(x) = ^dx' 
•'^ ^ dx' 

The corresponding pairing (• , •) of bases is written in this notation as 

Here Sj is the Kronecker delta, which equals unity for i = j and vanishes otherwise. That is, defining 
T*M requires a pairing (-,■): T*M x TM R. 

(Different pairings exist for curvilinear coordinates, Riemannian manifolds, etc.) 

6 Derivatives of different iable maps — the tangent lift 

We next define derivatives of differentiable maps between manifolds (tangent lifts). 

We expect that a smooth map f : U ^ V from a chart U C M to a chart V C N, will lift to a 
map between the tangent bundles TM and TN so as to make sense from the viewpoint of ordinary- 
calculus, 

U X R™ c TM — > y X R" C TAf 
{q\...,r;X\...,X"')^{Q\...,Q-;Y\...,Y^) 

Namely, the relations between the vector field components should be obtained from the diff'crcntial of 
the map f : U —>■ V. Perhaps not unexpectedly, these vector field components will be related by 

in which the quantity called the tangent lift 

of the function / arises from the chain rule and is equal to the Jacobian for the transformation 
Tf : TM ^ TN. 

The dual of the tangent lift is the cotangent lift, explained later in section 10.1. Roughly speaking, 
the cotangent lift of the function /, 

T*f = ^ 
^ dQ 

arises from 

PidQ"- = ajdq\ so f3i = aj—^ 
and T*/ : T*N ^ T*M. Note the directions of these maps: 

Tf : q,X eTM^Q,Y eTN 

/: qeM^QeN 

T*f : Q,j3&T*N<-^q,aeT*M (map goes the other way, see the picture 
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6.1 Summary remarks about derivatives on manifolds 

Definition 6.1 (DifFerentiable map) A map f : M ^ N from manifold M to manifold N is said 
to be differentiable (resp. C^) if it is represented in local coordinates on M and N by differentiable 
(resp. C^) functions. 

Definition 6.2 (Derivative of a differentiable map) The derivative of a differentiable map 

f : M 

at a point x E M is defined to be the linear map 

T^f : T^M ^ T^N 

constructed, as follows. For v G T^M , choose a curve c{t) that maps an open interval t G (— e, e) 
around the point t = to the manifold M 

c : (-e,e) — > M 
with c(0) = X 

dc 



and velocity vector c'(0) . 



= V . 

t=o 



Then T^f ■ v is the velocity vector at t = of the curve foe: M — >■ A^. That is, 



d 



T^f .v= —ficit)) =^—c(t) 



dt 



df d 



t=o 



Definition 6.3 The union Tf = Txf of the derivatives Txf : TxM — ^ TxN over points x E M is 
called the tangent lift of the map f : M —>■ N. 

Remark 6.4 The chain-rule definition of the derivative Txf of a differentiable map at a point x 

depends on the function f and the vector v. Other degrees of differentiability are possible. For example, 
if M and N are manifolds and f : M ^ N is of class C'^'^^, then the tangent lift (Jacobian) Txf : 
TxM TxN is C^. 

Exercise 6.5 Let (f)t : S"^ ^ S"^ rotate points on S"^ about a fixed axis through an angle ip{t). Show 
that (pt is the flow of a certain vector field on 5^. 

Exercise 6.6 Let f : S"^ ^ W be defined by f{x,y,z) = z. Compute df using spherical coordinates 

{e,c(>). 

Exercise 6.7 Compute the tangent lifts for the two stereographic projections of S'^ — ^ in example 
5.14- That is, assuming {x,y,z) depend smoothly on t, find 

L How {^n,Vn) depend on {x,y, z). Likewise, for {^s-i'Hs)- 

2. How {iN,f]N) depend on 

Hint: Recall (1 + z){l — z) = 1 — z"^ = x'^ + y'^ and use xx + yy + zz = when {x, y, z) is tangent to 
5^ at {x, y, z). 
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7 Lie groups and Lie algebras 
7.1 Matrix Lie groups 

Definition 7.1 A group is a set of elements with 

1. A binary product (multiplication), G x G ^ G 

The product of g and h is written gh 
The product is associative, {gh)k = g{hk) 

2. Identity element e : eg = g and ge = g, \/g G G 

3. Inverse operation G ^ G, so that gg~^ = g~^g = e 

Definition 7.2 A Lie group is a smooth manifold G which is a group and for which the group 
operations of multiplication, {g, h) — ^ gh for g,h E G, and inversion, g — g~^ with gg~^ = g''^g = e, 
are smooth. 

Definition 7.3 A matrix Lie group is a set of invertihle n x n matrices which is closed under 
matrix multiplication and which is a suhmanifold ofW^^"^. The conditions showing that a matrix Lie 
group is a Lie group are easily checked: 

A matrix Lie group is a manifold, because it is a submanifold ofW^^'^ 

Its group operations are smooth, since they are algebraic operations on the matrix entries. 

Example 7.4 (The general linesir group GL(n,M)) The matrix Lie group GL{n,M) is the group 
of linear isomorphisms o/M" to itself. The dimension of the matrices in GL{n,W) is n^. 

Proposition 7.5 Let K G GL(n,M) be a symmetric matrix, = K. Then the subgroup S of 
GL(n,M) defined by the mapping 

S = {U e GL{n,R)\ U^KU = K} 
is a submanifold o/M"^" of dimension n{n — l)/2. 

Remark 7.6 The subgroup S leaves invariant a certain symmetric quadratic form under linear trans- 
formations, 5 X M" — > M" given by x ^ Ux, since 

So the matrices U E S change the basis for this quadratic form, but they leave its value unchanged. 
Thus, S is the isotropy subgroup of the quadratic form associated with K. 

Proof. 

Is 5 a subgroup? We check the following three defining properties 

1. Identity: 

I e 5 because I^KI = K 

2. Inverse: 

U e S =^ U~'^ G S, because 
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3. Closed under multiplication: 
U,V e S ^UV e S, because 
{UV)'^KUV = V'^{U^KU)V = V'^{K)V = K 

Hence, 5 is a subgroup of GL(n,M). 

Is S* is a submanifold of M"^" of dimension n(n — l)/2? 

Indeed, S is the zero locus of the mapping U KU^ — K. This makes it a submanifold, because 
it turns out to be a submersion. 

For a submersion, the dimension of the level set is the dimension of the domain minus the 
dimension of the range space. In this case, this dimension is "n? — n{n + l)/2 = n{n — l)/2 

Exercise 7.7 Explain why one can conclude that the zero locus map for S is a submersion. In 
particular, pay close attention to establishing the constant rank condition for the linearization of this 
map. 

Solution Here is why 5 is a submanifold of i?"^". 

First, S is the zero locus of the mapping 

U U^KU - K , (locus map) 

Let U E S, and let 5U be an arbitrary element of i?"^". Then linearize to find 

{U + 5UfK{U + 5U)-K = U^KU -K + SU^KU + U^K5U + 0{5Uf, . 

We may conclude that 5 is a submanifold of i?"^" if we can show that the linearization of the locus 
map, namely the linear mapping defined by 

L = 5U ^ SU^KU + U^KSU , i?"""" ^ i?"''" 

has constant rank for all U & S. 

Lemma 7.8 The linearization map L is onto the space ofnxn of symmetric matrices and hence the 
original map is a submersion. 

Proof that L is onto. 

• Both the original locus map and the image of L lie in the subspace of n x n symmetric matrices. 

• Indeed, given U and any symmetric matrix S we can find SU such that 

SU^KU + U^KSU = S . 

Namely 

6U = K-^U-'^S/2. 

• Thus, the linearization map L is onto the space of n x n of symmetric matrices and the original 
locus map U ^ U KU^ — K to the space of symmetric matrices is a submersion. 

For a submersion, the dimension of the level set is the dimension of the domain minus the dimension 
of the range space. In this case, this dimension is — n(n + l)/2 = n(n — l)/2. 
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Corollary 7.9 (5 is a matrix Lie group) S is both a subgroup and a submanifold of the general 
linear group GL{n,W). Thus, by definition 7.3, S is a matrix Lie group. 

Exercise 7.10 What is the tangent space to S at the identity, TjS? 

Exercise 7.11 Show that for any pair of matrices A,BE: TjS, the matrix commutator [A,B] = 
AB-BAe TjS. 

Proposition 7.12 The linear space of matrices A satisfying 

A^K + KA = Q 

defines TjS, the tangent space at the identity of the matrix Lie group S defined in Proposition 7.5. 

Proof. Near the identity the defining condition for S expands to 

(7 + eA^ + 0{e^))K{I + eA + 0{e^)) = K, for e < 1 . 
At Hnear order 0(e) one finds, 

A'^K + KA = 0. 

This relation defines the linear space of matrices A G TjS. 
lfA,Be TiS, does it follow that [A, B] € TjSl 

Using [^,5]-^ = [B'^,A'^], we check closure by a direct computation, 

[B'^ , A^]K + K[A, B] = B'^A^K-A^B'^K + KAB-KBA 

= B^A^K - A^B'^K - A^KB + B'^KA = . 

Hence, the tangent space of S at the identity TjS is closed under the matrix commutator [• , •]. 

Remark 7.13 In a moment, we will show that the matrix commutator for TjS also satisfies the Jacobi 
identity. This will imply that the condition A^K + KA = defines a matrix Lie algebra. 

7.2 Defining Matrix Lie Algebras 

We are ready to prove the following, in preparation for defining matrix Lie algebras. 

Proposition 7.14 Let S be a matrix Lie group, and let A,BE: TjS (the tangent space to S at the 
identity element). Then AB - BA e T/S. 

The proof makes use of a lemma. 

Lemma 7.15 Let R be an arbitrary element of a matrix Lie group S, and let B G TjS. Then 
RBR-^ G TiS. 

Proof of lemma. Let Rsit) be a curve in S such that Rb{0) = I and R'{0) = B. Define S{t) = 
RRB{t)R-^ e TjS for all t. Then S{0) = I and S'{0) = RBR'^. Hence, S'{0) € TjS, thereby proving 
the lemma. 
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Proof of Proposition 7.14. Let Ra{s) be a curve in S such that Ra{0) = I and R'j^{0) = A. 
Define S{t) = RA{t)BRA{t)-'^ G TiS. Then the lemma imphes that S{t) G TjS for every t. Hence, 
S'{t) e TjS, and in particular, S"(0) = AB-BAe TjS . 

Definition 7.16 (Matrix commutator) For any pair of n x n matrices A,B, the matrix com- 
mutator is defined as [A, B] = AB — BA. 

Proposition 7.17 (Properties of the matrix commutator) The matrix commutator has the fol- 
lowing two properties: 

(i) Any two n x n matrices A and B satisfy 

[B,A] = -[A,B] 

(This is the property of skew- symmetry.) 

(a) Any three n x n matrices A, B and C satisfy 

[[A,B],C]-h[[B,C],A]-^[[C,A],B]=0 

(This is known as the Jacobi identity.) 

Definition 7.18 (Matrix Lie algebra) A matrix Lie algebra q is a set of nxn matrices which is a 
vector space with respect to the usual operations of matrix addition and multiplication by real numbers 
(scalars) and which is closed under the matrix commutator [• , •]. 

Proposition 7.19 For any matrix Lie group S, the tangent space at the identity TjS is a matrix Lie 
algebra. 

Proof. This follows by proposition 7.14 and because TjS is a vector space. 
7.3 Examples of matrix Lie groups 

Example 7.20 (The Orthogonal Group 0{n)) The mapping condition U^KU = K in Proposi- 
tion 7.5 specializes for K = I to U'^U = I, which defines the orthogonal group. Thus, in this case, S 
specializes to 0{n), the group ofn x n orthogonal matrices. The orthogonal group is of special interest 
in mechanics. 

Corollary 7.21 (0(n) is a matrix Lie group) By Proposition 7.5 the orthogonal group 0{n) is 
both a subgroup and a submanifold of the general linear group GL(n,M). Thus, by definition 7.3, the 
orthogonal group 0{n) is a matrix Lie group. 

Example 7.22 (The Special Linear Group SL{n,R)) The subgroup ofGL{n,R) withdet{U) = 1 
is called SL{n,M.). 

Example 7.23 (The Special Orthogonal Group SO{n)) The special case of S with dei{U) = 1 
and K = I is called SO{n). In this case, the mapping condition U^KU = K specializes to U^U = I 
with the extra condition det(?7) = 1. 
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Example 7.24 (The tangent space of SO{n) at the identity TjSO{n)) The special case with K 
I ofTiSO{n) yields, 

A^ + A = 0. 

These are antisymmetric matrices. Lying in the tangent space at the identity of a matrix Lie group, 
this linear vector space forms a matrix Lie algebra. 

Example 7.25 (The Symplectic Group) Suppose n = 21 (that is, let n be even) and consider the 
nonsingular skew-symmetric matrix 

J- f ° ^ 
[-1 

where I is the I x / identity matrix. One may verify that 

Sp{l) = {U e GL{2l,R)\ U^JU = J} 

is a group. This is called the symplectic group. Reasoning as before, the matrix algebra TjSp{l) is 
defined as the set of n x n matrices A satisfying JA^ + AJ = 0. This algebra is denoted as sp{l). 

Example 7.26 (The Special Euclidean Group) Consider the Lie group o/4 x 4 matrices of the 
form 

'" R V 



E{R,v) = 



1 



where R G SO{3) and v & M.^. This is the special Euclidean group, denoted SE{3). The special 
Euclidean group is of central interest in mechanics since it describes the set of rigid motions and 
coordinate transformations of three-dimensional space. 

Exercise 7.27 A point P in undergoes a rigid motion associated with E{Ri,vi) followed by a rigid 
motion associated with E{R2,V2). What matrix element of SE{3) is associated with the composition 
of these motions in the given order? 

Exercise 7.28 Multiply the special Euclidean matrices of SE{3). Investigate their matrix commuta- 
tors in their tangent space at the identity. (This is an example of a semidirect product Lie group.) 

Exercise 7.29 (Tripos question) When does a stone at the equator of the Earth weigh the most? 
Two hints: (1) Assume the Earth's orbit is a circle around the Sun and ignore the declination of the 
Earth's axis of rotation. (2) This is an exercise in using SE(2). 

Exercise 7.30 Suppose the nx n matrices A and M satisfy 

AM + MA^ = . 

Show that exp{At)Mexp{A^t) = M for all t. Hint: A'^M = M{-A^Y. This direct calculation shows 
that for A G so{n) or A & sp{l), we have exp(74i) G SO{n) or eyi'g{At) G Sp{l), respectively. 



7.4 Lie group actions 

The action of a Lie group G on a manifold M is a group of transformations of M associated to elements 
of the group G, whose composition acting on M is corresponds to group multiplication in G. 

Definition 7.31 Let M be a manifold and let G be a Lie group. A left action of a Lie group G on 
M is a smooth mapping ^ : G x M ^ M such that 
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(i) <I>(e, x) = X for all x £ M, 

(a) ^{g, ^{h, x)) = ^{gh, x) for all g,h e G and x e M, and 
(in) ^{g, •) is a diffeomorphism on M for each g E G. 

We often use the convenient notation gx for ^{g,x) and think of the group element g acting on the 
point X G M. The associativity condition (ii) above then simply reads {gh)x = g{hx). 

Similarly, one can define a right action, which is a map ^ : M x G ^ M satisfying e) = x and 
'^{'^{x,g),h) = ^{x,gh). The convenient notation for right action is xg for ^{x,g), the right action 
of a group element g on the point x G M. Associativity h) = ^{x,gh) is then be expressed 

conveniently as {xg)h = x{gh). 

Example 7.32 (Properties of group actions) The action ^ : G x M ^ M of a group G on a 
manifold M is said to be: 

1. transitive, if for every x,y E M there exists a g E G, such that gx = y; 

2. free, if it has no fixed points, that is, ^g{x) = x implies g = e; and 

3. proper, if whenever a convergent subsequence in M exists, and the mapping gnXn converges 
in M, then {gn} has a convergent subsequence in G. 

Orbits. Given a group action of G on M, for a given point x G M, the subset 

Orb a; = {gx\g£ G} C M , 

is called the group orbit through x. In finite dimensions, it can be shown that group orbits are always 
smooth (possibly immersed) manifolds. Group orbits generalize the notion of orbits of a dynamical 
system. 

Exercise 7.33 The flow of a vector field on M can be thought of as an action o/M on M. Show that 
in this case the general notion of group orbit reduces to the familiar notion of orbit used in dynamical 
systems. 

Theorem 7.34 Orbits of proper group actions are embedded submanifolds. 

This theorem is stated in Ghapter 9 of [MaRal994], who refer to [AbMal978] for the proof. 

Example 7.35 (Orbits of 50(3)) A simple example of a group orbit is the action of SO {3) on 
given by matrix multiplication: The action of A E SO{3) on a point x G is simply the product Ax. 
In this case, the orbit of the origin is a single point (the origin itself), while the orbit of any other 
point is the sphere through that point. 

Example 7.36 (Orbits of a Lie group acting on itself) The action of a group G on itself from 
either the left, or the right, also produces group orbits. This action sets the stage for discussing the 
tangent lifted action of a Lie group on its tangent bundle. 

Left and right translations on the group are denoted Lg and Rg, respectively. For example, 
Lg : G ^ G is the map given by h ^ gh, while Rg : G ^ G is the map given by h ^ hg, for g,h E G. 
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(a) Left translation Lg : G ^ G; h ^ gh defines a transitive and free action of G on itself. Right 

multiplication Rg : G ^ G; h ^ hg defines a right action, while h — hg~^ defines a left action 
of G on itself. 

(b) G acts on G by conjugation, g ^ Ig = Rg-i o Lg. The map Ig : G ^ G given by h ^ 9hg~^ is the 

inner automorphism associated with g. Orbits of this action are called conjugacy classes. 

(c) Differentiating conjugation at e gives the adjoint action of G on q: 

Adg := TJg : TeG = ^ T^G = g. 

Explicitly, the adjoint action of G on g is given by 

Ad : G X g ^ , Adg{0 = Te{Rg-i o Lg)^ 

We have already seen an example of adjoint action for matrix Lie groups acting on matrix Lie 
algebras, when we defined S{t) = RA{t)BRA{t)~^ G TjS as a key step in the proof of Proposition 
7.14. 

(d) The coadjoint action of G on g* , the dual of the Lie algebra g ofG, is defined as follows. Let 
Ad* : 0* ^ 0* be the dual of Adg, defined by 

{Ad*ga,0 = {a,AdgO 
for a G 0*; ^ G and pairing {■,■): g* x g ^ M.. Then the map 

: G X 0* — 0* given by {g,a) Ad*_iQ; 
is the coadjoint action of G on g* . 

7.5 Examples: SO {3), SE{3), etc. 

A basis for the matrix Lie algebra so(3) and a map to 

The Lie algebra of SO{n) is called so{n). A basis (ei, 62, 63) for so(3) when n = 3 is given by 



X 



-z y 

z —X 
-y X 



xei + 2/62 + zez 



Exercise 7.37 Show that [61,62] = 63 and cyclic permutations, while all other matrix commutators 
among the basis elements vanish. 

Example 7.38 (The isomorphism between so(3) and M^) The previous equation may be written 
equivalently by defining the hat-operation {•) as 

Xjj = EijkX^ , where {x^,x'^,x^) = (x, y, z) . 

Here £123 = 1 and €213 = —1, with cyclic perm,utations. The totally antisymmetric tensor €ijk = 
~ ^jik = ~ ^ikj olso defines the cross product of vectors in M^. Consequently, we may write, 

(x X Y)i = €ijkx^y^ = 1tijy> , that is, x x y = xy 

Exercise 7.39 What is the analog of the hat map so(3) 1— > for the three dimensional Lie algebras 
sp{2,R), so{2,l), su{l,l), orsl{2,R)? 

Background reading for this lecture is Chapter 9 of [MaRal994]. 
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Compute the Adjoint and adjoint operations by differentiation 

1. Differentiate Ig{h) wrt h at h = e to produce the Adjoint operation 

Ad : Gxg^g : Adg 77 = Telg ij 

2. Differentiate Adg rj wrt g at g = e in the direction ^ to get the Lie bracket [^,r]] : QX ^ Q and 

thereby to produce the adjoint operation 

Te{Adgri)^=[^,ri] = ad^ri 

Compute the co- Adjoint and coadjoint operations by taking duals 

1. Ad* : g* ^ 0*, tlie dual of Ad^, is defined by 

(Ad^a, = {a, AdgO 
for a G 0*, ^ G and pairing {■,■): g* x g ^ R. The map 

: G X 0* ^ 0* given by {g, a) Ad*_iQ; 
defines the co- Adjoint action of G on 0* . 

2. The pairing 

(ad|a,r/) = {a,ad^r]) 
defines the coadjoint action of on 0*, for a e g* and ^, 77 G 0. 

See Chapter 9 of [MaRal994] for more discussion of the Ad and ad operations. 

Example: the rotation group 5*0(3) 
The Lie algebra so(3) and its dual. 

The special orthogonal group is defined by 

SO{3) := {A \ A a3x3 orthogonal matrix, det(^) = 1} . 
Its Lie algebra so(3) is formed by 3 x 3 skew symmetric matrices, and its dual is denoted so(3)*. 

The Lie algebra isomorphism " : (<so(3), [•, •]) (M^, x) 

The Lie algebra (so(3), [•,•]), where [•, ■] is the commutator bracket of matrices, is isomorphic to the 
Lie algebra (M"^, x), where x denotes the vector product in M^, by the isomorphism 














u := {u^,u'^,u^) G M:' 


u := 







—u 












G so(3) , that is, Ujj := — eij^u 



k 



Equivalently, this isomorphism is given by 

uv = u X V for all u, v G M^. 

The following formulas for u,v,w G W"' may be easily verified: 

(uxv)" = [u,v] 
[u,v]w = (u X v) X w 

u-v = — ^trace(uv). 
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Ad action of 50(3) on so(3) The corresponding adjoint action of S0{3) on so(3) may be obtained 
as follows. For SO{3) we have Ia{B) = ABA'^. Differentiating B{t) at B{0) = Id gives 



AdAV = — 
at 



AB(t)A-^ = AvA-^ , with v = S'(0) 

t=o 



One calculates the pairing with a vector w G as 

AdAv(vif) = Av{A~^w) = A{v x A~^w) = Av x w = {AvYw 
where we have used a relation 

^(u X v) = X 

which holds for any u,v G and A G 50(3). 
Consequently, 

AdAV = (^v)^ 

Identifying so(3) ~ then gives 

Ad^v = Av. 

So (speaking prose all our lives) the adjoint action of 50(3) on so(3) may be identitified with multi- 
plication of a matrix in 50(3) times a vector in M^. 

ad-action of so(3) on so(3) Differentiating again gives the ad-action of the Lie algebra so(3) on 
itself: 

(e*"vY = (uv)^ = (u X v)' . 

So in this isomorphism the vector cross product is identified with the matrix commutator of skew 
symmetric matrices. 



[u,v] = adov 



Infinitesimal generator Likewise, the infinitesimal generator corresponding to u G has the 

expression 



d 

UM3(x) := - 



e*"x = ux = u X X. 



Exercise 7.40 What is the analog of the hat map so(3) i— ^ for the three dimensional Lie algebras 
sp{2,R), so{2,l), s'u(l,l), orsl{2,R)? 

The dual Lie algebra isomorphism ~ : so(3)* 

Coadjoint actions The dual so(3)* is identified with R'^ by the isomorphism 
n G n G so(3)* : n(u) := n • u for any u G . 

In terms of this isomorphism, the co-Adjoint action of 50(3) on so(3)* is given by 

Ad^_i n = (^n)~ 

and the coadjoint action of so(3) on so(3)* is given by 

adtn=(nxu)~. (7.1) 
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Computing the co-Adjoint action of SO{3) on so(3)* This is given by 

(^Ad^_i n) (u) = n • Ad^-i u = n • (^~^u)" = n • a^u 
= An-u = {Any (u), 

tliat is, the co-Adjoint action of 50(3) on so(3)* has the expression 

Ad^^-i n = (^n)~ , 

Therefore, the co-Adjoint orbit O = {AU \ A e SO{3)} C of SO{3) through n G 
2-sphere of radius ||n||. 



IS a 



Computing the coadjoint action of so(3) on so(3)* Let u, v G and note that 



ad^n,v^ = ^n, [u,v]^ = ^n,(u X v)-^ = n • (u X v) 
= (n X u) • V = (n X u)~,v) , 

which shows that ad^ II = (11 x u)"", thereby proving (7.1). Therefore, TjjO = {ll x u | u G 

since the plane perpendicular to 11, that is, the tangent space to the sphere centered at the origin of 
radius ||n||, is given by (lI x u | u G 



8 Lifted Actions 

Definition 8.1 Let ^ : G x M ^ M be a left action, and write ^g{x) = ^{g,x) for x G M. The 
tangent lift action of G on the tangent bundle TM is defined by gv = Tx^g{v) for every v G T^M. 

Remark 8.2 In standard calculus notation, the expression for tangent lift may be written as 
d 

T^^ ■ V = —^(c(t)) =—c'(t) =:D^(x)-v, with c(0) = x , c'(0) = v. 
dt t=o oc t=o 

Definition 8.3 If X is a vector field on M and (p is a differentiable map from, AI to itself then the 
push-forward of X by (p is the vector field (p^X defined by {(p^X) {(p{x)) = Tx(p{X{x)) . That is, the 
following diagram commutes: 



TcP 



TM 



TM 



X 



.X 



M 



M 



If (p is a diffeomorphism then the pull-back (p*X is also defined: {(p*X) (x) = Tj^MCp ^ {X {(p{x))) . 



Definition 8.4 Let ^ : G x M ^ M be a left action, and write ^g{m) = <^{g,m). Then G has a left 
action on X E X{M) (the set of vector fields on M) by the push-forward: gX = {^g)^ X. 
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Definition 8.5 Let G act on M on the left. A vector field X on M is invariant with respect to this 
action (we often say "G-invariant" if the action is understood) if gX = X for all g E G; equivalently 
(using all of the above definitions!) g {X{x)) = X{gx) for all g E G and all x E X. 

Definition 8.6 Consider the left action of G on itself by left multiplication, ^g{h) = Lg{h) = gh. A 
vector field on G that is invariant with respect to this action is called left-invariant. From Definition 
8.5, we see that X is left-invariant if and only if g{X{h)) = X{gh), which in less compact notation 
means ThLgX{h) = X{gh). The set of all such vector fields is written X^{G). 

Proposition 8.7 Given a ^ G TgG, define X^{g) = g^ (recall: g^ = TgLg^). Then X^ is the unique 
left-invariant vector field such that X^{e) = ^. 

Proof To show that is left-invariant, we need to show that g (Xp{h)) = XHgh) for every 



g,h e G. This follows from the definition of X^ and the associativity property of group actions: 

g {X^{h)) = g {hO = (gh) ^ = X^{gh) 

We repeat the last line in less compact notation: 

ThLg {X^{h)) = ThLg (hC) = T^Lghi = X^igh) 

For uniqueness, suppose X is left-invariant and -^(e) = ^. Then for any g e G, we have X{g) = 
g{X{e)) = gC = X^{g). U 

Remark 8.8 Note that the map ^ i— X^ is an vector space isomorphism from T^G to X^{G). 

All of the above definitions have analogues for right actions. The definitions of right-invariant, 
X^{G) and X^ use the right action of G on itself defined by ^{g, h) = Rg{h) = kg. 

Exercise 8.9 There is a left action of G on itself defined by ^g{h) = hg~^. 

We will use the map ^ i-^ X^ to relate the Lie bracket on g, defined as r/] = ad^ry, with the 
Jacobi-Lie bracket on vector fields. 

Definition 8.10 The Jacobi-Lie bracket on X{M) is defined in local coordinates by 

[X,Y]j_L = {DX)-Y-{DY)-X 

which, in finite dimensions, is equivalent to 

[X, Y]j_L = -{X-V)Y + {Y- V)X = - [X, Y] 

Theorem 8.11 [Properties of the Jacobi-Lie bracket] 
L The Jacobi-Lie bracket satisfies 

d 



t=0 



where ^> is the flow of X. (This is coordinate-free, and can be used as an alternative definition.) 
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2. This bracket makes X^{M) a Lie algebra with [X,Y]j-i = —[X,Y], where [X,Y] is the Lie 
algebra bracket on X{M). 

3. (p^X, Y] = [(p^X, (p^Y] for any differentiable (p-.M^M. 

Remark 8.12 The first property of the Jacobi-Lie bracket is proved for matrices in section 9. The 
other two properties are proved below for the case that M is the Lie group G. 

Theorem 8.13 X^{G) is a subalgebra ofX{G). 

Proof. Let X,Y G X^{G). Using the last item of the previous theorem, and then the G invariance of 
X and Y, gives the push-forward relations 

{Lgl [X, Y]j_L = [{Lg), X, {Lg\ Y]j_L 
for all g & G. Hence [X,Y]j^l G X^{G). This is the second property in Theorem 8.11. ■ 

Theorem 8.14 Set [X^,X^]j_L{e) = [^,?7] for every ^, G 0, where the bracket on the right is the 
Jacobi-Lie bracket. (We say: the Lie bracket on g is the pull-back of the Jacobi-Lie bracket by the map 

Proof. The proof of this theorem for matrix Lie algebras is relatively easy: we have already seen that 
adyii? = AB — BA. On the other hand, since X^{C) = CA for all C, and this is linear in C, we have 
DX^{I) ■A = AB, so 

[A, B] = [Xi, j_i(7) = DX^il) . - DX^il) ■ X^{I) 

= DX^{I) ■ A - DX\{I) ■ B = AB - BA 

This is the third property of the Jacobi-Lie bracket listed in Theorem 8.11. For the general proof, see 
Marsden and Ratiu [MaRal994], Proposition 9.14. ■ 

Remark 8.15 This theorem, together with Item 2 in Theorem 8.11, proves that the Jacobi-Lie bracket 
makes g into a Lie algebra. 

Remark 8.16 By Theorem 8.13, the vector field [X^,X!^\ is left-invariant. Since [X^ ,X^]j-L{e) = 
{i,rj\, it follows that 



[X^ ,X^ \ — X| 



L 



Definition 8.17 Let ^ : G x M ^ M be a left action, and let ^ G g. Let g{t) be a path in G such 
that y(0) = e and g'{0) = ^. Then the infinitesimal generator of the action in the ^ direction is the 
vector field on M defined by 

d 



t=0 



Remark 8.18 Note: this definition does not depend on the choice of g{t). For example, the choice in 
Marsden and Ratiu [MaRal994j is exp(f0, where exp denotes the exponentiation on Lie groups (not 
defined here). 

Exercise 8.19 Consider the action of SO{3) on the unit sphere 5^ around the origin, and let ^ = 
(0,0,lj~. Sketch the vector field ^m- (Hint: the vectors all point "Eastward.") 
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Theorem 8.20 For any left action of G, the Jacobi-Lie bracket of infinitesimal generators is related 
to the Lie bracket on Q as follows (note the minus sign): 

[^m,Vm] = - [^,v]m 
For a proof, see Marsden and Ratiu [MaRal994], Proposition 9.3.6. 

Exercise 8.21 Express the statements and formulas of this lecture for the case of SO{3) action on 
its Lie algebra so(3). (Hint: look at the previous lecture.) Wherever possible, translate these formulas 
to by using the " map: so(3) — M^. 

Write the Lie algebra for so(3) using the Jacobi-Lie bracket in terms of linear vector fields on M^. 
What are the characteristic curves of these linear vector fields ? 
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9 Handout: The Lie Derivative and the Jacobi-Lie Bracket 

Let X and Y be two vector fields on the same manifold M. 

Definition 9.1 The Lie derivative ofY with respect to X is jCxY = m^t^\t=Q' "^^^"""^ ^ 
flow of X. 

The Lie derivative Cx^ is "the derivative of Y in the direction given by X." Its definition is 
coordinate-independent. By contrast, DY ■ X (also written as -'^[i^]) is also "the derivative of Y in 
the X direction" , but the value of DY ■ X depends on the coordinate system, and in particular does 
not usually equal CxY in the chosen coordinate system. 

Theorem 9.2 CxY = [X^Y\, where the bracket on the right is the Jacobi-Lie bracket. 

Proof. In the following calculation, we assume that M is finite-dimensional, and we work in local 
coordinates. Thus we may consider everything as matrices, which allows us to use the product rule 
and the identities (M-^)' = -M'^M'M-^ and | {D^t{x)) = D (|$t) (x). 



^^xY{x) = j^^y{x) 

d 



t=0 



dt 



~{D^t{x)r^ i^j^D^t{x)^ {D^t{x)r^Y{^t{x)) 



+ {D^t{x))-^ ^Y [^t{x)) 



^D^t{x)\Y{x) + ^Y{^t{x)) 



-D 



t=o 



Y {x) + DY{x) 



t=0 



= -DX{x) ■Y{x) + DY{x) ■ X{x) 
= [X,Y]j_l{x) 



Therefore CxY = [X, Y]j^l ■ 



Vorticity dynamics 

The same formula applies in infinite dimensions, although the proof is more elaborate. For example, 
the equation for the vorticity dynamics of an Euler fluid with velocity u (with div u = 0) and vorticity 
Lo = curl u may be written as, 

dtu = — u • Vo; -I- oj ■ Vu 

= - [u, w] 
= — aduo; 

= -CuLO 
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All of these equations express the invariance of the vorticity vector field oj under the flow of its 
corresponding divergenceless velocity vector field u. This is also encapsulated in the language of fluid 
dynamics in characteristic form as 

d / d \ dx 1 

dtr-a^J=°' ^ = u(x,t)=curr 

Here, the curl-inverse operator is defined by the Biot-Savart Law, 

u = curP^w = curl(— A)~"'^t(; , 

which follows from the identity 

curl curlu = — Au + Vdiv u , 

and application of div u = 0. Thus, in coordinates, 
dx 

— = u(x,t) x(t, xo) = $tXo with (^Q = Id, that is x(0,xo) = xq at t = , 

^ _ A/ ^ ^ 

Consequently, 



and a;-^'(x(t,xo),t) ^^.^^ ^ = w^(xo)^^o$/. 



^t* a;^(x(t, xo), t) = a;^(xo)^^^^ =: D^t ■ u; . 



^0 

This is the Cauchy (1859) solution of Euler's equation for vorticity, 

^ = [uj, curr^w] . 

This type of equation will reappear several more times in the remaining lectures. In it, the vorticity 
uj evolves by the ad-action of the right-invariant vector field u = curl~^a;. That is, 

du 

The Cauchy solution is the tangent lift of this flow, namely, 

$t*a;($t(xo)) = rxo$t(a;(xo)) . 
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10 Handout: Summary of Euler's equations for incompressible flow 
• Euler's equation of incompressible fluid motion 



ut + u • Vu + Vp = 

du/dt along dx/dt=u 
(advective time derivative) 

where u : x M ^ satisfies div u = 0. 
Geometric dynamics of vorticity 

uj = curl u 

ojt = — u • Vcj + uj ■ Vu 

= -[u,w] 
= — aduu; 
= ->CuW 

In these equations, one denotes ^ = ^ + -^u and, hence, may write Euler vorticity dynamics 
equivalently in any of the following three forms 



as well as 

dt ' V"' a 



9 ^ \ f d 



or 



-(a;.V)=0 along - = u 
The last form is found using the chain rule as 

-(...V) = -.V + u;.-V=(-u;-a;.Vu).V = 0. 

Ertel's theorem [Er 1942] The operators d/dt and u; • V commute on solutions of Euler's fluid 
equations. That is. 



d 

— ,U3-V 

dt' 



, so that -^(uj ■ VA) = u ■ V-^A for all differentiable A 
dt dt 



when u; = curl u and u is a solution of Euler's equations for incompressible fluid flow. Conse- 
quently, one finds the following infinite set of conservation laws: 



dt 



0, then J ^{u} ■ VA) d^x = const for all differentiable $ 
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Ohktani's formula [Ohkl993] 

(fu d 



, cj • Vu) 

dt^ dt ' ' 

„du 

= uj ■ V — 
dt 

= -u> ■ VVp 
= -Pa; 



where 



^ij = ^ . ( "Hessian" of pressure) 



In addition, one has the relations 



p = -A"4r(Vu^ • Vu) 

S = ^(Vu + Vu^) (strain rate tensor) 



so that, the following system of equations results. 



-Fuj 



duj 

~dt 

Kelvin (1890's) circulation theorem 

d d 

UJ • = OJ r 

du ^ 

where divu = 0, or cquivalently w' ■ = in index notation. The motion equation may be 
rewritten equivalently as a 1-form relation, 

dx^ = —dp = ViP dx* along = u 
dt dt 

— (uidx^) — Ui -rdx^ = -dp 
dt dt 



=Widu'=(i|up/2 



Consequently, 



which becomes 



dt 



(u-dx) = -d(^p- ^) 



^ u-(ix = — ® d { p 
dt Jc(u) Jc{u) 



u2 
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upon integrating around a closed loop C(u) moving with velocity u The 1-form relation above 
may be rewritten as 

{dt + Cu){u ■ dx) = -d\p - — 

whose exterior derivative yields using = 

{dt + C^){uj-dS) = Q 
where w • dS = curl u • dS = d (u • dx) . 

For these geometric quantities, one sees that the characteristic, or advective derivative is equiv- 
alent to a Lie derivative. Namely, 

d 



dt 



advect 



fluids geometry 

Stokes theorem 

The classical theorem due to Stokes 



Jd 



u • dx = / / curl u • dS 
dS J Js 



shows that Kelvin's circulation theorem is equivalent to conservation of flux of vorticity 

j^JJ a; • dS = , with dS = C(u) 

through any surface comoving with the flow. 
Recall the definition, 

d , 
(jO- — -i d'^x = u) ■ dS 
ox 



One may check this formula directly, by computing 

d d d \ 

ox^ ox/ ox-^ J ' 



= oj^dx^ A dx^ + J^dx^ A dx^udx^ A dx^ 
= a;-dS 

One may then use the vorticity equation in vector-field form, 

(9t + £u)u;-— =0 
ax 

to prove that the flux of vorticity through any comoving surface is conserved, as follows. 

= = div u d-^x = 

That is, as computed above using the exterior derivative 

{dt + Ji:u)i^-dS = 0. 
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• Momentum conservation 

From Euler's fluid equation dui/dt + Vip = with u-'^ = one finds, 

J {dfUi + u^djUi + dip)d\ = 

Ui d\ + j dj {uiU^ + pSf) d\ 



d 
dt 

d^ 
dt' 

= 0, if fi • u = 



-Mi + jnj{uiU^+p5i)dS 



Local conservation of fluid momentum is expressed using differentiation by parts as 

dtUi = -djT^ 
where := UiU^ +pSj is the fluid stress tensor. 

Moreover, each component of the total momentum Mj = Juid\ for i = 1,2,3, is conserved 
for an incompressible Euler flow, provided the flow is tangential to any fixed boundaries, i.e., 
n • u = 0. 

Mass conservation 

For mass density Z)(x, i) with total mass / I?(x,i) d^x, along dsi/dt = u{x,t) one finds, 

d\ = {dt + C^){Dd\) = {dtD + divDu) d\ = 
at •- ' 

continuity eqn 

The solution of this equation is written in Lagrangian form as 

{Dd\) ■ g-'^{t) = D{yio)d\ 
For incompressible flow, this becomes 

1 (J ^ ^"'^ I 
D dxQ d^Q 

Likewise, in the Eulerian representation one flnds the equivalent relations, 

5,i^ + divPu)=0 }^^^^" = ° 

Energy conservation 

Euler's fluid equation for incompressible flow divu = 

5tu + u • Vu + Vp = 
conserves the total kinetic energy, deflned by 

KE= I ^\u\'^d^x 
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The vector calculus identity 

u • Vu = — u X curlu + ^V|up 

recasts Euler's equation as 

dfU — uxcurlu + V p + — 



So that 

d lul 



dt 2 

Consequently, 



-T- / ^-^d^ = - i (p+ u • dS = 



since u • dS = u • n dS = on any fixed boundary and one finds 

KE = J ^\u\^d^x = const 

for Euler fluid motion. 
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11 Lie group action on its tangent bundle 



Definition 11.1 A Lie group G acts on its tangent bundle TG by tangent lifts.^ Given X G T^G we 
can consider the action of G on X by either left or right translations, denoted as T^LgX or T^RgX, 
respectively. These expressions may be abbreviated as 

ThLgX = L;X = gX and T^RgX = R*gX = Xg . 

Left action of a Lie group G on its tangent bundle TG is illustrated in the figure below. 



TL, 



TG 



TG 



X 



G 



G 



For matrix Lie groups, this action is just multiplication on the left or right, respectively. 

Left- and Right-Invariant Vector Fields. A vector field X on G is called left-invariant, if for 
every g & G one has L*X = X, that is, if 

inLg)X{h) = X{gh) , 

for every h E G. The commutative diagram for a left-invariant vector field is illustrated in the figure 
below. 

TLg 

TG - TG 



X 



Ln 



X 



G 



G 



Proposition 11.2 The set Xl{G) of left invariant vector fields on the Lie group G is a subalgebra of 
X{G) the set of all vector fields on G. 

Proof. If X, y e Xl{G) and g £ G, then 

l;[x,y] = [l;x,l;y] = [x,y], 

Consequently, the Lie bracket [^, y] € Xl{G). Therefore, Xl{G) is a subalgebra of X{G), the set of 
all vector fields on G. 

Proposition 11.3 The linear maps Xl{G) and TgG are isomorphic as vector spaces. 



^Recall definition 6.3 of tangent lifts of a differentiable manifold. 
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Demonstration of proposition. For each ^ G T^G, define a vector field on G by letting 
X^{g) = TeLg{i) . Then 

X^{gh) = TeLgh{0 = Te{LgoLhm 

= ThLginihiO) = ThLg{X^{h)) , 

which shows that X^ is left invariant. (This proposition is stated in Chapter 9 of [MaRal994], who 
refer to [AbMal978] for the full proof) 

Definition 11.4 [Jacobi-Lie bracket of vector fields] Let g{t) and h{s) be curves in G with 
g(0) = e, /i(0) = e and define vector fields at the identity of G by the tangent vectors g'{0) = ^, 
h'{0) = rj. Compute the linearization of the Adjoint action of G on TeG as 

■= jti'^'^^'^'^'^'^'i.o,^^^ 

This is the Jacobi-Lie bracket of the vector fields ^ and rj. 
Definition 11.5 The Lie bracket in TeG is defined by 

:= [X5,X,](e), 

for r] G TgG and for [X^,Xj^] the Jacobi-Lie bracket of vector fields. This makes T^G into a Lie 
algebra. Note that 

for all ^,rie TeG. 

Definition 11.6 The vector space TgG with this Lie algebra structure is called the Lie algebra of 
G and is denoted by g. 

If we let ^l{9) = TeLg^, then the Jacobi-Lie bracket of two such left-invariant vector fields in fact 
gives the Lie algebra bracket: 

For the right-invariant case, the right hand side obtains a minus sign, 

The relative minus sign arises because of the difference in action (xh~^)g~^ = x{gh)~^ on the right 
versus {gh)x = g{hx) on the left. 

Infinitesimal Generator. In mechanics, group actions often appear as symmetry transformations, 
which arise through their infinitesimal generators, defined as follows. 

Definition 11.7 Suppose ^ : G x M ^ M is an action. For C £ 0; ^^{t,x) : M. x M ^ M defined 
by = $(expt^, x) = $expt{(-^) is an M.— action on M. In other words, ^expte, —>■ M is a flow on 

M. The vector field on M defined by^ 

d 



t=0 ^ ''^ ' 



is called the infinitesimal generator of the action corresponding to ^. 

The Jacobi-Lie bracket of infinitesimal generators is related to the Lie algebra bracket as follows: 

[^m,Vm] = -[^,v]m- 
See, for example. Chapter 9 of [MaRal994] for the proof 



Recall Definition 5.22 of vector fields. 
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12 Lie algebras as vector fields 

Definition 12.1 (The ad-operation) For A ^ g we define the operator adA to be the operator 
ad: 9 X Q ^ Q that maps B & g to [A, B]. We write adAB = [A, B]. 

Definition 12.2 A representation of a Lie algebra q on a vector space V is a mapping p from g to 
the linear transformations of V such that for A,B ^ q and any constant scalar c, 

{i) p{A + cB) = p{A) + cp{B) 

{ii) p{[A,B]) = p{A)p{B)-p{B)p{A). 

If the map p is 1-1 the representation is said to faithful. 

Exercise 12.3 For a Lie algebra g, show that the map A — > adA is a representation of the Lie algebra 
g, with g itself the vector space of the representation. This is called the adjoint representation. 

Example 12.4 (Vector field representations of Lie algebras) The Jacobi-Lie bracket of the vec- 
tor fields ^ and rj in definition 9.2 may be represented in coordinate charts as, 



dx 
^ ds 



dx 

= v(x) , and 6 = — = u(x) . 
s=o ^ ^' ^ dt t=o ^ ' 



The Jacobi-Lie bracket of these two vector fields yields a third vector field, 



drj 




d^ 




dv dx 




du dx 


dt 


t=o 


ds 


s=0 


dx dt 


t=o 


dx ds 



dv du ^ ^ 

= — ■ u — ■ V = u ■ \v — V ■ \u . 

s=o dx dx 

Thus, the Jacobi-Lie bracket of vector fields at the tangent space of the identity T(.G is closed and may 
be represented in coordinate charts by the Lie bracket (commutator of vector fields) 

[C> »?] := - = u ■ Vv - V ■ Vu =: [u, v] . 

This example also proves the following 

Proposition 12.5 Let X(M") be the set of vector fields defined on M". A Lie algebra g may be 
represented on coordinate charts by vector fields = G X(M"') for each element ^ G 0. This 

vector field representation satisfies 

where [^,rj] E g is the Lie algebra product and [X^,X^] is the vector field commutator. 

13 Lagrangian and Hamiltonian Formulations 

13.1 Newton's equations for particle motion in Euclidean space 
Newton's equations 

ruiiii = Fj , i = 1, . . . , N , (no sum on i) (13.1) 
describe the accelerations qj of particles with 

Masses , i = 1, . . . , N , 



Euclidean positions q := (qi, . . . , qjv) € 



r,3N 



D. D. Holm 



Imperial College London 



Applications of Poisson Geometry 49 



in response to prescribed forces, 

F = (Fi,...,F^), 

acting on these particles. Suppose the forces arise from a potential. That is, let 

F,(q) = - ^^:^1M) , V-.R^^'^R, (13.2) 
dcii 

where dV/dq^i denotes the gradient of the potential with respect to the variable q^. Then Newton's 
equations (13.1) become 

dV 

mm = - — , i = l,...,N. (13.3) 

Remark 13.1 Newton (1620) introduced the gravitational potential for celestial mechanics, now called 
the Newtonian potential, 

^( W) = E • (13.4) 

13.2 Equivalence Theorem 

Theorem 13.2 (Lagrangian and Hamiltonian formulations) Newton's equations in potential form, 

dV 

mhi = -^, i = l,...,N, (13.5) 

for particle motion in Euclidean space R^^ are equivalent to the following four statements: 
(i) The Euler- Lagrange equations 

d ( dL\ dL ^ , ,^ 

hold for the Lagrangian L : R^^ = {(q, q) | q, q G M^^} — R , defined by 

N 

i^(q,q):=ETll^^ll'"^^^^' ^^^"^^ 



i=l 



with ||qi|| = qi • qi = QiQfSjk (no sum on i). 

(ii) Hamilton's principle of stationary action, SS = 0, holds for the action functional ( drop- 

ping i's) 

5[q(-)] := / L{c^{t)A{t))dt. (13.8) 

(iii) Hamilton's equations of motion, 

. dH . dH 

hold for the Hamiltonian resulting from the Legendre transform, 

^(q, P) := P • q(q, p) - i^(q, q(q, p)) , (13.10) 
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where q(q, p) solves for q from the definition p := (?L(q, q)/(?q. 

In the case of Newton's equations in potential form (13.5), the Lagrangian in equation (13.7) 
yields Pi = miC\i and the resulting Hamiltonian is (restoring i's) 

^ 1 

V ' Potential 

Kinetic energy 

(iv) Hamilton's equations in their Poisson bracket formulation, 

F = {F, H] for all F G jr(p) , (13.11) 
hold with Poisson bracket defined by 

(-.o}^-E(i;-g-g-g) /»^«« (13.12) 

Wc will prove this theorem by proving a chain of linked equivalence relations: (13.5) <S=> (i) <S=> (ii) 
(iii) 4^ (iv) as propositions. (The symbol 4^ means "equivalent to".) 



Step I. Proof that Newton's equations (13.5) are equivalent to (i). 

Check by direct verification. 



Step II. Proof that (i) 4^ (ii): The Euler- Lagrange equations (13.6) are equivalent to 
Hamilton's principle of stationary action. 



To simplify notation, we momentarily suppress the particle index i. 



We need to prove the solutions of (13.6) are critical points (5<S = of the action functional 

rb 



5[q(-)] := / L(q(t),q(t))dt, (13.13) 

J a 

the space of smooth trajectories 



(where q = dq{t)/dt) with respect to variations on C°^{[a,b], 
q : [a,b] — M.^^ with fixed endpoints qa, q^. 

In C°°{[a, 6], R^^) consider a deformation q(t, s), s G (— e, e), e > 0, with fixed endpoints q„, q^, 
of a curve qo(i), that is, q(t, 0) = qo(0 for all t G [a, b] and q(a, s) = qo(a) = qa, <lib, s) = q,o{b) = q^ 
for all s G (— e, e). 

Define a variation of the curve qo(-) in C°°([a, 6],M^^) by 



ds 



s=0 



q(-,.)Grq„(.)C-([a,6],: 



and define the first variation of S at qo(t) to be the following derivative: 

d 



55:=D5[qo(-)](^q(-)) ■■= 



ds 



5[q(-,s)]. 



(13.14) 



s=0 



D. D. Holm 



Imperial College London 



Applications of Poisson Geometry 



51 



Note that (5q(a) = 5q(6) = 0. With these notations, Hamilton's principle of stationary action 
states that the curve qo(i) satisfies the Euler-Lagrange equations (13.6) if and only if qo(-) is a critical 
point of the action functional, that is, D5[qo(-)] = 0. Indeed, using the equality of mixed partials, 
integrating by parts, and taking into account that (^q(a) = 5q(6) = 0, one finds 

/ L(q(i,s),q(i,s))di 
3=0-' a 



i=l 




■ dq^i dt = 



for all smooth 5qi(t) satisfying (5qi(a) = 5qi(5) = 0. This proves the equivalence of (i) and (ii), upon 
restoring particle index i in the last two lines. 



dS := D5[qo(-)](<^q(-)) 



ds 



.5=0 



ds 



dL_ 
dcii 



5<ii{t,s) 



dL 



■ 



dt 



Step III. Proof that (ii) <^ (Hi): Hamilton's principle of stationary action is equivalent 
to Hamilton's canonical equations. 

Definition 13.3 The conjugate momenta for the Lagrangian in (13.7) are defined as 

dL 

Pi := — ^ = m^qj G M^, i = l,...,N, (no sum on i) (13.15) 
c'qi 

Definition 13.4 The Hamiltonian is defined via the change of variables (q, q) i—>- (q, p), called the 
Legendre transform, 



H{q., p) : = p • q(q, p) - L(q, q(q, p)) 

N 



N 



i=i ^ 

• ^ ■' Potential 



Kinetic energy 



Remark 13.5 The value of the Hamiltonian coincides with the total energy of the system. This value 
will be shown momentarily to remain constant under the evolution of Euler-Lagrange equations (13.6). 

Remark 13.6 The Hamiltonian H may be obtained from the Legendre transformation as a function 
of the variables (q, p), provided one may soZue /or q(q,p), which requires the Lagrangian to be regular, 
e.g., 

d'^L 

det ^ . ^ . ^ (no sum on i) . 

dqidqi 

Lagrangian (13.7) is regular and the derivatives of the Hamiltonian may be shown to satisfy, 

dH 1 . dqi , dH dV dL 

opi mi dt dqi dqi dqi 



D. D. Holm 



Imperial College London 



Applications of Poisson Geometry 



52 



Consequently, the Euler-Lagrange equations (13.6) imply 

dpi d ( dL \ dL dH 



dt dt \ dcii J d(ii dc\i 

These calculations show that the Euler-Lagrange equations (13.6) are equivalent to Hamilton's 
canonical equations 

dH , dH , , 

qi = ^ , Pi = - ^ , 13.17 

where dH/dqi,dH/dpi € are the gradients of H with respect to qj,Pi G M^, respectively. This 
proves the equivalence of (ii) and (iii). 

Remark 13.7 The Euler-Lagrange equations are second order and they determine curves in con- 
figuration space G C°°{[a,b],M.^^). In contrast, Hamilton's equations are first order and they 
determine curves in phase space (qj,Pi) G C°°{[a,b],'R^^), a space whose dimension is twice the 
dimension of the configuration space. 

Step IV. Proof that (iii) <^=> (iv): Hamilton's canonical equations may be written using 
a Poisson bracket. 

By the chain rule and (13.17) any F G -F(P) satisfies, 

dF ^ fdF . dF 



\ ^ ; ur , ur 

dt ^ \dqi dpi 

1=1 

^ 'dF dH dF dH 



^d(\i dpi dpi dq.. 

1=1 

This finishes the proof of the theorem, by proving the equivalence of (iii) and (iv). 

Remark 13.8 (Energy conservation) Since the Poisson bracket is skew symmetric, {H, F} = 
— {F,H}, one finds that H = {H,H} = 0. Consequently, the value of the Hamiltonian is preserved 
by the evolution. Thus, the Hamiltonian is said to be a constant of the motion. 

Exercise 13.9 Show that the Poisson bracket is bilinear, skew symmetric, satisfies the Jacobi identity 
and acts as derivation on products of functions in phase space. 

Exercise 13.10 Civen two constants of motion, what does the Jacobi identity imply about additional 
constants of motion? 

Exercise 13.11 Compute the Poisson brackets among 

Ji — ^ijkPjQk 

in Euclidean space. What Lie algebra do these Poisson brackets recall to you? 



Exercise 13.12 Verify that Hamiltons equations determined by the function {J{z),^) = C • (q x p) 
give infinitesimal rotations about the ^— axis. 
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14 Hamilton's principle on manifolds 

Theorem 14.1 (Hamilton's Principle of Stationary Action) Let the smooth function L : TQ 
M 6e a Lagrangian on TQ. A curve c : [a, 6] — > Q joining qa = c{a) to qb = c{b) satisfies the Euler- 
Lagrange equations if and only if 



6 j L{c{t),c{t))dt = 0. 

J a 



Proof. The meaning of the variational derivative in the statement is the following. Consider a family 
of curves c{t,s) for \s\ < e satisfying co(t) = c{t), c{a,s) = qa, and c(6, s) = qb for all s G {—e,e). 
Then 

S L{c{t),c{t))dt := — / L{c{t,s),c{t,s))dt. 

J a ds g_Q J a 

Differentiating under the integral sign, working in local coordinates (covering the curve c{t) by a finite 
number of coordinate charts), integrating by parts, denoting 



, , d 



c{t,s), 

s=0 



and taking into account that v{a) = v{b) = 0, yields 



/ 

J a 



, ^'^v' + ]dt= I ( ^ - ) v'dt. 

\ oq^ oq^ J Ja \ 'J't dt oq^ J 



This vanishes for any function v{t) if and only if the Euler-Lagrange equations hold. 

Remark 14.2 The integral appearing in this theorem 

rb 



S{c{-)) := [ L{c{t),c{t))dt 
J a 



is called the action integral. It is defined on curves c : [a, 6] — > Q with fixed endpoints, c(a) = qa 
and c{b) = qi,. 

Remark 14.3 Variational derivatives of functionals vs Lie derivatives of functions. The 

variational derivative of a functional S[u] is defined as the linearization 



S[u + ev] - S[u] _ d_ 
e de 



lim 



S[u + ev] = ,v 

e=0 \ OV 



5S 



Compare this to the expression for the Lie derivative of a function. If f is a real valued function on a 
manifold M and X is a vector field on M, the Lie derivative of f along X is defined as the directional 
derivative 

Cxf = X{f) :=d/-X. 

// M is finite-dimensional, this is 

Cxf = X[f] ■.= df.X=^X^ = ]im /(^ + - /(^) . 

The similarity is suggestive: Namely, the Lie derivative of a function and the variational derivative 
of a functional are both defined as linearizations of smooth maps in certain directions. 
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The next theorem emphasizes the role of Lagrangian one-forms and two-forms in the variational 
principle. The following is a direct corollary of the previous theorem. 



Theorem 14.4 Given a C^ Lagrangian L : TQ 
£C{L) : Q T*Q, where 



for k > 2, there exists a unique C^ ^ 



map 



Q:-- 



t=0 



e T{TQ) q{t) is a curve in Q 



is a suhmxmifold ofT{TQ), and a unique C^ ^ one- form 0l G A^{TQ), such that for all variations 
q{t,s) (defined on a fixed t-interval) of q{t,0) = qQ{t) := q{t), we have 



SS 



d_ 
ds 



s=0 



S[c{-,s)\=T>S[q{-)\-5q{- 



£C{L){q,q,q)-Sqdt + 0L{q,q)-6q 



(14.1) 



where 



Sq{t) 



ds 



cf. Noether Thm 

q{t, s) . 



s=0 
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15 Summary Handout for Differential Forms 

Vector fields and l-forms 

Let M be a manifold. In what follows, all maps may be assumed to be C°°, although that's not 
necessary. 

A vector field on M is a map X : M ^ TM such that X{x) G T^M for every x e M. The set of 
all smooth vector fields on M is written X{M). ("Smooth" means differentiable or C" for some r < oo, 
depending on context.) 

A (difFerential) 1-form on M is a map 9 : M ^T*M such that 9{x) eT*M for every x e M. 

More generally, if tt M is a bundle, then a section of the bundle is a map (p : M ^ E such 

that TT o ip(x) = X for all x e M. So a vector field is a section of the tangent bundle, while a 1-form is 
section of the cotangent bundle. 

Vector fields can added and also multiplied by scalar functions /e : M — > M, as follows: {Xi + X2) (x) ■ 
Xi{x) + X2{x), {kX) {x) = k{x)X{x). 

Differential forms can added and also multiplied by scalar functions k : M ^ M., as follows: 
(a + p) {x) = a{x) + (3{x), [kO) {x) = k{x)e{x). 

We have already defined the push-forward and pull-back of a vector field. The pull-back of a 
1-form 6 an. N hy a map </? : M — > A?" is the 1-form ip*6 on M defined by 

{^*e) {x)-v = e{ip{x))-T^{v) 

The push-forward of a 1-form a on M by a diffeomorphism i/j : M ^ N is the pull-back of a by 

A vector field can be contracted with a differential form, using the pairing between tangent and 

cotangent vectors: (X _l 9) {x) = 9{x) • X{x). Note that X -\9 \s a map from M to M. Many books 
write ix9 in place of X _l ^, and the contraction operation is often called interior product. 

The differential of / : M ^ M is a 1-form df on M defined by 

df{x) ■ V = ^ f{c{t) 

t=0 

for any x G M, any v G T^M and any path c{t) in M such that c(0) = and c'(0) = v. The left hand 
side, df{x) ■ v, means the pairing between cotangent and tangent vectors, which could also be written 
df{x){v) or {df{x),v) . 
Note: 

X ^df = £xf = X[f] 

Remark 15.1 df is very similar to Tf, but Tf is defined for all differentiable f : M ^ N, whereas 
df is only defined when X = M (in this course, anyway). In this case, Tf is a map from TM to TM, 
and Tf{v) = df{x) ■ v G 7/(a,)]R for every v G T^M (we have identified Tf{x)R with M..) 

In coordinates... Let M be n-dimensional, and let x^, . . . be differentiable local coordinates for 
M. This means that there's an open subset U of M and an open subset V of M" such that the map 
ip : U ^ V defined by ip{x) = [x^{x), . .., x"'{x)) is a diffeomorphism. In particular, each x^ is a map 
from M to R, so the differential dx^ is defined. There is also a vector field ^ for every i, which is 
defined by ^{x) = ^<P~^ ifi^) + t^i)\^^Q , where e^ is the i*^ standard basis vector. 

Exercise 15.2 Verify that 

(where = means the left hand side is a constant function with value SV 
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Remark 15.3 Of course, given a coordinate system if = (x^, . . . , x**) , it is usual to write x = 
[x^, . . . jX") , which means x is identified with (x-'-(x), . . . ,x"'(x)) = <^(x). 

For every x G M, the vectors ^(x) form a basis for TxM, so every v € T'^M can be uniquely ex- 
pressed as f = v^^{x). This expression defines the tangent- lifted coordinates x^, . . . ,x"',v^, . . .v"^ 
on TM (they are local coordinates, defined on TU C TM). 

For every x € M, the covectors dx^{x) form a basis for T*M, so every a G T^M can be 
uniquely expressed as a = aj(ix*(x). This expression defines the cotangent-lifted coordinates 
x^, . . . , x*^, ai, . . . on T*M (they are local coordinates, defined on T*U C T*M). 

Note that the basis (^) is dual to the basis (dx^, . . . , dx") , by the previous exercise. It follows 
that, 

d 



(ojcZx*) 



v'-^ \ = ay 
ax* 



(we have used the summation convention) . 

In mechanics, the configuration space is often called Q, and the lifted coordinates are written: 
(on TQ) and g^, . . . . . . ,p„ (on T*Q). 

Why the distinction betw^een subscripts and superscripts? This is to keep track of how 
quantities vary if coordinates are changed (see next exercise) . One benefit is that using the summation 
convention gives coordinate-independent answers. 

Exercise 15.4 Consider two sets of local coordinates and s* on Q, related by (s^,...,s"^) = 
■0 ((?^, . . . , q") . Verify that the corresponding tangent lifted coordinates and are related by 

' = ■ 

Note that the last equation can be written as s = Dip{q)ci, where s is the column vector {s^, . . . s"), 
and similarly /or q. 

Do the corresponding calculation on the cotangent bundle side. See Definition 16.11. 

The next level: TTQ,T*T*Q, et cetera 

Since TQ is a manifold, we can consider vector fields on it, which arc sections of TiTQ). In 
coordinates, every vector field on TTQ has the form X = o' l^ + -j^; where the a' and 6* are 

functions of q and q. Note that the same symbol g* has two interpretations: as a coordinate on TQ 
and as a coordinate on Q, so can mean a vector field TQ (as above) or on Q. 

The tangent lift of the bundle projection r : TQ — > Q is a map Tt : TTQ — > TQ. If X is written in 
coordinates as above, then TtoX = o*^- A vector field X on TTQ is second order if TtoX{v) = v; 

in coordinates, a' = q^. The name comes from the process of reducing of second order equations to 

first order ones by introducing new variables = — — . 

at 

One may also consider T*TQ,TT*Q and T*T*Q. However, the subscript /superscript distinction 
is problematic here. 

1-forms 

The 1-forms on T*Q are sections of T*T*Q. Given cotangent-lifted local coordinates 



{q^,...,q'^,pi,...,Pn) 
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on T*Q, the general 1-form on T*Q has the form aidq^ + bidpi, where and bi are functions of {q,p). 
The canonical 1-form on T*Q is 

= Pidq\ 

d d 

also written in the short form p dq. Pairing 9(q,p) with an arbitrary tangent vector v = a*— -^ + 6*—^ G 

og* op^ 

T^g,p)T*Q gives 

= {p,Tt*{v)), 

where r* : T*Q Q ]s projection. In the last line we have interpreted q^ as a coordinate on Q, which 
implies that Pidq^ = p, by definition of the coordinates pi. Note that the last line is coordinate- free. 

2- forms 

Recall that a 1-form on M, evaluated at a point x G M, is a linear map from TxM to M. 
A 2-form on M, evaluated at a point x G M, is a skew-symmetric bilinear form on TxM; and the 
bilinear form has to vary smoothly as x changes. (Confusingly, bilinear forms can be skew-symmetric, 
symmetric or neither; differential forms are assumed to be skew-symmetric.) 

The pull-back of a 2-form a; on by a map : AI ^ N is the 2-form (p*LU on M defined by 

{ip*Lu) {x) {v,w) = e{ip{x)) {Tip{v),T^{w)) 
The push-forwEird of a 2-form a; on M by a diffeomorphism ij; : M ^ N is the pull-back of u by 

A vector field X can be contracted with a 2-form uj to get a 1-form X _l u; defined by 

{X _i u) {x){v) = u{x) {X{x),v) 

for any v G TxM. A shorthand for this is {X _l a;) (v) = co{X, v), or just X = io{X, •). 
The tensor product of two 1-forms a and /? is the 2-form a (g) /3 defined by 

(a (8) /?) {v, w) = a{v)P{w) 

for all V, we T*M. 

The wedge product of two 1-forms a and (5 is the skew-symmetric 2-form a A /? defined by 

(a A /3) {v, w) = a{v)(3{w) — a{w)P{v) . 



Exterior derivative 

The differential df of a real- valued function is also called the exterior derivative of /. In this context, 
real-valued functions can be called 0-forms. The exterior derivative is a linear operation from 0-forms 
to 1-forms that satisfies the Leibniz identity, a.k.a. the product rule. 



d{fg) = f dg + gdf 
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The exterior derivative of a 1-form is an alternating 2-form, defined as follows: 

d (ttidx'^) = ^^dx^ A dx*. 

Exterior derivative is a linear operation from 1-forms to 2 forms. The following identity is easily 
checked: 

d{df) = 

for all scalar functions /. 
n-forms 

See [MaRal994], or [Le2003], or [AbMal978]. Unless otherwise specified, n-forms are assumed to 
be alternating. Wedge products and contractions generalise. 

It is a fact that all n-forms are linear combinations of wedge products of 1-forms. Thus we can 
define exterior derivative recursively by the properties 

d{aA/3) = daAl3 + i-'\.)^a A dp, 

for all fc-forms a and all forms and 

dod = 

In local coordinates, if a = ai^...ii^dx^^ A • • • A (sum over alHi < • • • < ik), then 

da = ^^^dx^ A cZx^i A • • • A dx'" 
ox^ 

The Lie derivative of an n-form 6 in the direction of the vector field X is defined as 

where is the flow of X. 

Pull-back commutes with the operations d, _l , A and Lie derivative. 
Csirtan's magic formula: 

£xa = d{X ^a) + X ^da 
This looks even more magic when written using the notation ixo. = X _\a : 

£x = dix + ixd 

An n-form a is closed if da = 0, and exact \i a = djS for some (3. All exact forms are closed 
(since do d = 0), but the converse is false. It is true that all closed forms are locally exact; this is the 
PoincEire Lemma. 

Remark 15.5 For a survey of the basic definitions, properties, and operations on differential forms, 
as well as useful of tables of relations between differential calculus and vector calculus, see, e.g.. Chapter 
2 of [B12004]. 
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16 Euler-Lagrange equations of manifolds 

In Theorem 14.4, 



6S 



d 
ds 



where 



S[c{-,s)]=IiS[q{-)]-6q{-) 

s=0 

f SC{L){q,q,q)-5qdt + @Liq,q)-dq'' (16.1) 
Ja V 

cf. Noether Thm 

d 



q{t, s) , 
s=0 



the map SC : Q ^ T*Q is called the Euler-Lagrange operator and its expression in local coordinates 

One understands that the formal time derivative is taken in the second summand and everything is 
expressed as a function of {q,q,q)- 

Theorem 16.1 Noether (1918) Symmetries and Conservation Laws If the action variation in 
(16.1) vanishes 5S = Q because of a symmetry transformation which does not preserve the end points 
and, the Euler-Lagrange equations hold, then the term marked cf. Noether Thm must also vanish. 
However, vanishing of this term now is interpreted as a constant of motion. Namely, the term, 

dL 

Aiy^w) := (¥L{v),w) , or, in coordinates A{q,q,5q) = -gTjSq^ , 

is constant for solutions of the Euler-Lagrange equations. This result first appeared in Noether [Nol918j. 
In fact, the result in [Nol918] is more general than this. In particular, in the PDE (Partial Differential 
Equation) setting one must also include the transformation of the volume element in the action prin- 
ciple. See e.g. [O12000] for good discussions of the history, framework and applications of Noether' s 
theorem. 

Exercise 16.2 Show that conservation of energy results from Noether 's theorem if, in Hamil- 
tons principle, the variations are chosen as 



q{t, s) , 

s=0 



corresponding to symmetry of the Lagrangian under reparametrizations of time along the given curve 
q{t)^q{T{t,s)). 

The canonical Lagrangian one- form and two-form. The one-form Ql, whose existence and 
uniqueness is guaranteed by Theorem 14.4, appears as the boundary term of the derivative of the 
action integral, when the endpoints of the curves on the configuration manifold are free. In finite 
dimensions, its local expression is 

dL 



&L {q, q) ■■= ^ dg* = Pi {q, q) dg*) . 
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The corresponding closed two- form = dOi, obtained by taking its exterior derivative may be 
expressed as 

Ql ■■= -dOL = Q^'^ft ^ + Q^'^'l' A d^^ = dpi {q, q) A dq'j . 
These coefficients may be written as the 2n x 2n skew-symmetric matrix 



A 



^L=[ , (16.2) 

where A is the skew-symmetric n x n matrix ^ g^iggj j ~ ( d^idqj ) • Non-degeneracy of J^l is equivalent 
to the invertibility of the matrix ^ g^f^j 



*dqJ 

Definition 16.3 The Legendre transformation ¥L : TQ T*Q is the smooth map near the 
identity defined by 

{¥L{v,),w,) := ^ 



L{Vq + SWq) . 

s=0 



In the finite dimensional case, the local expression of ¥L is 

¥Liq\e) = (^q\^^ = iq\piiq,q)). 

If the skew-symmetric matrix (16.2) is invertiblc, the Lagrangian L is said to be regular. In this 
case, by the implicit function theorem, ¥L is locally invertible. If ¥L is a diffeomorphism, L is called 
hyperregular. 

Definition 16.4 Given a Lagrangian L, the action of L is the map A : TQ — M given by 

dL 

A{v) := {¥L{v),v) , or, in coordinates A{q,q) = 'Q^iq^ , (16.3) 
and the energy of L is 

dL 

E{v) := A{v) — L{v) , or, in coordinates E{q,q) = -QT^q^ — L{q,q) . (16.4) 

16.1 Lagrangian vector fields and conservation laws 

Definition 16.5 A vector field Z on TQ is called a Lagrangian vector field if 

nL{v){Z{v),w) = {dE{v),w), 

for all V e TqQ, w e T^{TQ). 

Proposition 16.6 The energy is conserved along the flow of a Lagrangian vector field Z. 
Proof. Let v{t) G TQ be an integral curve of Z. Skew-symmetry of implies 

jEHt)) = {dE{v{t)),v{t)) = {dE{v{t)),Z{vm 
= nL{v{t)){Z{v{t)),Z{v{t))) = Q. 

Thus, E{v{t)) is constant in t. m 
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16.2 Equivalence of dynamics for hyperregular Lagrangians and Hamiltonians 

Recall that a Lagrangian L is said to be hyperregular if its Legendre transformation ¥L : TQ T*Q 
is a diffeomorphism. 

The equivalence between the Lagrangian and Hamiltonian formulations for hyperregular Lagrangians 
and Hamiltonians is summarized below, following [MaRal994]. 

(a) Let L be a hyperregular Lagrangian on TQ and H = E o (FL)~^, where E is the energy of L 
and (FL)^^ : T*Q ^ TQ is the inverse of the Legendre transformation. Then the Lagrangian 
vector field Z on TQ and the Hamiltonian vector field Xh on T*Q are related by the identity 

{¥L)*Xh = Z. 

Furthermore, if c(t) is an integral curve of Z and d{t) an integral curve of Xh with FL(c(0)) = 
d{0), then ¥L{c{t)) = d{t) and their integral curves coincide on the manifold Q. That is, 
= ■^Q(<^(i)) = 7(*)) where tq : TQ — Q and ttq : T*Q — Q are the canonical bundle 
projections. 

In particular, the pull back of the inverse Legendre transformation ¥L~^ induces a one-form 
and a closed two-form O on T*(5 by 

@ = {¥L-^)*@L, n = - d@ = {¥L-^)*nL . 

In coordinates, these are the canonical presymplectic and symplectic forms, respectively, 

@ = Pidq\ n = - dQ = dpiA dg' . 



(b) A Hamiltonian H : T*Q ^ R is said to be hyperregular if the smooth map ¥H : T*Q ^ TQ, 
defined by 

H{ag + sPg), ag,PgeT*Q, 

s=0 



{¥H{ag),/3g) := 



ds 



is a diffeomorphism. Define the action of by G := (O , Xh)- If if is a hyperregular Hamil- 
tonian then the energies of L and H and the actions of L and H are related by 

E = H o {¥H)-'^, A = Go{¥H)-'^. 

Also, the Lagrangian L = A — E is hyperregular and ¥L = ¥H~^. 

(c) These constructions define a bijective correspondence between hyperregular Lagrangians and 
Hamiltonians. 

Remark 16.7 For thorough discussions of many additional results arising from the Hamilton's prin- 
ciple for hyperregular Lagrangians see, e.g. Chapters 7 and 8 of [MaRal994]. 

Exercise 16.8 (Spherical pendulum) A particle rolling on the interior of a spherical surface un- 
der gravity is called a spherical pendulum. Write down the Lagrangian and the equations of motion 
for a spherical pendulum with S'^ as its configuration space. Show explicitly that the Lagrangian is 

hyperregular. Use the Legendre transformation to convert the equations to Hamiltonian form. Find 
the conservation law corresponding to angular momentum about the axis of gravity by "bare hands" 
methods. 
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Exercise 16.9 (Differentially rotating frames) The Lagrangian for a free particle of unit mass 
relative to a moving frame is obtained by setting 

^(q,q,t) = ^||qf + q-R(q,i) 

for a function R(q, t) which prescribes the space and time dependence of the moving frame veloc- 
ity. For example, a frame rotating with time-dependent frequency Q{t) about the vertical axis z is 
obtained by choosing R(q, t) = q x n{t)z. Calculate Ql (QtQ) {q,q), the Euler-Lagrange operator 
£C{L) {q,q,q), the Hamiltonian and its corresponding canonical equations. 

Exercise 16.10 Calculate the action and the energy for the Lagrangian in Exercise 16.9. 

Definition 16.11 (Cotangent lift) Given two manifolds Q and S related by a diffeomorphism f : 
S, the cotangent lift T* f : T*S ^ T*Q of f is defined by 

{T*f{a),v) = {a,Tf{v)) (16.5) 

where 

a £ T*S , V & TgQ , and s = f{q) . 

As explained in Chapter 6 of [MaRal994], cotangent lifts preserve the action of the Lagrangian L, 
which we write as 

(p,q) = (a,s), (16.6) 

where p = T* f{a) is the cotangent lift of a under the diffeomorphism f and s = Tf{q) is the tangent 
lift ofq under the function f , which is written in Euclidean coordinate components as q^ ^ = f^{c\). 
Preservation of the action in (16.6) yields the coordinate relations, 

df^ ■ 

(Tangent lift in coordinates) = -q-jQ' =^ 

Qfk 

Pi = '^k-Q-j (Cotangent lift in coordinates) 

Thus, in coordinates, the cotangent lift is the inverse transpose of the tangent lift. 
Remark 16.12 The cotangent lift of a function preserves the induced action one-form, 

(P , dq) = (a , ds) , 
so it is a source of (pre-)symplectic transformations. 



16.3 The cleissic Euler-Lagrange example: Geodesic flow 

An important example of a Lagrangian vector field is the geodesic spray of a Riemannian metric. A 
Riemannian manifold is a smooth manifold Q endowed with a symmetric nondegenerate covariant 
tensor g, which is positive definite. Thus, on each tangent space TqQ there is a nondegenerate definite 
inner product defined by pairing with g{q). 

If {Q, g) is a Riemannian manifold, there is a natural Lagrangian on it given by the kinetic energy 
K of the metric g, namely, 

K{.v) ■= ^g{q){vg,Vq), 



D. D. Holm Imperial College London Applications of Poisson Geometry 63 



for q £ Q and Vq G TqQ. In finite dimensions, in a local chart, 

Kiq,q) = ^gij{q)q'q^- 

The Legendre transformation is in this case ¥K{vq) = g{q){vq, •), for Vq G TqQ. In coordinates, this is 



¥K{q,q) = 




= {q\gij{q)q^) =■ {q\Pi)- 



The Euler-Lagrange equations become the geodesic equations for the metric g, given (for finite 
dimensional Q in a local chart) by 

q' + ryH'' = 0, i = l,...n, 

where the three-index quantities 

are the Christoffel symbols of the Levi-Civita connection on {Q,g). 

Exercise 16.13 Explicitly compute the geodesic equation as an Euler-Lagrange equation for the ki- 
netic energy Lagrangian K{q,q) = \gij{q)q^q-^ ■ 

Exercise 16.14 For kinetic energy Lagrangian K{q, q) = \gij{q)q^q^ with i,j = l,2,..., N: 

• Compute the momentum pi canonical to if for geodesic motion. 

• Perform the Legendre transformation to obtain the Hamiltonian for geodesic motion. 

• Write out the geodesic equations in terms of q^ and its canonical momentum pi . 

• Check directly that Hamilton's equations are satisfied. 

Remark 16.15 A classic problem is to determine the metric tensors gij{q) for which these geodesic 
equations admit enough additional conservation laws to be integrable. 

Exercise 16.16 Consider the Lagrangian 

Mq,q) = ^l|qf-^(l-||q||Y 

for a particle in W^. Let '^fc{t) be the curve in obtained by solving the Euler-Lagrange equations for 
Le with the initial conditions qo = 7e(0), qo = 7e(0). Show that 

limjeit) 

is a great circle on the two-sphere S"^, provided that qo has unit length and the initial conditions satisfy 

qo • qo = 0. 

Remark 16.17 The Lagrangian vector field associated to K is called the geodesic spray. Since 
the Legendre transformation is a diffeomorphism ( in finite dimensions or in infinite dimensions if the 
metric is assumed to be strong), the geodesic spray is always a second order equation. 
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16.4 Covariant derivative 

The variational approach to geodesies recovers the classical formulation using covariant derivatives, 
as follows. Let X{Q) denote the set of vector fields on the manifold Q. The covariant derivative 

V : X{Q) X XiQ) ^ (X, Y) ^ Vx(^) , 

of the Levi-Civita connection on {Q,g) is given in local charts by 

VxiY) = T^.X'Y^^ + X'^^ . 
^ ^ »J dqk Qqi Qqk 

If c(t) is a curve on Q and Y £ the covariant derivative of Y along c{t) is defined by 

or locally, 

/ Jiv \ ^ 

TtAc{t)ym^{c{t)) + ^Y\c{t)). 



(^) = ^%mym^{c{t)) + ^^^^ 



A vector field is said to be parallel transported along c{t) if 



Dt 

Thus c{t) is parallel transported along c{t) if and only if 

c^ + rjitc^c'^ = 0. 

In classical differential geometry a geodesic is defined to be a curve c(t) in Q whose tangent vector 
c{t) is parallel transported along c{t). As the expression above shows, geodesies are integral curves of 
the Lagrangian vector field defined by the kinetic energy of g. 

Definition 16.18 A classical mechanical system is given by a Lagrangian of the form L{vq) = 
K{vq) — V{q), for Vq G TqQ. The smooth function V : Q R is called the potential energy. The 
total energy of this system is given hy E = K + V and the Euler- Lagrange equations ( which are always 
second order for a hyperregular Lagrangian) are 

BV 

where g^^ are the entries of the inverse matrix of {gij). 

Definition 16.19 If Q = M.^ and the metric is given by gij = 5ij, these equations are Newton's 
equations of motion (13.3) of a particle in a potential field which launched our discussion in Lecture 
10. 

Exercise 16.20 [Gauge invariance] Show that the Euler- Lagrange equations are unchanged under 

Lidit), q(t)) ^L' = L+ A^(q(t), q(i)) , (16.7) 
for any function 7 : M.^^ = {(q, q) | q, q € 
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Exercise 16.21 [Generalized coordinate theorem] Show that the Euler- Lagrange equations are un- 
changed in form under any smooth invertible mapping / : {q i— > s}. That is, with 

L(q(t),q(i)) = L(s(t),s(t)), (16.8) 

show that 

d (dL\ dL d ( dL\ dl ^ , ^ 

Exercise 16.22 How do the Euler-Lagrange equations transform under q(t) = r{t) + s{t) ? 

Exercise 16.23 (Other example Lagrangians) Write the Euler-Lagrange equations, then apply 
the Legendre transformation to determine the Hamiltonian and Hamilton's canonical equations for the 
following Lagrangians. Determine which of them are hyperregular. 

( ■ 

• L((i, q) = [gij{q)q^q-' ) (Is it possible to assume that L{q, q) = 1? Why?) 



l{q,q) = -(i-q-q) 



1/2 



• L{q,q) = • q + ^q • A(q), for constants m, c and prescribed function A(q). How do the 
Euler-Lagrange equations for this Lagrangian differ from free motion in a moving frame with 
velocity ^A(q) ? 

Example: Charged particle in a magnetic field. Consider a particle of charge e and mass m 
moving in a magnetic field B, where B = V x A is a given magnetic field on M^. The Lagrangian for 
the motion is given by the "minimal coupling" prescription (jay-dot-ay) 

^(q,q) = ^l|qf + ^A(q)-q, 

in which the constant c is the speed of light. The derivatives of this Lagrangian are 

dL • e . , dL .rp , 

— - = mq -I- -A =: p and — = - VA • q 
aq c oq c 

Hence, the Euler-Lagrange equations for this system are 

mq= ^(VA^-q- VA-q) = ^qx B 

c c 

(Newton's equations for the Lorentz force) 
The Lagrangian L is hyperregular, because 

P = ff'i(q, q) = mq + ^ A(q) 



has the inverse 



q = Fi/(q,p) = -(p--A(q)) . 

m \ c / 



m \ c 

The corresponding Hamiltonian is given by the invertible change of variables. 



1 



2 



J^(q,p) = P-q-Mq,q) = ^ p--a . (16.10) 

The Hamiltonian H is hyperregular since 

1 / e \ e 

q = ¥H(q, p) = — I p A ) has the inverse p = FL(q, q) = mq -|- - A . 

m \ c y c 

The canonical equations for this Hamiltonian recover Newton's equations for the Lorentz force law. 



D. D. Holm 



Imperial College London 



Applications of Poisson Geometry 



66 



Example: Charged particle in a magnetic field by the Kaluza-Klein construction. Al- 
though the minimal-couphng Lagrangian is not expressed as the kinetic energy of a metric, Newton's 
equations for the Lorentz force law may still be obtained as geodesic equations. This is accomplished 
by suspending them in a higher dimensional space via the Kaluza-Klein construction, which pro- 
ceeds as follows. 

Let Qkk be the manifold x with variables (q, 9). On Qkk introduce the one- form A + d9 
(which defines a connection one-form on the trivial circle bundle x 5^ — >■ M^) and introduce the 
Kaluza-Klein Lagrangian Lkk '■ TQkk — TM.^ x TS^ i-^ M as 



LKKici,e,ci,e) = ^m\\cif + ^ (A + de,ici,ci,e,e) 



2 



1 „-„2 

2'^\m + 



^(A-q+^)'. (16.11) 



The Lagrangian Lkk is positive definite in (q, ^); so it may be regarded as the kinetic energy of 
a metric, the Kaluza-Klein metric on TQkk- (This construction fits the idea of U{1) gauge 
symmetry for electromagnetic fields in R^. It can be generalized to a principal bundle with compact 
structure group endowed with a connection. The Kaluza-Klein Lagrangian in this generalization leads 
to Wong's equations for a color-charged particle moving in a classical Yang-Mills field.) The Legendre 
transformation for Lkk gives the momenta 

p = ynq-h(A-q-F^)A and 7r = A-q-|-^. (16.12) 

Since Lkk does not depend on 6, the Euler-Lagrange equation 

d OLkk QLkk 



dt de 99 



0, 



shows that vr = dLKK/d9 is conserved. The charge is now defined by e := ctt. The Hamiltonian 
Hkk associated to Lkk by the Legendre transformation (23.12) is 

Hkk{ci, 0, p, tt) = p • q -I- tt^ - Lkk{<1, q, 9, 6) 

= P ■ ^ (P ~ ttA) -|- 7r(7r — A ■ q) 
m 

1 II-II2 1 2 

— -m q — -TT 

2 2 

= P • — (P - 7I"A) + -TT 

m 2 

— ttA • — (p - ttA) - - — Hp - ttAII^ 

m 2m 

= i^||p-.Af + i.^. (16.13) 

On the constant level set vr = e/c, the Kaluza-Klein Hamiltonian Hkk is a function of only the 
variables (q, p) and is equal to the Hamiltonian (16.10) for charged particle motion under the Lorentz 
force up to an additive constant. This example provides an easy but fundamental illustration of 
the geometry of (Lagrangian) reduction by symmetry. The canonical equations for the Kaluza-Klein 
Hamiltonian Hkk now reproduce Newton's equations for the Lorentz force law. 
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17 The rigid body in three dimensions 



In the absence of external torques, Euler's equations for rigid body motion are: 



m2 



ih 



72)01^2, 



(17.1) 



or, equivalently, 

lO = lO X O , 

where ft = (fii, $12, ^^3) is the body angular velocity vector and h^h^h are the moments of inertia 
of the rigid body. 

Question 17.1 Can these equations - as they are written - he cast into Lagrangian or Hamiltonian 
form in any sense? (Since there are an odd number of equations, they cannot he put into canonical 
Hamiltonian form.) 

We could reformulate them as: 

Euler-Lagrange equations on rS0(3) or 

Canonical Hamiltonian equations on r*S0(3), 

by using Euler angles and their velocities, or their conjugate momenta. However, these reformulations 
on TS0(3) or T*S0(3) would answer a different question for a six dimensional system. We are 
interested in these structures for the equations as given above. 

Answer 17.2 (Lagrangian formulation) The Lagrangian answer is this: These equations may he 
expressed in Euler-Poincare form on the Lie algehra using the Lagrangian 



which is the (rotational) kinetic energy of the rigid hody. 

The Hamiltonian answer to this question will he discussed later. 

Proposition 17.3 The Euler rigid hody equations are equivalent to the rigid body action principle 
for a reduced action 



i{n) = -{ml + ml + isOi) = -fi^ . m 



(17.2) 




(17.3) 



where variations of fl are restricted to he of the form 



= E + X E 



(17.4) 



in which S(f) is a curve in that vanishes at the endpoints in time. 
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Proof. Since = and I is symmetric, we obtain 

/•b pb 

5 / i{n)dt = / {m,6ft)dt 

J a J a 



j (in, s + n X s) 

J a 

j 

J a 



a 



^m,5] ) + (m,n X x;) 

^ ' d 
— m + in X n, s; ) dt, 

dt 



dt 



upon integrating by parts and using the endpoint conditions, = S(a) = 0. Since S is otherwise 
arbitrary, (23.3) is equivalent to 

-^(m) + inxn = o, 

which are Euler's equations (17.1). ■ 

Let's derive this variational principle from the standard Hamilton's principle. 

17.1 Hamilton's principle for rigid body motion on TS0(3) 

An clement R G S0(3) gives the configuration of the body as a map of a reference configuration 
B C to the current configuration R(B); the map R takes a reference or label point X G B to a 
current point x = R(-'^) € R(B). 

When the rigid body is in motion, the matrix R is time-dependent. Thus, 

x{t) = K{t)X 

with R(t) a curve parameterized by time in S0(3). The velocity of a point of the body is 

x{t) = R(t)X = RR-^(t)x(i) . 

Since R is an orthogonal matrix, R~^R and RR~^ are skew matrices. Consequently, we can write 

(recall the hat map) 

i = RR~^x = w X X. (17.5) 

This formula defines the spatial angular velocity vector co. Thus, is essentially given by right 
translation of R to the identity. That is, the vector 

uj = {RR-^y. 

The corresponding body angular velocity is defined by 

ft = R-^a;, (17.6) 

so that fl is the angular velocity relative to a body fixed frame. Notice that 

R-^RX = R-^RR-^x = R"^(u; x x) 

= R-^u; X R-^x = flxX, (17.7) 
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so that fl is given by left translation of R to the identity. That is, the vector 

O = (R-^R)". 

The kinetic energy is obtained by summing up m\x\^ /2 (where | • | denotes the Euclidean norm) 
over the body. This yields 

K=l [ piX)\RXfd^X, (17.8) 



2JB 

in which p is a given mass density in the reference configuration. Since 

|RX| = |u; X x| = |R"^(a; x x)\ = \fl x X\, 
K is a quadratic function of ft. Writing 

K=^ft'^-Ift (17.9) 

defines the moment of inertia tensor I, which, provided the body docs not degenerate to a line, is 
a positive-definite (3 x 3) matrix, or better, a quadratic form. This quadratic form can be diagonalized 
by a change of basis; thereby defining the principal axes and moments of inertia. In this basis, we 
write I = diag(7i, /2, h). 

The function K is taken to be the Lagrangian of the system on TS0(3) (and by means of the 
Legendre transformation we obtain the corresponding Hamiltonian description on T*S0(3)). Notice 
that K in equation (17.8) is left (not right) invariant on rS0(3), since 

ft = {R-^Ay. 

It follows that the corresponding Hamiltonian will also be left invariant. 

In the framework of Hamilton's principle, the relation between motion in R space and motion in 
body angular velocity (or ft) space is as follows. 

Proposition 17.4 The curve R{t) G S0(3) satisfies the Euler-Lagrange equations for 

L(R,-R) = l f p{X)\-RX\^(fX, (17.10) 

if and only if ft{t) defined by R^-'^Rv = ft x v for all v G satisfies Euler's equations 

lft = lftxft. (17.11) 

The proof of this relation will illustrate how to reduce variational principles using their symmetry 
groups. By Hamilton's principle, R(i) satisfies the Euler-Lagrange equations, if and only if 

5 J L(R,R)dt = 0. 

Let l{ft) = ^{Ift) ■ ft, so that l{ft) = L(R, R) where the matrix R and the vector ft are related by 
the hat map, ft = (R~^R)". Thus, the Lagrangian L is left S0(3)-invariant. That is, 

l{ft) = L(R, R) = L(e, R-^R) . 
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To see how we should use this left-invariance to transform Hamilton's principle, define the skew matrix 
by fiv = fi X V for any v G M^. 

We differentiate the relation R~^R = ft with respect to R to get 

-R-\5R)R-^R + R-\5R) = m. (17.12) 

Let the skew matrix 'S be defined by 

S = R-^SR, (17.13) 

and define the vector S by 

Sv = S X V. (17.14) 

Note that 

t = - R-iRR-i^R + R-^SR, 

so 

R-i(5R = I: + r-^rI: . (17.15) 

Substituting (17.15) and (17.13) into (17.12) gives 

-tn + t + Cit = m, 

that is, 

5f2 = l: + [S7,i:]. (17.16) 

The identity [O, S] = (fi x E)" holds by Jacobi's identity for the cross product and so 

(5n = S + nxS. (17.17) 

These calculations prove the following: 

Theorem 17.5 For a Lagrangian which is left-invariant under S0(3), Hamilton's variational princi- 
ple 

5S = S [ L{R,R)dt = (17.18) 



on rS0(3) is equivalent to the reduced variational principle 

SSred = S [ /(fi) dt = (17.19) 

J a 

with Q = (R~^R)'' on where the variations 6Q are of the form 

(5n = s + n X s , 

with E(a) = S;(6) = 0. 
Recall Theorem 17.5: 

For a Lagrangian which is left-invariant under S0(3), Hamilton's variational principle 



6S 



= 5 j L{R,R)dt = (17.20) 

Ja 
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on rS0(3) is equivalent to the reduced variational principle 

dSred = S I l{n)dt = (17.21) 

J a 

with ft = (R^^R)' on where the variations SCt are of the form 

dn = t + Q.xT,, 

with X;(a) = = 0. 

Reconstruction of R(t) G 50(3). In Theorem 17.5, Euler's equations for the rigid body 
follow from the reduced variational principle (17.19) for the Lagrangian 

= ^{m) ■ n , (17.22) 

which is expressed in terms of the left- invariant time-dependent angular velocity in the body, fl G so(3). 
The body angular velocity ft{t) yields the tangent vector R{t) G Tj^(j')50(3) along the integral curve 
in the rotation group R(t) G 5*0(3) by the relation, 

R(t) = R{t)n{t) . 

This relation provides the reconstruction formula. It's solution as a linear differential equation 
with time-dependent coefficients yields the integral curve R(t) G 50(3) for the orientation of the rigid 
body, once the time dependence of 0(i) is determined from the Euler equations. 

17.2 Hamiltonian Form of rigid body motion. 

A dynamical system on a manifold M 

±{t) = F(x) , X G M 

is said to be in Hamiltonian form, if it can be expressed as 

x{t) = {x, iJ} , for H : M^R, 

in terms of a Poisson bracket operation, 

{•,•}: J^(M) X J^(M) ^ J^(M) , 

which is bilinear, skew-symmetric and satisfies the Jacobi identity and (usually) the Leibnitz rule. 

As wc shall explain, reduced equations arising from group-invariant Hamilton's principles on Lie 
groups are naturally Hamiltonian. If we Legendre transform our reduced Lagrangian for the 50(3) 
left invariant variational principle (17.21) for rigid body dynamics, then its simple, beautiful and 
well-known Hamiltonian formulation emerges. 

Definition 17.6 The Legendre transformation ¥1 : so(3) — so(3)* is defined by 

The Legendre transformation defines the body angular momentum by the variations of the rigid- 
body's reduced Lagrangian with respect to the body angular velocity. For the Lagrangian in (17.22), 
the components of the body angular momentum are 

Ui = Iini = —, i = 1,2,3. (17.23) 
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17.3 Lie-Poisson Hamiltonian formulation of rigid body dynamics. 

Let 

h{u) := {n,n)-i{n), 

where the pairing (•,•): so(3)* x so(3) — >■ M is understood in components as the vector dot product 
on 

{u,n) -n ^l. 

Hence, one finds the expected expression for the rigid-body Hamiltonian 

-1 

The Legendre transform ¥1 for this case is a diffeomorphism, so we may solve for 

dh ^ I dn\ / dl dn\ ^ 

In coordinates, this relation expresses the body angular velocity as the derivative of the reduced 
Hamiltonian with respect to the body angular momentum, namely (introducing grad- notation). 

Hence, the reduced Euler-Lagrange equations for I may be expressed equivalently in angular momen- 
tum vector components in M'^ and Hamiltonian h as: 

d 

— (m) = m X o <^ n = n X Vuh ■- {n, h} . 

This expression suggests we introduce the following rigid body Poisson bracket on functions of the 
n's: 

{/,^}(n):=-n-(Vn/x Vn/t). (17.25) 

For the Hamiltonian (17.24), one checks that the Euler equations in terms of the rigid-body angular 
momenta. 

Hi = ^-^f-^U2U3, 

n2 = ^^Vr^naHi, (17.26) 

il3 = ^^^^HiH2, 

that is, 

n = nxn. (17.27) 

are equivalent to 

f = {f,h}, with / = n. 

The Poisson bracket proposed in (17.25) is an example of a Lie Poisson bracket, which we will show 
separately satisfies the defining relations to be a Poisson bracket. 
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17.4 Poisson bracket. 

The rigid body Poisson bracket (17.25) is a special case of the Poisson bracket for functions on R^, 

{f,h} = -Vc-Vf xVh (17.28) 

This bracket generates the motion 

X = {x, /i} = Vc X V/i (17.29) 

For this bracket the motion takes place along the intersections of level surfaces of the functions c 
and /i in M^. In particular, for the rigid body, the motion takes place along intersections of angular 
momentum spheres c = ||x|p/2 and energy ellipsoids /i = x • Ix. (See the cover illustration of Marsden 
and Ratiu [2003].) 

Exercise 17.7 Consider the Poisson bracket 

{/, /j} = - Vc • V/ X V/i (17.30) 

Let c = x-^ • Cx be a quadratic form on W^, and let C be the associated symmetric 3x3 matrix. 
Determine the conditions on the quadratic function c(x) so that this Poisson bracket will satisfy the 
Jacobi identity. 

Exercise 17.8 Find the general conditions on the function c(x) so that the bracket 

{/, h} = -Vc-Vfx\7h 

satisfies the defining properties of a Poisson bracket. Is this bracket also a derivation satisfying 
the Leibnitz relation for a product of functions on ? If so, why? 

Exercise 17.9 How is the bracket related to the canonical Poisson bracket? Hint: restrict to level 
surfaces of the function c(x) . 

Exercise 17.10 (Casimirs of the bracket) The Casimirs (or distinguished functions, as Lie 
called them) of a Poisson bracket satisfy 

{c,/i}(x) = 0, V/i(x) 

Suppose the function c(x) is chosen so that the bracket (17.28) satisfies the defining properties of 
a Poisson bracket. What are the Casimirs for the bracket (17.28)? Why? 

Exercise 17.11 Show that the motion equation 

X = {x, h} 

for the bracket (17.28) is invariant under a certain linear combination of the functions c and h. 
Interpret this invariance geometrically. 

18 Momentum maps 
The Main Idea 

Symmetries are often associated with conserved quantities. For example, the flow of any 50(3)- 
invariant Hamiltonian vector field on T*R^ conserves angular momentum, q x p. More generally, 
given a Hamiltonian if on a phase space P, and a group action of G on P that conserves H, there is 
often an associated "momentum map" J : P ^ q* that is conserved by the flow of the Hamiltonian 
vector field. 

Note: all group actions in this section will be left actions until otherwise specified. 
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18.1 Hamiltonian systems on Poisson manifolds 

Definition 18.1 A Poisson bracket on a manifold P is a skew-symmetric bilinear operation on 
J^{P) := (P, M) satisfying the Jacobi identity and the Leibniz identity, 

{FG, H} = F{G, H} + {F, H}G 

The pair (P, {•, •}) is called a Poisson manifold. 

Remark 18.2 The Leibniz identity is sometimes not included in the definition. Note that bilinearity, 

skew-symmetry and the Jacobi identity are the axioms of a Lie algebra. In what follows, a Poisson 
bracket is a binary operation that makes J^{P) into a Lie algebra and also satisfies the Leibniz identity. 

Exercise 18.3 Show that the classical Poisson bracket, defined in cotangent-lifted coordinates 

{q^,...,q^,pi,...,PN) 

on an 2N -dimensional cotangent bundle T*Q by 

N 



' . \ dq^ dpi dpi dq^ J ' 



satisfies the axioms of a Poisson bracket. Show also that the definition of this bracket is independent 
of the choice of local coordinates (g^, . . . , q^^ . 

Definition 18.4 A Poisson map between two Poisson manifolds is a map ip : (Pi,{-,-}i) — ^ 
(P2,{-,-}2) that preserves the brackets, meaning 

{Foip,Go = {P, G}2 o if, for all F,G e T{P2) ■ 

Definition 18.5 An action ^ of G on a Poisson manifold (P, {,}) is canonical if is a Poisson 
map for every g, i.e. 

{F0^g,K0^g} = {F,K}0^g 

for every F,K e J^{P). 

Definition 18.6 Let (P, {•,•}) be a Poisson manifold, and let H : P ^ M be differentiable. The 
Hamiltonian vector field for H is the vector field Xh defined by 

Xh{F) = {P, H}, for any F G F{P) 

Remark 18.7 Xh is well-defined because of the Leibniz identity and the correspondance between 
vector fields and derivations (see [Le2003]). 

Remark 18.8 Xh(F) = £xh^ = -^^^ derivative of F along the flow of Xh- The equations 

F = {F,H}, 

called "Hamilton's equations" , have already appeared in Theorem 6.2, and are an equivalent definition 
ofXn. 

Exercise 18.9 Show that Hamilton's equations for the classical Poisson bracket are the canonical 
Hamilton's equations, 

dH . dH 



D. D. Holm 



Imperial College London 



Applications of Poisson Geometry 



75 



18.2 Infinitesimal invariance under Hamiltonian vector fields 

Let G act smoothly on P, and let ^ G g. Recall (from Lecture 9) that the infinitesimal generator is 
the vector field on P defined by 

for some path g{t) in G such that g{0) = e and ^'(0) = ^. 



t=o 



Remark 18.10 For matrix groups, we can take g{t) = exp(t^) . This works in general for the expo- 
nential map of an arbitrary Lie group. For matrix groups, 



(x) = ^ exp(tOx 



= (matrix multiplication) . 



Exercise 18.11 If H : P ~^ R is G-invariant, meaning that II{gx) = II{x) for all g e G and x e P, 
then £^pH = for all ^ G g. This property is called infinitesimal invariance. 

Example 18.12 (The momentum map for the rotation group) Consider the cotangent bundle 
of ordinary Euclidean space R^. This is the Poisson (symplectic) manifold with coordinates (q, p) € 
r*M^ ~ M®, equipped with the canonical Poisson bracket. An element g of the rotation group 50(3) 
acts on r*M^ according to 

5(q, P) = (^q, 9P) 

Set g{t) = exp(tA), so that J^|j^o^(*) ~ ^ ^'"'^ corresponding Hamiltonian vector field is 

^A = (q,p) = (^q,^p) 
where A G so(3) is a skew-symmetric matrix. The corresponding Hamiltonian equations read 

q = ^q = ^— , p = Ap = - — — . 

ap oq 

Hence, 

J A (q, p) = -Ap • q = ai€ijkPk<lj = a • q x p . 

for a vector a G with components ai, i = 1,2,3. So the momentum map for the rotation group is 
the angular momentum J = q x p. 

Example 18.13 Consider angular momentum J = q x p, defined on P = T*'M?. For every ^ G M^, 
define 

J( (q, P) := ? • (q X p) = p • (C x q) 

Using Exercise 18.9 and Example 18.10, 

dp ' 9q 
= (^ X q,^ X p) 

= ip (q, P) , 

where the last line is the infinitesimal generator corresponding to ^ ^ so(3). Now suppose H : P 
is SO (3) -invariant. From Exercise 18.11, we have £^H = 0. It follows that 

£xhJ^ = {J^, H} = — {H, J^} = —£xj^H = —£^pH = 0. 

Since this holds for all ^, we have shown that J is conserved by the Hamiltonian flow. 



Xj. (q, p) 



D. D. Holm 



Imperial College London 



Applications of Poisson Geometry 



76 



18.3 Defining Momentum Maps 

In order to generalise this example, we recast it using the hat map ' : ^ so(3) and the associated 
map~: (M^)* so{3)* , and the standard identification (M^)* = via the Euclidean dot product. 
We consider J as a function from P to so(3)* given by J (q, p) = (q x p) . For any ^ = v, we define 
(QjP) = ^(q ^ P) ) = (q X p) • V. As before, we find that Xj^ = for every ^, and J is 
conserved by the Hamiltonian flow. We take the first property, Xj^ = S,p, as the general definition of 
a momentum map. The conservation of J follows by the same Poisson bracket calculation as in the 
example; the result is Noether's Theorem. 

Definition 18.14 A momentum map for a canonical action of G on P is a map J : P ^ Q* such 
that, for every ^ € g, the map : P — > M defined by J^ip) = {J{p),0 satisfies 

Xj, = 

Theorem 18.15 (Noether's Theorem) Let G act canonically on (-P, {• , •}) with momentum map 
J. If H is G -invariant, then J is conserved by the flow of Xh- 

Proof For every ^ G 5, 

£xhJ^ = {J^, H} = — {H, Jjc} = —£xj^H = —£^pH = 0. ■ 

Exercise 18.16 Momentum maps are unique up to a choice of a constant element of g* on every 
connected component of M. 

Exercise 18.17 Show that the action on the torus := x given by a {6, 0) = (a + 0, (p) is 
canonical with respect to the classical bracket (with 6,(1) in place of q,p), but doesn't have a momentum 
map. 

Exercise 18.18 Show that the Petzval invariant for Fermat's principle in axisymmetric, translation- 
invariant media is a momentum map, r*M^ 1— >■ sp{2,M.)* taking (q, p) ^ {X,Y,Z). What is its 
symmetry? What is its Hamiltonian vector field? 

Theorem 18.19 (also due to Noether) Let G act on Q, and by cotangent lifts on T*Q. Then 
J : T*Q ^ Q* defined by, for every ^ G 0, 

{oiq) = {oiq, (q)) , for every aq G T*Q, 

is a momentum map (the "standard one") for the G action with respect to the classical Poisson bracket. 

(A proof using symplectic forms is given in Marsden and Ratiu [2003].) 

Proof We need to show that Xj^ = ^t*q^ for every ^ G g. From the definition of Hamiltonian 
vector fields, this is equivalent to showing that iT*Q[F\ = {F^ J^} for every F G J^{T*Q). We verify 
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this for finite-dimensional Q by using cotangent-lifted local coordinates. 



dp 



d 

Q I d 

Q-^{p,T^ieMtO)Qqi(l) 

d 

<T*Q{q,p),-Q-{q) 



t=o 



t=0 



{q,p) = -iT*Q{q,p) 



So for every F G J^{T*Q), 

d 



iT'Q[F\ = ^^F{exp{tOq,exp{tOp) 

OF ^ , , dF ^ , , 
= -Q<Q{.q) ^^T*Q[.q,p) 



t=0 



dt \dq 



t=0 



dt 



r*$(exp(te))P,7b(g) 



d_ 



t=o 



dF dJ^ dF dJj: 
dq dp dp dq 



Example 18.20 Let G C M„ (M) he a matrix group, with cotangent-lifted action on {q,p) G T*W^. 
For every g C M„(M), q gq. The cotangent-lifted action is {q,p) ^ {gq-,9~^p)- Thus, writing 
g = exp(t^), the linearization of this group action yields the vector field 

X^ = {Cq,-ep) 

The corresponding Hamiltonian equations read 

dJ^ J. dJ^ 

^'=^' ~^^=-^ 

This yields the momentum map J{q,p) given by 

Ji{q,p) = {J{q,p),0 = p^^Q (q) = p^^q ■ 

In coordinates, p^^q = PiCjq'^ > so J{q,p) = q^Pj- 

Exercise 18.21 Calculate the momentum map of the cotangent lifted action of the group of transla- 
tions ofR^. 

Solution 18.22 x G acts on q G by addition of vectors, 

X • (q) = q + X . 

The infinitesimal generator is linix^o Sc(^ + ~ Fhus, = Id and 

{Jk,0 = ((q>p)>?q) = {pJd) =PiSk = Pk 

This is also Hamiltonian with = p, so that {p, J^} = and {q, J^} = Id. 
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Example 18.23 Let G act on itself by left multiplication, and by cotangent lifts on T*G. We first 
note that the infinitesimal action on G is 

Let Jl be the momentum map for this action. Let ag G T*G. For every ^ & q, we have 

{JL{ag),0 = = {c^9,TRgO = {TR*gag,^) 

so Jl (dg) = TR*ag. Alternatively, writing ag = T*Lg-in for some fj, E Q* we have 

Jl {T*Lg-,ii) = TR*gT*Lg-.ii = Adl-,ii. 

Exercise 18.24 Show that the momentum map for the right multiplication action Rg{h) = kg is 
Jnio^g) = TL*gag. 

For matrix groups, the tangent lift of the left (or right) multiplication action is again matrix 
multiplication. Indeed, to compute TRg[A) for any A G TqSO{3), let B[t) be a path in SO{3) such 
that B{0) = Q and B'{0) = A. Then 



TRg{A) = j^B{t)G 



= AG. 



Similarly, TLg{A) = GA. To compute the cotangent lift similarly, we need to be able to consider 
elements of T*G as matrices. This can be done using any nondegenerate bilinear form on each tangent 
space TqG. We will use the pairing defined by 

{{A,B)) ■.= -\tr{A^B)=-\tr{AB^). 

(The equivalence of the two formulas follows from the properties tr{CD) = tr{DC) and tr{C'^) = 
tr{C)). 

Exercise 18.25 Check that this pairing, restricted to so(3), corresponds to the Euclidean inner product 
via the hat map. 

Example 18.26 Consider the previous example for a matrix group G. For any Q G G, the pairing 
given above allows use to consider any element P G TqG as a matrix. The natural pairing of TqG 
with TqG now has the formula, 

(P, A) = -]^tr {P'^A) , for all A G TqG. 

We compute the cotangent-lifts of the left and right multiplication actions: 

{T*Lq{P),A) = {P,TLq{A)) = {P,QA) 

= -hr [P^QA) = -\tr [{Q^ Pf a) = {Q^P,A) 
{T*Rq{P),A) = {P,TRq{A)) = {P,AQ) 



\tr {P^AQf) = -\tr {PQ^ A^) = {PQ^,A) 
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In summary, 

T*Lq{P) = Q^P and T*Rq{P) = PQ^ 

We thus compute the momentum maps as 

Jl{Q,P) = T*RqP = PQ^ 
Jr{Q,P) = T*LqP = Q^P 

In the special case of G = SO{3), these matrices PQ^ and Q^P are skew-symmetric, since they are 
elements of so{3). Therefore, 

Jl (Q, P) = T*RqP = ^ {PQ^ - QP^) 
Jr (g, P) = T*LqP = ^ {Q^P - P^Q) 

Exercise 18.27 Show that the cotangent lifted action on SO{n) is expressed as 

Q-P = Q^P 

as matrix multiplication. 

Definition 18.28 A momentum map is said to be equivariant when it is equivariant with respect to 
the given action on P and the coadjoint action on g* . That is, 

J{g ■p)= Ad;., J (p) 

for every g G G, p E P, where g ■ p denotes the action of g on the point p and where Ad denotes the 
adjoint action. 

Exercise 18.29 Show that the momentum map derived from the cotangent lift in Theorem 18.19 

is equivariant. 

Example 18.30 (Momentum map for symplectic representations) Let {V, O) be a symplectic 
vector space and let G be a Lie group acting linearly and symplectically on V . This action admits an 
equivariant momentum map J : F — > g given by 

J^{v) = {i{v),£,) = hl{i-v,v), 

where ^ ■ v denotes the Lie algebra representation of the element ^ £ q on the vector v £ V . To 
verify this, note that the infinitesimal generator ^v{v) = ^ ■ v, by the definition of the Lie algebra 
representation induced by the given Lie group representation, and that ■u.,v) = — r2(tt, ^ • v) for all 
u^v^V. Therefore 

dJ^{u){v) = ■u,v) + ■v,u) = • u, v). 

Equivariance of i follows from the obvious relation g^^ • • g • v = (Adg-i ^) • v for any 3 G G, ^ G g, 
and V eV. 
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Example 18.31 (Cayley-Klein parameters and the Hopf fibration) Consider the natural ac- 
tion of SU{2) on C^. Since this action is by isometrics of the Hermitian metric, it is automatically 
symplectic and therefore has a momentum map J : — ^ 5u(2)* given in example 18.30, that is, 

{3{z,w),^) = l-n{^ ■ {z,w),{z,w)), 



where z,w E C and ^ G su(2). Now the symplectic form on is given by minus the imaginary part of 

the Hermitian inner product. That is, C" has Hermitian inner product given by z ■ w := ^j^j^ 
where z = (zi, . . . , w = (wi, . . . ,Wn) G C". The Rym,plectic form is thus given by il(z.w) : = 
— Im(z • w) and it is identical to the one given before on M^" by identifying z = u + iv G C" with 
(u, v) G M^'* and w = u' + G C" with (u', v') G M^". 

The Lie algebra su(2) of SU{2) consists of 2x2 skew Hermitian m,atrices of trace zero. This Lie 
algebra is isomorphic to 50 (3) and therefore to (M^, x) by the isomorphism given by 



X = {x-^ , x'^ , x"^) G 



X 



—zx 
-ix^ + x^ 



-} 

ix"^ 



-IX — X 

r.3 



Gsu(2). 



Thus we have [x,y] = (x x y) for any x,y G M^. Other useful relations are det(2x) = ||x|p and 
trace(xy) = — jX • y. Identify su(2)* with by the map /x G su(2)* /i G defined by 



/i • X := 



-2(m,x) 



for any x G M^. With these notations, the momentum map J : 
coordinates: for any x G we have 



can be explicitly computed in 



3{z,w) • X = —2{3{z, w),x) 



= - Im 
2 



— IX 

—ix^ + x'^ 



ix" 



-IX — X 
n3 





z 




z 




w 




w 



1 



(2Yle{wz),21m.{wz), \z\ 



\w\ 



X. 



Therefore 



J{z, w) 



'-{2wz, \zf - \wf) G 



endowed 
has the — 



Thus, J is a Poisson map from C^, endowed with the canonical symplectic structure, to 
with the + Lie Poisson structure. Therefore, —J : — > «s a canonical map, if 1 
Lie-Poisson bracket relative to which the free rigid body equations are Hamiltonian. Pulling back the 
HamAltonian H(U) = U ■ I-'^U/2 to gives a Hamiltonian function (called collective) on C^. The 
classical Hamilton equations for this function are therefore projected by —J to the rigid body equations 
II = n X I~^n. In this context, the variables {z,w) are called the Cayley-Klein parameters. 



Exercise 18.32 Show that — JI53 : S S is the Hopf fibration. In other words, the momentum 
map of the SU{2)-action on C^, the Cayley-Klein parameters and the family of Hopf fibrations on 
concentric three-spheres in are all the same map. 



Exercise 18.33 Optical traveling wave pulses The equation for the evolution of the complex 
amplitude of a polarized optical traveling wave pulse in a material medium is given as 

. _ 1 dH 
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with Hamiltonian H : 



M defined by 




u = {u^,u^,u^) G u : 




,2 



—iu — w 



2 H 



Gsu(2). 



1. Prove that this isomorphism is an equivariant momentum map. 

2. Deduce the equations of motion for the Stokes vectors of this optical traveling wave and write it 
as a Lie Poisson Hamiltonian system. 

3. Determine how this system is related to the equations for an SO{3) rigid body. 

Exercise 18.34 The formula determining the momentum map for the cotangent-lifted action of a Lie 
group G on a smooth manifold Q may be expressed in terms of the pairing (■,■): g* x g i—> M as 



on the coordinate q. 

Define appropriate pairings and determine the momentum maps explicitly for the following actions, 

[a] £^q = ^xq for x 

[b] £^q = ad^q for ad-action ad : q x q q in a Lie algebra q 

[cj AqA~^ for A G GL(3, R) acting on q & GL{3, R) by matrix conjugation 
[dj Aq for left action of A e 50(3) on q e SO{3) 
[e] AqA^ for A G GL{3, R) acting on q E Sym(3), that is q = q^ . 
Answer 18.35 

[a] p-^xq = qxp-^^J = qxp. (The pairing is scalar product of vectors.) 

[b] {p , ad^q ) = — ( ad* p , C) =^ J = ad* p for the pairing (•,•): g* x g i-^ M 

[c] Compute T(,{AqA~^) = ~ = ol for £, = ^'(0) G 5/(3, R) acting on q £ GL{3, R) by matrix 

Lie bracket [■ , ■]. For the matrix pairing {A, B) = iTace{A^B), we have trace(p^[^, g]) = 
trace ((pg-^ — q^p)^Ci ^ J = pq^ — q^P- 

[d] Compute Te{Aq) = £^q for ^ = A'{0) G so{3) acting on q & SO{3) by left matrix multiplication. 

For the matrix pairing {A, B) = trace(^-^B), we have trace{p^^q) = tra£e{{pq^)^ ^) J = 
\{pq^ — q^p), where we have used antisymmetry of the matrix ^ G so{3). 

[e] Compute Te{AqA^) = ^q+q(^ for^ = ^'(0) G gl{3, R) acting onq £ Sym{3) . For the matrix pair- 

ing {A, B) = trace(^'^i?), we have trace{p'^ {^q -\- q^'^)) = trace(q'(p^ = trace(2qp)^^) =^ 
J = 2qp, where we have used symmetry of the matrix ^q + q^^ to choose p = p^ . (The momen- 
tum canonical to the symmetric matrix q = q^ should be symmetric to have the correct number 
of components!) 
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Equi variance 

Definition 18.36 A momentum map is Ad* -equivariant iff 

J {g ■ x) = Ad^.^J {x) 

for all g e G,x e P. 

Proposition 18.37 All cotangent-lifted actions are Ad* -equivariant. 

Proposition 18.38 Every Ad* -equivariant momentum map J : P q* is a Poisson map, with 
respect to the '-h' Lie-Poisson bracket on g*. 

19 Quick summary for momentum maps 

Let G be a Lie group, g its Lie algebra, and let g* be its dual. Suppose that G acts symplectically 
on a symplectic manifold P with symplectic form denoted by Q. Denote the infinitesimal generator 
associated with the Lie algebra element ^ by and let the Hamiltonian vector field associated to a 
function / : P — > M be denoted Xf, so that df = Xf _l fi. 

19.1 Definition, History and Overview 

A momentum map J : P — > g* is defined by the condition relating the infinitesimal generator of 
a symmetry to the vector field of its corresponding conservation law, (J, 0> 

for all ^ eg. Here (J, ^) : P — M is defined by the natural pointwise pairing. 

A momentum map is said to be equivariant when it is equivariant with respect to the given action 
on P and the coadjoint action on g*. That is, 

J{g-p)=Adl-,J{p) 

for every g G G, p E P, where g ■ p denotes the action of g on the point p and where Ad denotes the 
adjoint action. 

According to [Wel983], [Liel890] already knew that 

1 . An action of a Lie group G with Lie algebra g on a symplectic manifold P should be accompanied 
by such an equivariant momentum map J : P — ^ g* and 

2. The orbits of this action are themselves symplectic manifolds. 

The links with mechanics were developed in the work of Lagrange, Poisson, Jacobi and, later, Noether. 
In particular, Noether showed that a momentum map for the action of a group G that is a symmetry 
of the Hamiltonian for a given system is a conservation law for that system. 

In modern form, the momentum map and its equivariance were rediscovered in [Kol966] and 
[Sol970] in the general symplectic case, and in [Sml970] for the case of the lifted action from a 
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manifold Q to its cotangent bundle P = T*Q. In this case, the equivariant momentum map is given 
explicitly by 

(J(ag),0 = {oiq,^Qiq)) , 

where G T*Q, ^ G g, and where the angular brackets denote the natural pairing on the appropriate 
spaces. See [MaRal994] and [OrRa2004] for additional history and description of the momentum map 
and its properties. 

20 Rigid body equations on SO(n) 

Recall from [Manl976] and [RalQSO] that the left invariant generalized rigid body equations on SO(n) 
may be written as 

Q = Qn, 

M = Mn - QM =: [M, Q] , (RBn) 

where Q G SO(n) denotes the configuration space variable (the attitude of the body), $7 = Q~^Q G 
so{n) is the body angular velocity, and 

M := J{n) = D'^Vt + VLD'^ G so*{n) , 

is the body angular momentum. Here J : so{n) so{n)* is the symmetric (with respect to the above 
inner product) positive definite operator defined by 

J{n) = D^n + QD^, 

where is the square of the constant diagonal matrix D = diag {di, ^2, d^} satisfying df + djyO for 
alii ^ j. For n = 3 the elements of df are related to the standard diagonal moment of inertia tensor 
J by 

I = diag{/i,/2,/3}, Ii = dl + dl, l2 = dl + dl, h = dl + dl. 

The Euler equations for the SO{n) rigid body M = [M, fi] are readily checked to be the Euler- 
Lagrange equations on so{n) for the Lagrangian 

L(Q, Q) = i{n) = ^ {n, j{n)) , with n = q^q. 

The momentum is found via the Legendre transformation to be 

di 



and the corresponding Hamiltonian is 



r)l 1 
H{M) = — -^-m) = -{M,J-\M)). 

The quantity M is the angular momentum in the body frame. The corresponding angular momentum 
in space, 

m = QMQ^ , is conserved m = . 

Indeed, conservation of spatial angular momentum m implies Euler's equations for the body angular 
momentum M = Q^mQ = Adgm. 
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20.1 Implications of left invariance 

This Hamiltonian H{M) is invariant under the action of SO{n) from the left. The corresponding 
conserved momentum map under this symmetry is known from the previous lecture as 

Jl : T*SO{n) ^ so{n)* is Jl(Q, P) = PQ^ 

On the other hand, we know (from Lectures 18 & 19) that the momentum map for right action is 

Jr : T*SO{n) ^ so{nT , Jr{Q, P) = Q'^P 

Hence M = Q^P = Jr. Therefore, one computes 

H{Q,P) = H{Q,Q-M) = H{Id,M) (by left invariance) 

= H{M)=^-{M,J-\M)) 

= liQ^P^J-HQ^'P)) 

Hence, we may write the SO{n) rigid body Hamiltonian as 

H{Q,P) = ^{Q^P,n{Q,P)) 

Consequently, the variational derivatives of H{Q,P) = ^{Q'^P , ^{Q,P)) are 

6H = {Q^5P + 5Q^P,Q{Q,P)) 

= ti{dP^Qn) + ti{p'^6Qn) 

= ti{6P^Qn) + tr{6QnP^) 
= tr(5P^Q17) + tr((5Q'^P0'^) 

= {sp , Qn) - {SQ , pn) 

where skew symmetry of O is used in the last step, i.e., = —J). Thus, Hamilton's canonical 
equations take the form, 

Q = — = Qn, 

P=--=™. (20.1) 

Equations (20.1) are the symmetric generalized rigid body equations, derived earher in [BlCrl997] 
and [BlBrCrl997] from the viewpoint of optimal control. Combining them yields, 

Q-iQ = n = p-^P ^ {PQTy = , 

in agreement with conservation of the momentum map Jl{Q,P) = PQ^ corresponding to symmetry 
of the Hamiltonian under left action of SO{n). This momentum map is the angular momentum in 
space, which is related to the angular momentum in the body by PQ^ = m = QMQ'^. Thus, we 
recognize the canonical momentum as P = QM (see exercise 18.23), and the momentum maps for left 
and right actions as, 

Jl = m = PQ^ (spatial angular momentum) 
Jr = M = Q^P (body angular momentum) 

Thus, momentum maps TG* g* corresponding to symmetries of the Hamiltonian produce conser- 
vation laws; while momentum maps TG* i-^ q* which do not correspond to symmetries may be used 
to re-express the equations on g*, in terms of variables on TG*. 
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21 Manakov's formulation of the SO {4) rigid body 

The Euler equations on S'0(4) are 



dM 



MQ - nu = [M, n] 



(RBn) 



where and M are skew symmetric 4x4 matrices. The angular frequency is a linear function of 
the angular momentum, M. [Manl976] "deformed" these equations into 

d 



dt 



(M + XA) = [{M + XA), {n + XB)] , 



where A, B are also skew symmetric 4x4 matrices and A is a scalar constant parameter. For these 
equations to hold for any value of A, the coefficent of each power must vanish. 

• The coefficent of A^ is 

= [A,B] 

So A and B must commute. So, let them be constant and diagonal: 

Aij = diag(aj)^y , Bij = diag(6j)(5y 

• The coefficent of A is 



(no sum) 



dA 

"dt 



= ^ = [A,Q] + [M,5] 
Therefore, by antisymmetry of M and O, 

(ttj - aj)^ij = (pi - bj)Mij 



hj - bj 
a,- — ai 



Mi. 



(no sum) 



Finally, the coefficent of A° is the Euler equation, 

dM 



dt 



= [M,n], 



but now with the restriction that the moments of inertia are of the form, 



ni 



-Mi- 



(no sum) 



ai — aj 

which turns out to possess only 5 free parameters. 
With these conditions, Manakov's deformation of the S'0(4) rigid body implies for every power n that 



d^ 
dt 



(M + A^)" = [(M + XA)"", {n + XB)] 



Since the commutator is antisymmetric, its trace vanishes and one has 

^trace(M + XA)"" = 

after commuting the trace operation with time derivative. Consequently, 

trace(M + A^)" = constant 

for each power of A. That is, all the coefficients of each power of A are constant in time for the >S'0(4) 
rigid body. [Manl976] proved that these constants of motion are sufficient to completely determine 
the solution. 
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Remark 21.1 This result generalizes considerably. First, it holds for SO{n). Indeed, as as proven 
using the theory of algebraic varieties in [IIal984], Manakov's method captures all the algebraically 

integrable rigid bodies on SO{n) and the moments of inertia of these bodies possess only 2n — 3 param- 
eters. (Recall that in Manakov's case for S'0(4) the moment of inertia possesses only five parameters.) 
Moreover, [MiFol978] prove that every compact Lie group admits a family of left- invariant metrics 
with completely integrable geodesic flows. 

Exercise 21.2 Try computing the constants of motion trace(M + AA)" for the values n = 2, 3, 4. 
How many additional constants of motion are needed for integrability for these cases? How many for 
general n? Hint: keep in mind that M is a skew symmetric matrix, M'^ = —M, so the trace of the 
product of any diagonal matrix times an odd power of M vanishes. 

Answer 21.3 The traces of the powers trace(M + XA)'^ are given by 



n=2 



n=3 



n=4 



tr + 2Atr (AM) + A^tr 
tr + 3Atr (AM^) + 3A^tr A^M + Xhr A^ 

tr M"^ + 4Atr {AM^) + X'^{2ti A^M"^ + 4tr AM AM) + A^tr A^M + A^'tr 



The number of conserved quantities for n = 2,3,4 are, respectively, one (C\ = trM"^), one (I\ = 
trAM'^) and two (C2 = trM^ and h = 2tiA^M'^ + Air AM AM). The quantities Ci and C2 a,re 
Casimirs for the Lie-Poisson bracket for the rigid, body. Thus, {Ci. H} = = {C2,H} for any 
Hamiltonian H{M); so of course Ci and C2 are conserved. However, each Casimir only reduces the 
dimension of the system by one. The dimension of the original phase space is dimT* SO{n) = n(n— 1). 
This is reduced in half by left invariance of the Hamiltonian to the dimension of the dual Lie algebra 
dimso{n)* = n{n — l)/2. For n = 4, dimso{4)* = 6. One then subtracts the number of Casimirs 
(two) by passing to their level surfaces, which leaves four dimensions remaining in this case. The 
other two constants of motion /i and I2 turn out to be sufficient for integrability, because they are in 
involution {/i,/2} = and because the level surfaces of the Casimirs are symplectic manifolds, by the 
Marsden-Weinstein reduction theorem [MaWe74j. For more details, see [Ral980]. 

Exercise 21.4 How do the Euler equations look on so(4)* as a matrix equation? Is there an analog 
of the hat map for so(3)* ? Hint: the Lie algebra so(4) is locally isomorphic to so(3) x so(3). 

Exercise 21.5 Write Manakov's deformation of the rigid body equations in the symmetric form 
(20.1). 

22 Free ellipsoidal motion on GL(n) 

Riemann [Ril860] considered the deformation of a body in M" given by 

x{t,xo) = Q{t)xQ, (22.1) 

with a; , € M", Q{t) G G-L+(n, M) and x{tQ,XQ) = xq, so that Q{tQ) = Id. (The subscript + in 
GL_|-(n, R) means nxn matrices with positive determinant.) Thus, x(t,XQ) is the current (Eulerian) 
position at time t of a material parcel that was at (Lagrangian) position xo at time to. The "defor- 
mation gradient," that is, the Jacobian matrix Q = dx/dxo of this "Lagrange-to-Euler map," is a 
function of only time, t, 

dx/dxQ = Q{t) , with det Q > . 
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The velocity of such a motion is given by 

x{t, xo) = Q{t) xo = Q{t)Q-^ {t) X = u{t, x) . (22.2) 
The kinetic energy for such a body occupying a reference volume B defines the quadratic form, 

L=^j^p{xo)\x{t,xo)\''d^xo=\tr{Q{tflQ{t)) = \q\I^''Q'b- 

Here I is the constant symmetric tensor, 

= I p{xo)xoX^ d^xo, 
Jb 

which we will take as being proportional to the identity I^^ = C'^S^^ for the remainder of these 
considerations. This corresponds to taking an initially spherical reference configuration for the fiuid. 
Hence, we are dealing with the Lagrangian consisting only of kinetic energy,^ 



L 

2 



\tr(Q{tfQ{t)) . 

The Euler-Lagrange equations for this Lagrangian simply represent free motion on the group GL^(n, M), 

Q{t) = 0, 

which is immediately integrable as 

Q{t) = Q{0) + Q{0)t , 

where Q{0) and Q{0) are the values at the initial time t = 0. Legendre transforming this Lagrangian 
for free motion yields 

The corresponding Hamiltonian is expressed as 

H{Q,P) = \tr[p^p) = l\\Pf. 

The canonical equations for this Hamiltonian are simply 

Q = P, with P = . 

22.1 Polar decomposition of free motion on GL+(n, M) 

The deformation tensor Q{t) G GL+(n,M) for such a body may be decomposed as 

Q{t) = R-\t)D{t)S{t) . (22.3) 

This is the polar decomposition of a matrix in GL+(n,M). The interpretations of the various compo- 
nents of the motion can be seen from equation (22.1). Namely, 

• Re SO{n) rotates the x-coordinates. 



^[Ril860] considered the much more difficult problem of a self- gravitating ellipsoid deforming according to (22.1) in 
R^. See [Chl969] for the history of this problem. 
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• 5 G SO{n) rotates the xo-coordinates in the reference configuration^ and 

• D is a diagonal matrix which represents stretching deformations along the principal axes of the 
body. 

The two SO{n) rotations lead to their corresponding angular frequencies, defined by 

= RR-^ , A = SS-^ . (22.4) 

Rigid body motion will result, when S restricts to the identity matrix and -D is a constant diagonal 
matrix. 

Remark 22.1 The combined motion of a set of fluid parcels governed by (22.1) along the curve 
Q{t) € GL+(n,R) is called, "ellipsoidal," because it can be envisioned, in three d,im,ensions as a, fluid 
ellipsoid whose orientation in space is governed by R E SO{n), whose shape is determined by D 
consisting of its instantaneous principle axes lengths and whose internal circulation of material is 
described by S £ SO{n). In addition, fluid parcels initially arranged along a straight line within the 
ellipse will remain on a straight line. 



22.2 Euler-Poincare dynamics of free Riemann ellipsoids 

In Hamilton's principle, 6 J L dt = 0, we chose a Lagrangian L : TGL^(n,]S.) ^ M in the form 

LiQ,Q) = Tin,A,D,D), (22.5) 

in which the kinetic energy T is given by using the polar decomposition Q{t) = R~^{t)D{t)S{t) in 
(22.3), as follows. 

Q = R-^{-nD + D + DA)S. (22.6) 
Consequently, the kinetic energy for ellipsoidal motion becomes 



T = - trace 
2 



= - trace 
2 



-no^n - noD + qdad + ddq + 



-DA^D - DAD + DADn + DAD 

-n^D^ -A'^D'^ + 2nDAD +D' 
^ ^ ' 

Coriolis coupling 



(22.7) 



Remark 22.2 Note the discrete exchange symmetry of the kinetic energy: T is invariant under f2 -^-^ 

For A = and D constant expression (22.7) for T reduces to the usual kinetic energy for the 
rigid-body, 

T = -]-tv&ce\9.{D9. + 9.D)\ . (22.8) 

\=Q,D=const 4 L J 

This Lagrangian (22.5) is invariant under the right action, R — Rg and S — Sg, for g G SO{n). 
In taking variations we shall use the formulas* 



(50 = S + [S,J7] = E-adQE, T. = 5RR-^ 
SA = S + [S,A] = S-adAS, E = dSS-\ 



(22.9) 
(22.10) 



^This is the "particle relabeling map" for this class of motions. 
According to [Chl969] this discrete symmetry was first noticed by Riemann's friend, [Del860]. 
*These variational formulas are obtained directly from the definitions of CI and A. 
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in which the ad-operation is defined in terms of the Lie-algebra (matrix) commutator 
adnS = [f2,S]. Substituting these formulas into Hamilton's principle gives 



as, e.g. 



= djLdt = Ji 



dn OA dD dD 



= dt 



dL 
dt 





dL 






t, - ad^S 


OA ' 


E - adA^ 


+ 



+ 



d dL 
dtdn ~ 
d dL _dL 

dt~dD ~ dD 



6L 



6n 

6D, 



d dL 
didX 



dD 
dL 
dX 



d dL 
dtdb 



6D, 



(22.11) 



where, the operation ad^, for example, is defined by 

,,dL dL 



(22.12) 



and the dot ' ■ ' denotes pairing between the Lie algebra and its dual. This could also have been written 
in the notation using (•,•): g* x g — M as. 



dL 



, adnS ) = - 



dL 



dn' " / \dn 

The Euler-Poincare dynamics is given by the stationarity conditions for Hamilton's principle. 



5D 



d dL 

dtdn 
d dL 

dtdA 
d dL 

dt dD 



adj 

dL_ 

' dD 



dL 
'dA 

= 0. 



0, 



(22.13) 

(22.14) 
(22.15) 
(22.16) 



These are the Euler-Poincare equations for the ellipsoidal motions generated by Lagrangians of 
the form given in equation (22.5). For example, such Lagrangians determine the dynamics of the 
Riemann ellipsoids - circulating, rotating, self-gravitating fluid flows at constant density within an 
ellipsoidal boundary. 



22.3 Left and right momentum maps: Angular momentum versus circulation 

The Euler-Poincare equations (22.14-22.16) involve angular momenta defined in terms of the angular 
velocities n, A and the shape D by 



M = 
N = 



dT 

dn 

dT 
dX 



= -nD'^ - D^n + 2D AD , 
= -AL>2 _ jj2j^ _^ 2DnD . 



(22.17) 
(22.18) 



These angular momenta are related to the original deformation gradient Q = R ^DS in equation 
(22.1) by the two momentum maps from Example 18.26 



PQ' -QP' =QQ' -QQ' = R-^MR, 



P^Q-Q^P = Q^Q-Q^Q 



S-^NS. 



(22.19) 
(22.20) 
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To see that A'^ is related to the vorticity, we consider the exterior derivative of the circulation one-form 
u • dx defined as 

d{u ■ dx) = curlu • dS = ^{Q^Q- Q^Q)jk dx^ A dx!^ = {S-^NS)jkdxi A dx!^ . (22.21) 

Thus, S''^NS is the fluid vorticity, referred to the Lagrangian coordinate frame. For Euler's fluid 
equations, Kelvin's circulation theorem implies {S~^NS)' = 0. 

Likewise, M is related to the angular momentum by considering 

UiXj — UjXi = QikX^X^Q-^ — QikX^xlfQ^ . (22.22) 

For spherical symmetry, we may choose XqXq = 8^^ and, in this case, the previous expression becomes 

UiXj - UjXi = [QQ^ - QQ^]ij = [R-^MR]ij . (22.23) 

Thus, R~^MR is the angular momentum of the motion, referred to the Lagrangian coordinate frame 
for spherical symmetry. In this case, the angular momentum is conserved, so that (R^^MR)' = 0. 

In terms of these angular momenta, the Euler-Poincare-Lagrange equations (22.14-22.16) are ex- 
pressed as 

M = [n,M], (22.24) 
N = [A,N], (22.25) 
d fdL\ dL 



, , ■ , ■ (22.26) 

dt\dDj dD ^ ' 

Perhaps not unexpectedly, because of the combined symmetries of the kinetic-energy Lagrangian (22.5) 
under both left and right actions of SO{n), the first two equations are consistent with the conservation 
laws, 

[R-^MR)- = and {S'^NSy = , 

respectively. Thus, equation (22.24) is the angular momentum equation while (22.25) is the vorticity 
equation. (Fluids have both types of circulatory motions.) The remaining equation (22.26) for the 
diagonal matrix D determines the shape of the ellipsoid undergoing free motion on GL(n, M). 

22.4 Vector representation of free Riemann ellipsoids in 3D 

In three dimensions these expressions may be written in vector form by using the hat map, written 
now using upper and lower case Greek letters as, 

with ei23 = 1, and D = diag {di, d2, d^}. 

Exercise 22.3 What is the analog of the hat map in four dimensions? Hint: locally the Lie algebra 
so(4) is isomorphic to so(3) x so(3). 

Hence, the angular-motion terms in the kinetic energy may be rewritten as 

^ trace {U^D^) = j- \{dl + dl)ujl + {dj + dl)iof + {dj + d?)a;|l , (22.27) 



2 ^ ' 21 

1 '-2^2, 1 



- -trace {A'D') = ^idi + di)Xi + {di + d^)Af + {di + df)A^ , (22.28) 
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and 

- -trace (nDAD) = did2{uJzh) + d2d3(wiAi) + d:idi{uj2\2) ■ (22.29) 

On comparing equations (22.8) and (22.27) for the kinetic energy of the rigid body part of the motion, 
we identify the usual moments of inertia as 

Ik = d^ + dj , with i, j, k cyclic. 

The antisymmetric matrices M and N have vector representations in 3D given by 

dT 



Mk = ^ = idj + d^AuJk-2didj\k, (22.30) 



dT 

^fe = = {d^ + d''j)Xk - 2didjUk , (22.31) 

again with i, j, k cyclic permutations of {1, 2, 3}. 

Vector representation in 3D In terms of their 3D vector representations of the angular momenta 
in equations (22.30) and (22.31), the two equations (22.24) and (22.25) become 

M = {RR-^)M = J^M = a;xM, N = {SS-^)N = AN = A x N . (22.32) 

Relative to the Lagrangian fluid frame of reference, these equations become 

{R-'^My = i?"^(M - X M) = 0, (22.33) 
(5-iN)- = S-\N - XxN) = 0. (22.34) 

So each of these degrees of freedom represents a rotating, deforming body, whose ellipsoidal shape is 
governed by the Euler-Lagrange equations (22.26) for the lengths of its three principal axes. 

Exercise 22.4 (Elliptical motions vi^ith potential energy on GL(2,R)) Compute equations (22.24- 
22.26) for elliptical motion in the plane. Find what potentials V{D) are solvable for L = T{il, A, D, D) — 
V{D) by reducing these equations to the separated Newtonian forms, 

£^2 ^Y^^^ £^ dW{a) 



dt^ dr^ ' dt^ da ' 

for = + (i| and a = tan~^(d2/c?i) with di{t) and d2{t) in two dimensions. Hint: consider the 
potential energy, 

V{D) = v(trD^,det{D)j , 
for which the equations become homogeneous in r^(i). 

Exercise 22.5 (Ellipsoidal motions with potential energy on GL(3,M)) Choose the Lagrangian 

in 3D, 



L = -ti[Q'Qj -V[ti{Q'Q),dct{Q)j , 

where Q(t) G GL{3, M) is a 3 x 3 matrix function of time and the potential energy V is an arbitrary 
function of ti (Q^Q) and det(Q). 

1. Legendre transform this Lagrangian. That is, find the momenta Pij canonically conjugate to Qij, 
construct the Hamiltonian H{Q, P) and write Hamilton's canonical equations of motion for this 
problem. 
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2. Show that the Hamiltonian is invariant under Q — > OQ where O G SO{3). Construct the 
cotangent lift of this action on P. Hence, construct the momentum map of this action. 

3. Construct another distinct action of S0{3) on this system which also leaves its Hamiltonian 
H{Q, P) invariant. Construct its momentum map. Do the two momentum maps Poisson com- 
mute? Why? 

4. How are these two momentum maps related to the angular momentum and circulation in equa- 
tions (22.17) and (22.18)? 

5. How does the 2D restriction of this problem inform the previous one? 
Exercise 22.6 (GL(n, M)—inv£iriant motions) Begin with the Lagrangian 



where S is an n x n symmetric matrix and q G is an n— component column vector. 

1. Legendre transform to construct the corresponding Hamiltonian and canonical equations. 

2. Show that the system is invariant under the group action 



for any constant invertible n x n matrix, A. 

3. Compute the infinitesimal generator for this group action and construct its corresponding mo- 
m,entum map. Is this momentum map equivariant? 

4. Verify directly that this momentum map is a conserved n x n matrix quantity by using the 
equations of motion. 

5. Is this system completely integrable for any value of n > 2? 

23 Heavy top equations 
23.1 Introduction and definitions 

A top is a rigid body of mass m rotating with a fixed point of support in a constant gravitational 
field of acceleration —gz pointing vertically downward. The orientation of the body relative to the 
vertical axis z is defined by the unit vector F = R~^(i)z for a curve R{t) € SO{3). According to its 
definition, the unit vector T represents the motion of the vertical direction as seen from the rotating 
body. Consequently, it satisfies the auxiliary motion equation. 




q ^ and S ASA^ 



f = -R"^R(t)r = r X n. 



Here the rotation matrix R(t) G 50(3), the skew matrix S7 = R ^R G so(3) and the body angular 
frequency vector fi G are related by the hat map, Q, = (R~^R)", where " : (so(3), [•, ■]) — >■ (M^, x) 
with Civ = 17 X V for any v G M^. 



The motion of a top is determined from Euler's equations in vector form. 



m = mx^l-\-mgT XX, (23.1) 
f = rxO, (23.2) 

where O, T, % G are vectors in the rotating body frame. Here 
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• Q, = (Oi, 1^2, ^^3) is the body angular velocity vector. 

• I = diag(/i, /2, /a) is the moment of inertia tensor, diagonalized in the body principle axes. 

• r = R~^{t)z represents the motion of the unit vector along the vertical axis, as seen from the 
body. 

• X is the constant vector in the body from the point of support to the body's center of mass. 

• m is the total mass of the body and g is the constant acceleration of gravity. 

23.2 Heavy top action principle 

Proposition 23.1 The heavy top equations are equivalent to the heavy top action principle for a 
reduced action 

SSred = , with S'red = / K^, T) dt = [ J ( lO , O) - (mi? X , T) dt, (23.3) 

J a J a ^ 

where variations of fl and T are restricted to be of the form 

(5n = i; + n X E and 5T = T x , (23.4) 

arising from variations of the definitions ft = (R^^R)' and F = R~^(t)z in which 5](t) = (R~"^5R)'~ 
is a curve in that vanishes at the endpoints in time. 

Proof. Since I is symmetric and x is constant, we obtain the variation, 

rb 



sf l{n,r)dt= j {m,8Vl) -{mgx,5'^)dt 

J a Ja 

b 



I 

J a 

a- 



(m, S + O X S) - (m£?x, r x 

a 

b 



, x; ) + ( m , n X x;) - (mc/ x , r X T.)dt 
-^m + m X fi + mor xx,'^)dt, 

dt 



upon integrating by parts and using the endpoint conditions, S(6) = S(a) = 0. Since S is otherwise 
arbitrary, (23.3) is equivalent to 

- + I^lx^l + mg^ xx = o, 

dt 

which is Euler's motion equation for the heavy top (23.1). This motion equation is completed by the 
auxiliary equation F = F x fi in (23.2) arising from the definition of F. ■ 

The Legendre transformation for l(fl,T) gives the body angular momentum 
The well known energy Hamiltonian for the heavy top then emerges as 

h{n,T) = n-^i-i{n,r) = ^{n,r'n) + {mgx,r), (23.5) 

which is the sum of the kinetic and potential energies of the top. 
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The Lie-Poisson Equations. Let f,h : — > M be real- valued functions on the dual space q* . 
Denoting elements of g* by /x, the functional derivative of / at /it is defined as the unique element 
5f /SjJL oi Q defined by 

lim ^[/(/x + e5ii) - /(/x)] = (^^M, ^) , (23.6) 
for all G g*, where (• , •) denotes the pairing between g* and g. 



Definition 23.2 (Lie-Poisson brackets & Lie-Poisson equations) 
The (±) Lie-Poisson brackets are defined by 



{/,/i}±(m) = ±(/x, 



51 5h 

The corresponding Lie-Poisson equations, determined by / = {/, /i} read 



Sf 



where one defines the ad* operation in terms of the pairing (• , •), by 



(23.7) 



(23.8) 



{/, h} = [iJ,, adsh/si^ 



5/ 



The Lie-Poisson setting of mechanics is a special case of the general theory of systems on Poisson 
manifolds, for which there is now an extensive theoretical development. (See Marsden and Ratiu 
[2003] for a start on this literature.) 



23.3 Lie-Poisson brackets and momentum maps. 

An important feature of the rigid body bracket carries over to general Lie algebras. Namely, Lie- 
Poisson brackets on g* arise from canonical brackets on the cotangent bundle (phase space) T*G 
associated with a Lie group G which has g as its associated Lie algebra. Thus, the process by which 
the Lie-Poisson brackets arise is the momentum map 

T*G ^ g* . 

For example, a rigid body is free to rotate about its center of mass and G is the (proper) rotation 
group SO (3). The choice of T*G as the primitive phase space is made according to the classical 
procedures of mechanics described earlier. For the description using Lagrangian mechanics, one forms 
the velocity phase space TG. The Hamiltonian description on T*G is then obtained by standard 
procedures: Legendre transforms, etc. 

The passage from T*G to the space of II's (body angular momentum space) is determined by left 
translation on the group. This mapping is an example of a momentum map; that is, a mapping whose 
components arc the "Noether quantities" associated with a symmetry group. The map from T*G to 
g* being a Poisson map is a general fact about momentum maps. The Hamiltonian point of view of 
all this is a standard subject. 

Remark 23.3 (Lie-Poisson description of the heavy top) As it turns out, the underlying Lie 

algebra for the Lie-Poisson description of the heavy top consists of the Lie algebra se(3,M) of in- 
finitesimal Euclidean motions in . This is a bit surprising, because heavy top motion itself does not 
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actually arise through actions of the Euclidean group of rotations and translations on the body, since 
the body has a fixed point! Instead, the Lie algebra se(3, M) arises for another reason associated with 
the breaking of the S0(3) isotropy by the presence of the gravitational field. This symmetry breaking 
introduces a semidirect-product Lie-Poisson structure which happens to coincide with the dual of the 
Lie algebra se(3, M) in the case of the heavy top. As we shall see later, a close parallel exists between 
this case and the Lie-Poisson structure for compressible fluids. 

23.4 The heavy top Lie-Poisson brackets 

The Lie algebra of the special Euclidean group in 3D is se(3) = x with the Lie bracket 

[it u), (t7, v)] = X r/, ^ X V - r? X u) . (23.9) 

Wc identify tlic dual space with pairs (11, F); the corresponding (— ) Lie-Poisson bracket called the 
heavy top bracket is 

{/ , /i}(n, r) = - n • Vn/ xVnh-r- (Vn/ X Vrh - Vnh x Vr/) • (23.10) 

This Lie-Poisson bracket and the Hamiltonian (23.5) recover the equations (23.1) and (23.2) for the 
heavy top, as 

n = {n,h} = n X Vn^ + r X Vr/i = n X i~^n-i-r X m^x, 

t = {T,h} = T xVuh = T xr^U. 

Remark 23.4 (Semidirect products and symmetry breaking) The Lie algebra of the Euclidean 
group has a structure which is a special case of what is called a semidirect product. Here, it is the 

semidirect product action so(3)(DM^ of the Lie algebra of rotations so(3) acting on the infinitesimal 
translations M^, which happens to coincide with se(3,M). In general, the Lie bracket for semidirect 
product action q(M)V of a Lie algebra g on a vector space V is given by 



(X,a),(X,a)] = ([X,X],X{a)-X{a)') 



in which X,X £ q and a,a G V. Here, the action of the Lie algebra on the vector space is denoted, 
e.g., X{a). Usually, this action would be the Lie derivative. 

Lie-Poisson brackets defined on the dual spaces of semidirect product Lie algebras tend to occur 
under rather general circumstances when the symmetry in T*G is broken, e.g., reduced to an isotropy 
subgroup of a set of parameters. In particular, there are similarities in structure between the Poisson 
bracket for compressible flow and that for the heavy top. In the latter case, the vertical direction of 
gravity breaks isotropy of from SO (3) to 5*0(2). The general theory for semidirect products is 
reviewed in a variety of places, including [MaRaWeS^a, MaRaWeS^b]. Many interesting examples of 
Lie-Poisson brackets on semidirect products exist for fluid dynamics. These semidirect-product Lie- 
Poisson Hamiltonian theories range from simple fluids, to charged fluid plasmas, to magnetized fluids, 
to multiphase fluids, to super fluids, to Yang-Mills fluids, relativistic, or not, and to liquid crystals. 
See, for example, [GiHoKul982], [HoKul982], [HoKul983], [HoKul988]. For discussions of many of 
these theories from the Euler-Poincare viewpoint, see [HoMaRal998a] and [Ho2002a]. 

23.5 The heavy top formulation by the Kaluza-Klein construction 

The Lagrangian in the heavy top action principle (23.3) may be transformed into a quadratic form. 
This is accomplished by suspending the system in a higher dimensional space via the Kaluza-Klein 
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construction. This construction proceeds for the heavy top as a slight modification of the well-known 
Kaluza- Klein construction for a charged particle in a prescribed magnetic field. 

Let Qkk be the manifold S0{3) xM^ with variables (R, q). On Qkk introduce the Kaluza-Klein 
Lagrangian Lkk ■ TQkk - TS0{3) x TM^ i-^ M as 

LkkCR, q, R, q; z) = LKKi^, T, q, q) = ^ ( m , ) + ^|r + q|2 , (23.11) 

with ft = (R~^R)' and F = R~^z;. The Lagrangian Lkk is positive definite in (fi,r,q); so it may 
be regarded as the kinetic energy of a metric, the Kaluza-Klein metric on TQxk- 
The Legendre transformation for Lkk gives the momenta 

n = m and p = r + q. (23.12) 

Since Lkk does not depend on q, the Euler-Lagrange equation 

d OLkk _ QLkk _ Q 
dt 9q 9q ' 

shows that p = dLKK/dc{ is conserved. The constant vector p is now identified as the vector in 
the body, 

p = r + q = - mgx- 

After this identification, the heavy top action principle in Proposition 23.1 with the Kaluza-Klein 
Lagrangian returns Euler's motion equation for the heavy top (23.1). 

The Hamiltonian Hkk associated to Lkk by the Legendre transformation (23.12) is 

HKK{n,T,(i,p) = n-0 + p-q-LKi^(n,r,q,q) 

= ^n-r^n-p-r+^lpp 
= ln.i-in + ^|p-rp-l|r|2. 

Recall that F is a unit vector. On the constant level set |Fp = 1, the Kaluza-Klein Hamiltonian 
Hkk is a positive quadratic function, shifted by a constant. Likewise, on the constant level set 
p = —mgx, the Kaluza-Klein Hamiltonian Hkk is a function of only the variables (n,F) and is 
equal to the Hamiltonian (23.5) for the heavy top up to an additive constant. Consequently, the 
Lie-Poisson equations for the Kaluza-Klein Hamiltonian Hkk now reproduce Euler's motion equation 
for the heavy top (23.1). 

Exercise 23.5 Write the Kaluza-Klein construction on SE{3) = S'0(3)(S)M^. 

24 Euler-Poincare (EP) reduction theorem 

Remark 24.1 (Geodesic motion) As emphasized by [Arl966], in many interesting cases, the Euler- 
Poincare equations on the dual of a Lie algebra g* correspond to geodesic motion on the corresponding 
group G. The relationship between the equations on g* and on G is the content of the basic Euler- 
Poincare theorem discussed later. Similarly, on the Hamiltonian side, the preceding paragraphs de- 
scribed the relation between the Hamiltonian equations on T*G and the Lie-Poisson equations on g*. 
The issue of geodesic motion is especially simple: if either the Lagrangian on g or the Hamiltonian on 
g* is purely quadratic, then the corresponding motion on the group is geodesic motion. 
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24.1 We were already speaking prose (EP) 

Many of our previous considerations may be recast immediately as Euler-Poincare equations. 

• Rigid bodies ~ (^EP S'O(n)) , 

• Deforming bodies ~^EPGL_|_(n, M)j , 

• Heavy tops ~ (f,PSO{3) x M^) , 

• EPDiff 

24.2 Euler-Poinccire Reduction 

This lecture applies reduction by symmetry to Hamilton's principle. For a G— invariant Lagrangian 

defined on TG, this reduction takes Hamilton's principle from TG to TG/G ~ g. Stationarity of the 
symmetry-reduced Hamilton's principle yields the Euler-Poincare equations on g*. The corresponding 
reduced Legendre transformation yields the Lie-Poisson Hamiltonian formulation of these equations. 

Euler-Poincare Reduction starts with a right (respectively, left) invariant Lagrangian L : 
TG M on the tangent bundle of a Lie group G. This means that L(TgRh{v)) = L{v), respec- 
tively LiTgLhiy)) = L{v), for all g,h e G and all v G TgG. In shorter notation, right invariance of the 
Lagrangian may be written as 

L{g{t),g{t)) = L{9m,g{t)h), 

for all heG. 

Theorem 24.2 (Euler-Poincare Reduction) Let G be a Lie group, L : TG — > M a right-invariant 
Lagrangian, and I := L\g : q be its restriction to g. For a curve g{t) G G, let 

m = m ■ 9{tr^ ■= Tg^t^Rg^t^igit) e g . 

Then the following four statements are equivalent: 

(i) g{t) satisfies the Euler-Lagrange equations for Lagrangian L defined on G. 

(ii) The variational principle 

5 [ L{g{t),g{t))dt = 

J a 

holds, for variations with fixed endpoints. 

(iii) The (right invariant) Euler-Poincare equations hold: 

d SI _ SI 

dts^- "^U^- 

(iv) The variational principle 

S [ lim)dt = 

J a 

holds on Q, using variations of the form S^ = rj— ij], where r]{t) is an arbitrary path in g which 
vanishes at the endpoints, i.e., r}{a) = r]{b) = 0. 

Proof. 
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Step I. Proof that (i) <^=^ (ii): This is Hamilton's principle: the Euler-Lagrange equations follow 
from stationary action for variations 5g which vanish at the endpoints. (See Lecture 10.) 

Step II. Proof that (ii) <;=^ (iv): Proving equivalence of the variational principles (ii) on TG 
and (iv) on g for a right-invariant Lagrangian requires calculation of the variations of ^ = gg~^ 
induced by 5g. To simplify the exposition, the calculation will be done first for matrix Lie groups, 
then generalized to arbitrary Lie groups. 

Step IIA. Proof that (ii) <^=^ (iv) for a matrix Lie group. For ^ = gg~^, define ge{t) to be a 
family of curves in G such that go{t) = g{t) and denote 



5g := 

The variation of ^ is computed in terms of 5g as 



dg,{t) 



de 



e=0 



d_ 
de 



{9eg:^) 



<fg 



e=0 



dtde 



- 99-\Sg)g-K 



(24.1) 



e=0 



Set 7] := g ^Sg. That is, r/(t) is an arbitrary curve in q which vanishes at the endpoints. The time 
derivative of r/ is computed as 



r] 



drj d 
di ~ di 



d 
d~e 



9e 9 



e=0 



d'g 



dtde 



9-' - mg-^gg-^ . 



(24.2) 



€=0 



Taking the difference of (24.1) and (24.2) implies 

- = -99''^iSg)g-^ + {5g)g'^gg-^ = - ^j] + r]^ = - %r]] . 
That is, for matrix Lie algebras, 

where rj] is the matrix commutator. Next, we notice that right invariance of L allows one to change 
variables in the Lagrangian by applying g~^{t) from the right, as 

L{g{t),g{t)) = L{e,g{t)g-\t)) =: Z(^(t)) . 

Combining this definition of the symmetry-reduced Lagrangian Z : g — ^ M together with the formula 
for variations S( just deduced proves the equivalence of (ii) and (iv) for matrix Lie groups. 

Step IIB. Proof that (ii) <^=^ (iv) for an arbitrary Lie group. The same proof extends to 
any Lie group G by using the following lemma. 



Lemma 24.3 Let g : U C 



G be a smooth map and denote its partial derivatives by 



<:(+ .^ . T 7? dg{t,e) ,^ . T r %(^^ ^) 



Then 



(24.3) 

(24.4) 



5^ df] 
de~ dt ' 

where rj] is the Lie algebra bracket on q. Conversely, ifU C is simply connected and ^,r] : U q 
are smooth functions satisfying (24.4), then there exists a smooth function g : U ^ G such that (24.3) 
holds. 
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Proof. Write ^ = gg ^ and r] = g'g ^ in natural notation and express the partial derivatives 
g = dg/dt and g' = dg/de using the right translations as 

g = ^o g and g' = r]o g. 

By the chain rule, these definitions have mixed partial derivatives 

g' = = V^-ri and g' = ri = Vri-^. 

The difference of the mixed partial derivatives implies the desired formula (24.4), 

^' - ?7 = • ?7 - Vt? • ^ = - 77] = - ad^r? . 

(Note the minus sign in the last two terms.) ■ 

Step III. Proof of equivalence (Hi) <^=^ (iv) Let us show that the reduced variational principle 
produces the Euler-Poincare equations. We write the functional derivative of the reduced action 
Sred = JaKOdt with Lagrangian Z(^) in terms of the natural pairing (• , •) between q* and g as 

The last equality follows from integration by parts and vanishing of the variation •q{t) at the endpoints. 
Thus, stationarity 5 J^l{^{t))dt = for any ri{t) that vanishes at the endpoints is equivalent to 

d SI _ 61 



which are the Euler-Poincare equations. 



Remark 24.4 (Left-invariant Euler-Poincare equations) The same theorem holds for left in- 
variant Lagrangians on TG, except for a sign in the Euler-Poincare equations, 

d 51 _ 51 
JtTr^^^W 

which arises because left-invariant variations satisfy 6^ = rj + r]] (with the opposite sign). 

Exercise 24.5 Write out the corresponding proof of the Euler-Poincare reduction theorem for left- 
invariant Lagrangians defined on the tangent space TG of a group G. 

Reconstruction. The procedure for reconstructing the solution v{t) G Tg(j.^G of the Euler-Lagrange 
equations with initial conditions g{G) = go and g(0) = vq starting from the solution of the Euler- 
Poincare equations is as follows. First, solve the initial value problem for the right-invariant Euler- 
Poincare equations: 

d 51 ,^ 51 . . _i 

= - ad^^ with ^(0) = ^0 := vogo ■ 
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Then from the solution for ^(i) reconstruct the curve g{t) on the group by solving the "linear differential 
equation with time-dependent coefficients" 

9{t) = ^{t)g{t) with g{0)=go. 

The Eulcr-Poincare reduction theorem guarantees then that v{t) = g{t) = ^{t) ■ g{t) is a solution of 
the Euler-Lagrange equations with initial condition vq = ^ogo- 

Remark 24.6 A similar statement holds, with obvious changes for left-invariant Lagrangian systems 
on TG. 



24.3 Reduced Legendre transformation 

As in the equivalence relation between the Lagrangian and Hamiltonian formulations discussed earlier, 
the relationship between symmetry-reduced Euler-Poincare and Lie-Poisson formulations is determined 
by the Legendre transformation. 

Definition 24.7 The Legendre transformation ¥1 : q ^ q* is defined by 

Lie-Poisson Hamiltonian formulation. Let ^(/x) := — l{^)- Assuming that ¥1 is a diffeo- 

morphism yields 

tj./ '>(\ c 



So the Euler-Poincare equations for I are equivalent to the Lie-Poisson equations for h: 

d {5l\ ^Jl 

dt[5-d=~'"^'6-^^^ = -'^''^/'^^- 
The Lie-Poisson equations may be written in the Poisson bracket form 

f = {f,h}, (24.5) 

where / : g* — > M is an arbitrary smooth function and the bracket is the (right) Lie-Poisson bracket 
given by 



{f,h}{^i) = {^, 



6f Sh 



In the important case when i is quadratic, the Lagrangian L is the quadratic form associated to 
a right invariant Riemannian metric on G. In this case, the Euler-Lagrange equations for L on G 
describe geodesic motion relative to this metric and these geodesies are then equivalently described 
by either the Euler-Poincare, or the Lie-Poisson equations. 

Exercise 24.8 Exercise 22. 6 requires an extension of the pure EP reduction theorem for a Lagrangian 
L : {TG X TQ) — > M. Following the proof of the EP reduction theorem, make this extension. 



Exercise 24.9 Compute the pure EP equations for geodesic motion on SE{3) . These equations turn 
out to be applicable to the motion of an ellipsoidal body through a fluid. 
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25 EPDiff: the Euler-Poincare equation on the diffeomorphisms 
25.1 The n— dimensional EPDifF equation and its properties 

Eulerian geodesic motion of a fluid in n— dimensions is generated as an EP equation via Hamilton's 
principle, when the Lagrangian is given by the kinetic energy. The kinetic energy defines a norm ||u|p 
for the Eulerian fluid velocity, taken as u(x, t) : i?" x ^ K^. The choice of the kinetic energy as a 
positive functional of fluid velocity u is a modeling step that depends upon the physics of the problem 
being studied. We shall choose the Lagrangian, 

2 _ / „ /O „ .in / „ ™ jn 



||uf = J u-QopUd''x = J u md^'x, (25.1) 

so that the positive-definite, symmetric, operator Qop defines the norm ||u||, for appropriate (homo- 
geneous, say, or periodic) boundary conditions. The EPDiff equation is the Euler-Poincare equation 
for this Eulerian geodesic motion of a fluid. Namely, 

S^+'<l^=°- 4ul = ^l|uP. (25.2) 

Here ad* is the dual of the vector-fleld ad-operation (the commutator) under the natural pairing 
(• , •) induced by the variational derivative S£[u] = {61/ 6u , 6u). This pairing provides the definition 
of ad*, 

(ad* m , v) = — (m , aduv) , (25.3) 

where u and v are vector fields, aduV = [u,v] is the commutator, i.e., the Lie bracket given in 
components by (summing on repeated indices) 

^"'^^'^^^Ox^^^'^aa;^' °^ [u,v] = u- Vv-v- Vu. (25.4) 

The notation adu v := [u, v] formally denotes the adjoint action of the right Lie algebra of Diff(P) 
on itself, and m. = 51/ 5u is the fluid momentum, a one- form density whose co- vector components 
are also denoted as m. 

If u = uW/dx^ , m = rriidx^^dV, then the preceding formula for ad* (m(g)(iy) has the coordinate 
expression in M", 

(\ ■ f d ■ du'^ \ 

ad* mj dx* (g) = l-^{u^mi) + mj-^jdx'0dV. (25.5) 

In this notation, the abstract EPDiff equation (25.2) may be written explicitly in Euclidean coor- 
dinates as a partial differential equation for a co- vector function m(x, t) : i?" x i?-*^ — i?". Namely, 



— m + u • Vm -I- Vu"^ • m + m(divu) =0, with m= — = QopU. (25.6) 
ot ^ ' ' ^ ' ou 

Convection Stretching Expansion 



To explain the terms in underbraces, we rewrite EPDiff as preservation of the one-form density of 
momentum along the characteristic curves of the velocity. Namely, 

f m • dx (g) dv] = along ^ = u = G * m . (25.7) 
dt\ / dt 

This form of the EPDiff equation also emphasizes its nonlocality, since the velocity is obtained from 
the momentum density by convolution against the Green's function G of the operator Qop. Thus, 
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u = G * m with QopG = 6{x), the Dirac measure. We may check that this "characteristic form" of 
EPDiff recovers its Eulerian form by computing directly, 

fm • dx (g) dv) = ^ • dx (g) + m • (g) dF + m • dx (g) (^dv\ along ^ = u = G * m 
dt\ J dt dt \dt / dt 

= + u • Vm + Vu^ ■ m + m(div u)^ • dx (g = . 

Exercise 25.1 Show that EPDiff may be written as 

+ -Cu) (m • dx ® dy) = , (25.8) 

where Cu is the Lie derivative with respect to the vector field with components u = G*in. Hint: How 
does this Lie-derivative form of EPDiff in (25.8) differ from its characteristic form (25.7)9 

EPDiff may also be written equivalently in terms of the operators div, grad and curl in 2D and 
3D as, 

d 

— m — u X curl m + V(u • m) + m(div u) = . (25.9) 

Thus, for example, its numerical solution would require an algorithm which has the capability to deal 
with the distinctions and relationships among the operators div, grad and curl. 



25.2 Derivation of the n— dimensional EPDiff equation as geodesic flow 

Let's derive the EPDiff equation (25.6) by following the proof of the EP reduction theorem leading to 
the Euler-Poincare equations for right invariance in the form (25.2). Following this calculation for the 
present case yields 

aduV^ dt 

where (■ , •) is the pairing between elements of the Lie algebra and its dual. In our case, this is the 
pairing, e.g.. 

This pairing allows us to compute the coordinate form of the EPDiff equation explicitly, as 

f ( d SI d ( 51 A 51 dun 

Substituting m. = Sl/du now recovers the coordinate forms for the coadjoint action of vector fields in 
(25.5) and the EPDiff equation itself in (25.6). When £[u] = EPDiff describes geodesic motion 

on the diffeomorphisms with respect to the norm ||u||. 



J a 



)dt = 
^ I 51 



5l_ 
~5^ 



Su ) dt 



51 



5l_ 



^ , aduV ) dt 
Su / 
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Lemma 25.2 In Step IIB of the proof of the Euler-Poincare reduction theorem that (ii) <^=^ (iv) 
for an arbitrary Lie group, a certain formula for the variations for time-dependent vector fields was 
employed. That formula was employed again in the calculation above as, 

<5u = v-aduV. (25.10) 

This formula may be rederived as follows in the present context. We write u = gg^^ and v = g'g~^ 
in natural notation and express the partial derivatives g = dg/dt and g' = dg/de using the right 
translations as 

g = uo g and g' = v o g . 

To compute the mixed partials, consider the chain rule for say u{g{t, e)xo) and set x(t, e) = g{t, e) • xq. 
Then, 

, du dx. du . 5u , _i du , , 

The chain rule for v gives a similar formula with u and v exchanged. Thus, the chain rule gives two 
expressions for the mixed partial derivative g' as 

g' = VL = Vu • V and g' = v = Vv • u . 

The difference of the mixed partial derivatives then implies the desired formula (25.10), since 

u' — V = Vu • V — Vv • u = — [u, v] = — aduV . 

26 EPDiff: the Euler-Poincare equation on the diffeomorphisms 

In this lecture, we shall discuss the solutions of EPDiff for prcssureless compressible geodesic motion 
in one spatial dimension. This is the EPDiff equation in ID,^ 

dtm + ad*m = , or, equivalently, (26.1) 
dtm + urux + 2uxm = , with m = QopU ■ (26.2) 

• The EPDiff equation describes geodesic motion on the diffeomorphism group with respect to a 
family of metrics for the fluid velocity u{t,x), with notation, 

S£ 

m = — = QopU for a kinetic-energy Lagrangian (26.3) 
du 

£{u) = - uQopudx=-\\uf. (26.4) 

• In one dimension, Qop in equation (26.3) is a positive, symmetric operator that defines the kinetic 
energy metric for the velocity. 

• The EPDiff equation (26.2) is written in terms of the variable m = Si/Su. It is appropriate to 
call this variational derivative m, because it is the momentum density associated with the fluid 
velocity u. 

®A one-form density in ID takes the form m {dxf and the EP equation is given by 

'^(midx)^) = ^^(dx)^ -\-2rn{du)(dx) = Q with dx = du = Uxdx and u = G*m, 
dt dt dt 

where G * m denotes convolution with a function G on the real line. 
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• Physically, the first nonlinear term in the EPDiff equation (26.2) is fluid transport. 

• The coefficient 2 arises in the second nonlinear term, because, in one dimension, two of the 
summands in ad* m = unix + 2uxm are the same, cf. equation (25.5). 

• The momentum is expressed in terms of the velocity by m = Si/5u = QopU- Equivalently, for 
solutions that vanish at spatial inflnity, one may think of the velocity as being obtained from 
the convolution, 

u{x) = G * m{x) = J G{x — y)m{y) dy , (26.5) 
where G is the Green's function for the operator Qop on the real line. 

• The operator Qop and its Green's function G are chosen to be even under reflection, G{—x) = 
G{x), so that u and m have the same parity. Moreover, the EPDiff equation (26.2) conserves 
the total momentum M = J m{y) dy, for any even Green's function. 

Exercise 26.1 Show that equation (26.2) conserves M = J m{y)dy for any even Green's func- 
tion G{—x) = G{x), for either periodic, or homogeneous boundary conditions. 



• The traveling wave solutions of ID EPDiff when the Green's function G is chosen to be even 
under reflection are the "pulsons," 

u{x,t) = cG{x — ct) . 

Exercise 26.2 Prove this statement, that the traveling wave solutions of ID EPDiff are pulsons 
when the Green's function is even. What role is played in the solution by the Green's function 
being even? Hint: Evaluate the derivative of an even function at x = 0. 

• See Fringer and Holm FrHo2001 and references therein for further discussions and numerical 
simulations of the pulson solutions of the ID EPDiff equation. 



26.1 Pulsons 

The EPDiff equation (26.2) on the real line has the remarkable property that its solutions collec- 
tivize^^ into the flnite dimensional solutions of the "iV— pulson" form that was discovered for a special 
form of G in Camassa and Holm CaHol993, then was extended for any even G in Pringer and Holm 
FrHo2001, 

N 

u{x,t)=Y,Pi{t)G{x-qi{t)). (26.6) 

i=l 

Since G{x) is the Green's function for the operator Qop, the corresponding solution for the momentum 
m = QopU is given by a sum of delta functions, 

N 

m{x,t) = Y,Pi{t)^{x-<li{t))- (26.7) 

i=l 



^"Sce [GuStl984] for discussions of the concept of collective variables for Hamiltonian theories. We will discuss the 
collectivization for the EPDiff equation later from the viewpoint of momentum maps. 
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Thus, the time-dependent "coUective coordinates" qi{t) and pi{t) are the positions and velocities of the 
N pulses in this solution. These parameters satisfy the finite dimensional geodesic motion equations 
obtained as canonical Hamiltonian equations 

= ^ = J2PjGiQi-Qj), (26.8) 

Pi = - = -Pi^ Pj G'iQi- Qj) , (26.9) 

in which the Hamiltonian is given by the quadratic form, 

1 ^ 

Hn=^J2 PiPj ^("li - • (26.10) 

Remark 26.3 In a certain sense, equations (26.8-26.9) comprise the analog for the peakon momen- 
tum relation (26.7) of the "symmetric generalized rigid body equations" in (20.1). 

Thus, the canonical equations for the Hamiltonian Hj^j describe the nonlinear collective interactions 
of the A^— pulson solutions of the EPDiff equation (26.2) as finite-dimensional geodesic motion of a 
particle on an A^— dimensional surface whose co-metric is 

G'^{q)=G{qi-qj). (26.11) 

Pringer and Holm PrHo2001 showed numerically that the A?^— pulson solutions describe the emergent 
patterns in the solution of the initial value problem for EPDiff equation (26.2) with spatially confined 
initial conditions. 

Exercise 26.4 Equations (26.8-26.9) describe geodesic motion. 

1. Write the Lagrangian and Euler- Lagrange equations for this motion. 

2. Solve equations (26.8-26.9) for N = 2 when lim|j.|_>3o G{x) = 0. 

(a) Why should the solution be described as exchange of momentum in elastic collisions? 

(b) Consider both head-on and overtaking collisions. 

(c) Consider the antisymmetric case, when the total momentum vanishes. 



Integrability Calogero and Francoise [Cal995], [CaFrl996] found that for any finite number N 
the Hamiltonian equations for Hn in (26.10) are completely integrable in the Liouville sense^^ for 
G = Gi{x) = A -|- lucos(^'x) -|- jUi sin(i/|x|) and G = G2{x) = a-\- P\x\ -\- 70;^, with A, fi, /ii, v, and a, P, 
7 being arbitrary constants, such that A and fi are real and and u both real or both imaginary.^^ 
Particular cases of Gi and G2 are the peakons Gi{x) = e"'^'/" of [CaHol993] and the compactons 
G2{x) = max(l — 0) of the Hunter-Saxton equation, [HuZhl994]. The latter is the EPDiff equation 
(26.2), with i{u) = \ j u^dx and thus m = —Uxx- 

^^A Hamiltonian system is integrable in the Liouville sense, if the number of independent constants of motion in 
involution is the same as the number of its degrees of freedom. 
^^This choice of the constants keeps Hn real in (26.10). 
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Lie-Poisson Hamiltonian form of EPDiff In terms of m, the conserved energy Hamiltonian for 
the EPDiff equation (26.2) is obtained by Legendre transforming the kinetic energy Lagrangian, as 



Thus, the Hamiltonian depends on m, as 



h{m) 



1 



m{x)G{x — y)m{y) dxdy , 



which also reveals the geodesic nature of the EPDiff equation (26.2) and the role of G{x) in the kinetic 
energy metric on the Hamiltonian side. 

The corresponding Lie-Poisson bracket for EPDiff as a Hamiltonian evolution equation is given 

by, 

6h , 6h 



dtm = {m, = — {dm, + md) 



- — and - — = u , 
dm dm 



which recovers the starting equation and indicates some of its connections with fluid equations on the 
Hamiltonian side. For any two smooth functionals f,h of m in the space for which the solutions of 
EPDiff exist, this Lie-Poisson bracket may be expressed as. 



m 



5l_ Sh^ 
6m ' Sm 



dx 



where [• , •] denotes the Lie algebra bracket of vector fields. That is. 



Sf^ 5h_ 
dm ' dm 



Sm Sm Sm Sm 



Exercise 26.5 What is the Casimir for this Lie Poisson bracket? What does it mean from the view- 
point of coadjoint orbits ? 



26.2 Peakons 

The case G{x) = e~l^l/" with a constant lengthscalc a is the Green's function for which the operator 
in the kinetic energy Lagrangian (26.3) is Qop = 1 — ct^d^- For this (Helmholtz) operator Qgp, the 
Lagrangian and corresponding kinetic energy norm are given by, 

= :^ = t: I u QonU dx = - u + a u„ dx , for lim u = . 

2" " 2 J ^ 2 J ^ \x\^oo 

This Lagrangian is the H ^ norm of the velocity in one dimension. In this case, the EPDiff equation 
(26.2) is also the zero-dispersion limit of the completely integrable CH equation for unidirectional 
shallow water waves first derived in Camassa and Holm CaHol993, 

mt + umx + 2mUx = —cqUx + JUxxx > m = u — ct^Uxx ■ (26.12) 

This equation describes shallow water dynamics as completely integrable soliton motion at quadratic 
order in the asymptotic expansion for unidirectional shallow water waves on a free surface under 
gravity. See Dullin, Gottwald and Holm DGH[2001, 2003, 2004] for more details and explanations of 
this asymptotic expansion for unidirectional shallow water waves to quadratic order. 



D. D. Holm 



Imperial College London 



Applications of Poisson Geometry 



107 



Because of the relation m = u — a^Uxx, equation (26.12) is nonlocal. In other words, it is an 
integral-partial differential equation. In fact, after writing equation (26.12) in the equivalent form, 

2 

(1 - a^d^){ut + uUx) = - div?' + y-^x) - cq'^x + lu^xx , (26.13) 
one sees the interplay between local and nonlocal linear dispersion in its phase velocity relation. 



u cq — 'J k'^ 



k 1 + q2A;2' 



(26.14) 



for waves with frequency lo and wave number k linearized around u = 0. For 7/co < 0, short waves 
and long waves travel in the same direction. Long waves travel faster than short ones (as required in 
shallow water) provided 7/co > — a^. Then the phase velocity lies in the interval uj/k G (— 7/a^,co]. 
The famous Korteweg-de Vries (KdV) soliton equation, 

Ut + SuUx = -CqUx + ^Uxxx , (26.15) 

emerges at linear order in the asymptotic expansion for shallow water waves, in which one takes ^ 
in (26.13) and (26.14). In KdV, the parameters cq and 7 are seen as deformations of the Riemann 
equation, 

Ut + 3uux = . 

The parameters cq and 7 represent linear wave dispersion, which modifies and eventually balances the 

tendency for nonlinear waves to steepen and break. The parameter a, which introduces nonlocality, 
also regularizes this nonlinear tendency, even in the absence of cq and 7. 



27 DifFeons — singular momentum solutions of the EPDifF equation 
for geodesic motion in higher dimensions 

As an example of the EP theory in higher dimensions, wc shall generalize the one-dimensional pulson 
solutions of the previous section to ra— dimensions. The corresponding singular momentum solutions 
of the EPDiff equation in higher dimensions are called "diffeons." 



27.1 n— dimensional EPDifF equation 

Eulerian geodesic motion of a fluid in n— dimensions is generated as an EP equation via Hamilton's 
principle, when the Lagrangian is given by the kinetic energy. The kinetic energy defines a norm ||u|p 
for the Eulerian fluid velocity, u(x, t) : R"^ x i?^ — > R^. As mentioned earlier, the choice of the kinetic 
energy as a positive functional of fluid velocity u is a modeling step that depends upon the physics 
of the problem being studied. Following our earlier procedure, as in equations (25.1) and (25.2), we 
shall choose the Lagrangian, 

||uf = j u-Qopud''x = Ju-md^-x, (27.1) 

so that the positive-definite, symmetric, operator Qop defines the norm ||u||, for appropriate boundary 
conditions and the EPDiff equation for Eulerian geodesic motion of a fluid emerges. 



(27.2) 
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Legendre transforming to the Hamiltonian side The corresponding Legendre transform yields 
the following invertible relations between momentum and velocity, 

m = <3opU and u = G * m , (27.3) 

where G is the Green's function for the operator Qop, assuming appropriate boundary conditions 
(on u) that allow inversion of the operator Qop to determine u from m. 
The corresponding Hamiltonian is, 

^[m] = (m, u) - ^||uf = ^ /"m-G*m J||mf , (27.4) 

Z Z J z 

which also defines a norm ||m|| via a convolution kernel G that is symmetric and positive, when the 
Lagrangian £[u] is a norm. As expected, the norm ||m|| given by the Hamiltonian h[m] specifies the 
velocity u in terms of its Legendre-dual momentum m by the variational operation. 



_ 6h^ 
Sm 



G * m = y G(x - y) m(y) d "j/ . (27.5) 



We shall choose the kernel G(x — y) to be translation-invariant (so Noether's theorem implies that 
total momentum M = J m d is conserved) and symmetric under spatial reflections (so that u and 
m have the same parity). 

After the Legendre transformation (27.4), the EPDiff equation (27.2) appears in its equivalent 
Lie-Poisson Hamiltonian form, 

d 

— m = {m,/t} = -ad|^/5„m. (27.6) 

Here the operation {• , • } denotes the Lie-Poisson bracket dual to the (right) action of vector fields 
amongst themselves by vector-field commutation 



{/,^} = -(m 



5f 5h' 
Sm ' 5m 

For more details and additional background concerning the relation of classical EP theory to Lie- 
Poisson Hamiltonian equations, see [MaRal994, HoMaRal998a]. 

In a moment we will also consider the momentum maps for EPDiff. 



27.2 DifFeons: n— dimensional Analogs of Pulsons for the EPDiff equation 

The momentum for the one-dimensional pulson solutions (??) on the real line is supported at points 
via the Dirac delta measures in its solution ansatz, 

N 

m{x,t) = ^Pi{t)6{x - qi{t)) , meR^. (27.7) 

i=l 

We shall develop n— dimensional analogs of these one-dimensional pulson solutions for the Euler- 
Poincare equation (25.9) by generalizing this solution ansatz to allow measure-valued n— dimensional 
vector solutions m G i?" for which the Euler-Poincare momentum is supported on co-dimension— /c 
subspaces R^~^ with integer k G [1,?^]- For example, one may consider the two-dimensional vector 
momentum m G i?^ in the plane that is supported on one-dimensional curves (momentum fronts). 
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Likewise, in three dimensions, one could consider two-dimensional momentum surfaces (sheets), one- 
dimensional momentum filaments, etc. The corresponding vector momentum ansatz that we shall use 
is the following, cf. the pulson solutions (27.7), 

N 

m(x,t) = ^ / Pi(s,t)(5(x-Qj(s,t))ds, m G i?" . (27.8) 

i=l 

Here, Pt,Qi G i?" for i = l,2,...,iV. For example, when n — k = 1, so that s G i?^ is one- 
dimensional, the delta function in solution (27.8) supports an evolving family of vector-valued curves, 
called momentum filaments. (For simplicity of notation, we suppress the implied subscript i in 
the arclength s for each Pj and Qj.) The Legendrc-dual relations (27.3) imply that the velocity 
corresponding to the momentum filament ansatz (27.8) is, 

N 

u(x,i) = G*m = ^ Pj{s',t)G{x-Qj{s',t))ds' . (27.9) 

Just as for the ID case of the pulsons, wc shall show that substitution of the n— D solution ansatz 
(27.8) and (27.9) into the EPDiff equation (25.6) produces canonical geodesic Hamiltonian equations 
for the n— dimensional vector parameters Qj(s, t) and Pi(s, t), i = 1,2, . . . ,N. 



27.2.1 Canonical Hamiltonian dynamics of difFeon momentum filaments in 

For definiteness in what follows, wc shall consider the example of momentum filaments m G 
supported on one-dimensional space curves in i?", so s G i?^ is the arclength parameter of one of these 
curves. This solution ansatz is reminiscent of the Biot-Savart Law for vortex filaments, although the 
flow is not incompressible. The dynamics of momentum surfaces, for s e R'^ with k < n, follow a 
similar analysis. 

Substituting the momentum filament ansatz (27.8) for s G R^ and its corresponding velocity (27.9) 
into the Euler-Poincare equation (25.6), then integrating against a smooth test function (/>(x) implies 
the following canonical equations (denoting explicit summation on j G 1, 2, . . . TV), 



d_ 
dt 



Gli{s,t) =Y,jv^{s',t)G{Qli{s,t)-Qlj{s',t)))ds' 

(27.10) 



5Pz 

N 



^^i{s,t) = -E/ {^i{s,t)-P^{s',t))-^^^^G{Cli{s,t)-Cl^{s',t))ds' 

, (sum on _7, no sum on z) . (27.11) 



The dot product Pj • Pj denotes the inner, or scalar, product of the two vectors Pj and Pj in i?". Thus, 
the solution ansatz (27.8) yields a closed set of integro-partial-differential equations (IPDEs) 

given by (27.10) and (27.11) for the vector parameters Qj(s, t) and Pi(s, t) with i = 1,2 . . . N . These 
equations are generated canonically by the following Hamiltonian function '■ {R"' x i^")®^ 



N 

^^=\jjll {^i{s,t)-Pj{s',t))G(Cti{s,t)-Ct^{s',t))dsds' . (27.12) 
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This Hamiltonian arises by substituting the momentum ansatz (27.8) into the Hamiltonian (27.4) 
obtained from the Legendre transformation of the Lagrangian corresponding to the kinetic energy norm 
of the fluid velocity. Thus, the evolutionary IPDE system (27.10) and (27.11) represents canonically 
Hamiltonian geodesic motion on the space of curves in with respect to the co-metric given on these 
curves m (27.12). The Hamiltonian = in (27.12) defines the norm ||P|| in terms of this 

co-metric that combines convolution using the Green's function G and sum over filaments with the 
scalar product of momentum vectors in i?". 

Remark 27.1 Note the Lagrangian property of the s coordinate, since 

^Qi{s,t) = u{Qi{s,t),t). 

28 Singular solution momentum map Jsmg for diffeons 

The diffcon momentum filament ansatz (27.8) reduces, and collectivizes the solution of the geodesic 
EP PDE (25.6) in n + 1 dimensions into the system (27.10) and (27.11) of 2N canonical evolutionary 
IPDEs. One can summarize the mechanism by which this process occurs, by saying that the map 
that implements the canonical (Q, P) variables in terms of singular solutions is a (cotangent bundle) 
momentum map. Such momentum maps are Poisson maps; so the canonical Hamiltonian nature of 
the dynamical equations for (Q,P) fits into a general theory which also provides a framework for 
suggesting other avenues of investigation. 

Theorem 28.1 The momentum a,nsatz (27.8) for measure-valued solutions of the EPDiff equation 
(25.6), defines an equivariant momentum map 

Jsing : r*Emb(S,M") ^ X(M")* 

that is called the singular solution momentum map in [HoMa2004]- 

We shall explain the notation used in the theorem's statement in the course of its proof. Right 
away, however, we note that the sense of "defines" is that the momentum solution ansatz (27.8) 

expressing m (a vector function of spatial position x) in terms of Q, P (which are functions of s) can 
be regarded as a map from the space of (Q(s), P(s)) to the space of m's. This will turn out to be the 
Lagrange-to-Euler map for the fluid description of the singular solutions. 

Following [HoMa2004], we shall give two proofs of this result from two rather different viewpoints. 
The first proof below uses the formula for a momentum map for a cotangent lifted action, while the 
second proof focuses on a Poisson bracket computation. Each proof also explains the context in which 
one has a momentum map. (See [MaRal994] for general background on momentum maps.) 

First Proof. For simplicity and without loss of generality, let us take = 1 and so suppress the 
index a. That is, we shall take the case of an isolated singular solution. As the proof will show, this 
is not a real restriction. 

To set the notation, fix a /c-dimensional manifold S with a given volume element and whose points 

arc denoted s G S. Let Emb(S', M") denote the set of smooth cmbcddings Q : S — M". (If the 
EPDiff equations are taken on a manifold A4, replace M" with M.) Under appropriate technical 
conditions, which we shall just treat formally here, Emb(S', M**) is a smooth manifold. (See, for 
example, [EbMal970] and [MaHul983] for a discussion and references.) 

The tangent space rQEmb(S',M") to Emb(5,M") at the point Q G Emb(5,M") is given by the 
space of material velocity fields, namely the linear space of maps V : S ^ that are vector 
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fields over the map Q. The dual space to this space will be identified with the space of one-form 
densities over Q, which we shall regard as maps P : S — (M")*. In summary, the cotangent bundle 
T* Emb(S',M"') is identified with the space of pairs of maps (Q,P). 

These give us the domain space for the singular solution momentum map. Now we consider the 
action of the symmetry group. Consider the group (5 = Diff of diffeomorphisms of the space & in 
which the EPDiff equations are operating, concretely in our case M". Let it act on © by composition 
on the left. Namely for rj G Diff(M"), we let 

ri-Q = rioQ. (28.1) 

Now lift this action to the cotangent bundle T* Emb(5', R") in the standard way (see, for instance, 
MaRal994 for this construction). This lifted action is a symplectic (and hence Poisson) action and has 
an equivariant momentum map. We claim that this momentum map is precisely given by the ansatz 
(27.8). 

To see this, one only needs to recall and then apply the general formula for the momentum map 
associated with an action of a general Lie group on a configuration manifold Q and cotangent lifted 
to T*Q. 

First let us recall the general formula. Namely, the momentum map is the map 3 ■.T*Q ^ q* {g* 
denotes the dual of the Lie algebra g of 0) defined by 

JK)-e=K,CQ((Z)>, (28.2) 

where aq G T*Q and ^ G 0, where (^q is the infinitesimal generator of the action of © on Q associated 
to the Lie algebra element ^, and where {ceq,^Q{q)) is the natural pairing of an element of T*Q with 
an element oi TqQ. 

Now we apply this formula to the special case in which the group <S is the diffeomorphism group 
Diff(M"'), the manifold Q is Emb(S', M") and where the action of the group on Emb(S', M") is given 
by (28.1). The sense in which the Lie algebra of = Diff is the space q = X oi vector fields is well- 
understood. Hence, its dual is naturally regarded as the space of one-form densities. The momentum 
map is thus a map J : T* Emb(S',M") — X*. 

With J given by (28.2), we only need to work out this formula. First, we shall work out the 

infinitesimal generators. Let X G X be a Lie algebra element. By differentiating the action (28.1) 
with respect to r] in the direction of X at the identity element we find that the infinitesimal generator 
is given by 

-'^Emb(S,R")(Q) = -'^^ O Q- 

Thus, taking aq to be the cotangent vector (Q,P), equation (28.2) gives 

(J(Q,P),X) = ((Q,P),XoQ) 

= / Piis)X\Q{s))d''s. 
Js 

On the other hand, note that the right hand side of (27.8) (again with the index a suppressed, and 
with t suppressed as well), when paired with the Lie algebra element X is 

P(s)(5(x-Q(s))d'=s,X^ = j j (^Pi{s)S{yi-Q{s))d''s^X\yi)(r'x 

= / Pi{s)X'{Q{s)d''s. 
Js 



D. D. Holm 



Imperial College London 



Applications of Poisson Geometry 



112 



This shows that the expression given by (27.8) is equal to J and so the result is proved. ■ 

Second Proof. As is standard (see, for example, MaRal994), one can characterize momentum maps 
by means of the following relation, required to hold for all functions F on T* Emb(S', M"); that is, 
functions of Q and P: 

{F,{J,0} = ^P[F]. (28.3) 

In our case, we shall take J to be given by the solution ansatz and verify that it satisfies this relation. 
To do so, let ^ G j£ so that the left side of (28.3) becomes 

On the other hand, one can directly compute from the definitions that the infinitesimal generator 
of the action on the space T*Emb(S', R") corresponding to the vector field ^*(x)g^ (a Lie algebra 
element), is given by (see MaRal994, formula (12.1.14)): 

5Q = eoQ, 5P = -p^{s)^e{Q{s)), 

which verifies that (28.3) holds. 

An important element left out in this proof so far is that it does not make clear that the momentum 
map is equivariant, a condition needed for the momentum map to be Poisson. The first proof took 
care of this automatically since momentum maps for cotangent lifted actions are always 
equivariant and hence are Poisson. 

Thus, to complete the second proof, we need to check directly that the momentum map is equivari- 
ant. Actually, we shall only check that it is infinitesimally invariant by showing that it is a Poisson map 

from T* Emb(S', R") to the space of m's (the dual of the Lie algebra of X) with its Lie-Poisson bracket. 
This sort of approach to characterize equivariant momentum maps is discussed in an interesting way 
in [We2002]. 

The following direct computation shows that the singular solution momentum map (27.8) is Pois- 
son. This is accomplished by using the canonical Poisson brackets for {P}, {Q} and applying the 
chain rule to compute {mi(x), mj(y)}, with notation S'i^{y) = dd{y)/dy'^. We get 



SF 
5Q 



SF 6 
SPjSQ- 



{mi(x),mj(y)} 



N N 

Y^lds Ptis, t) 5(x - Q%s, *)) , E/ ds'P^is', t) 6{y - Q\s' , t)) 

a=l 6=1 



N 



a,fe=l 



E dsds' 



{Pt{s),P^{s')}6{^ - Q«(.)) S{y - Q\s')) 



- {P,ns),Qi{s')}P'l{s') 5(x - Q'^{s))6',{y - Q\s')) 



b( J\ 



- {Ql{s), P^{s')}Pf{s)5'^{^ - Q"(.)) 5{y - Ct\s')) 

+ {Ql{s), Q\{s')}P^{s)P^{s')5'^{^ - Ct{s)) 5',{y - (^\s')) 
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Substituting the canonical Poisson bracket relations 

{QUs),Q%s')} = 0, and 
{QUs),P^{s')} = S^%jS{s-s') 

into the preceding computation yields, 

{mi(x),mj(y)} 

. N N 

= dsPt{s,t)S{:K-Q'^{s,t)),J2 ds'P^{s',t)6{y-Q'{s',t)) 

^ a=l-' 6=1'^ 

N 

= ^ /d.Pf(.),5(x-Q»),5;(y-Q"(.)) 
0=1 

N 

dsPt{s)S'^{^-Q'^{s))6{y-Q^{s)) 

a=l 

Thus, 

{mi(x) , mj{y)} = - (mj(x)^ + ^ mi(x))<5(x - y) , (28.4) 

which is readily checked to be the Lie-Poisson bracket on the space of m's, restricted to their singular 
support. This completes the second proof of theorem. ■ 

Each of these proofs has shown the following basic fact. 

Corollary 28.2 The singular solution momentum map defined by the singular solution ansatz (27.8), 
namely, 

Jsing : r*Emb(5,M") ^ X(M")* 

is a Poisson map from the canonical Poisson structure on T* Emb(S', M") to the Lie-Poisson structure 
on 3e(R'")*. 

This is perhaps the most basic property of the singular solution momentum map. Some of its more 
sophisticated properties are outlined in [HoMa2004]. 

Pulling Back the Equations. Since the solution ansatz (27.8) has been shown in the preceding 

Corollary to be a Poisson map, the pull back of the Hamiltonian from X* to r*Emb(S', M") gives 
equations of motion on the latter space that project to the equations on X* . 

Thus, the basic fact that the momentum map Jsing is Poisson explains why the functions 
Q"(s, i) ond P"(s,t) satisfy canonical Hamiltonian equations. 

Note that the coordinate s G R*^ that labels these functions is a "Lagrangian coordinate" in the sense 
that it does not evolve in time but rather labels the solution. 
In terms of the pairing 

(•,•): 0*xg^M, (28.5) 
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between the Lie algebra g (vector fields in M") and its dual g* (one- form densities in M"), the following 
relation holds for measure- valued solutions under the momentum map (27.8), 

(m , u) = / m ■ u d ""x , pairing for m & u G M", 




= ((P,Q)), (28.6) 

which is the natural pairing between the points (Q.P) G T* Emb(S', R") and (Q, Q) G TEmb(S', M"). 
This corresponds to preservation of the action of the Lagrangian ^[u] under cotangent lift of Diff (M"). 

The pull-back of the Hamiltonian H[m\ defined on the dual of the Lie algebra g*, to T* Emb(S', M") 
is easily seen to be consistent with what we had before: 

H[in] ^ ^(m , G * m) = ^((P , G * P)) = Hm[P, Q] . (28.7) 

In summary, in concert with the Poisson nature of the singular solution momentum map, we see 
that the singular solutions in terms of Q and P satisfy Hamiltonian equations and also define an 
invariant solution set for the EPDiff equations. In fact, 

This invariant solution set is a special coadjoint orbit for the diffeomorphism group, as we 
shall discuss in the next section. 

29 The Geometry of the Momentum Map 

In this section we explore the geometry of the singular solution momentum map discussed earlier in a 
little more detail. The treatment is formal, in the sense that there are a number of technical issues in 
the infinite dimensional case that will be left open. We will mention a few of these as we proceed. 

29.1 Coadjoint Orbits. 

We claim that the image of the singular solution momentum map is a coadjoint orbit in X* . This 
means that (modulo some issues of connectedness and smoothness, which we do not consider here) 
the solution ansatz given by (27.8) defines a coadjoint orbit in the space of all one-form densities, 
regarded as the dual of the Lie algebra of the diffeomorphism group. These coadjoint orbits should 
be thought of as singular orbits — that is, due to their special nature, they are not generic. 

Recognizing them as coadjoint orbits is one way of gaining further insight into why the singular 
solutions form dynamically invariant sets — it is a general fact that coadjoint orbits in g* are symplectic 
submanifolds of the Lie-Poisson manifold g* (in our case j£(M"')*) and, correspondingly, are dynamically 
invariant for any Hamiltonian system on g*. 

The idea of the proof of our claim is simply this: whenever one has an equivariant momentum 
map 3 : P ^ Q* ioT the action of a group G on a symplectic or Poisson manifold P, and that action is 
transitive, then the image of J is an orbit (or at least a piece of an orbit). This general result, due to 
Kostant, is stated more precisely in [MaRal994], Theorem 14.4.5. Roughly speaking, the reason that 
transitivity holds in our case is because one can "move the images of the manifolds S around at will 
with arbitrary velocity fields" using diffeomorphisms of M". 
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29.2 The Momentum map J5 and the Kelvin circulation theorem. 

The momentum map Jsing involves Diff(M'^), the left action of the diffconiorphism group on the space 
of embeddings Emb(S', M") by smooth maps of the target space W\ namely, 

Difr(M") : Q-r7 = ?7oQ, (29.1) 



where, recall, Q : 5 — > M". As above, the cotangent bundle T* Emb(S',M'*) is identified with the space 
of pairs of maps (Q, P), with Q : S" ^ and P : 5" ^ r*R". 

However, there is another momentum map associated with the right action of the diffeomor- 
phism group of S on the embeddings Emb(S', M") by smooth maps of the "Lagrangian labels" S (fluid 
particle relabeling by rj : S ^ S). This action is given by 

Diff (5) : Q • ?7 = Q o . (29.2) 

The infinitesimal generator of this right action is 

^Emb(S,R") {Q) = f^^_Qor]t = TQoX. (29.3) 



t=0 



where X G X is tangent to the curve r/t at t = 0. Thus, again taking A'^ = 1 (so we suppress the 
index a) and also letting in the momentum map formula (28.2) be the cotangent vector (Q,P), 
one computes J^: 

{Js(Q,P),x) = {(Q,P),rQ.x) 

= J x(p{s)-dQ{s)^d''s 

= (^j ¥{s)-dCl{s)(S) d^s,X{s)^ 
= (P-dQ,X). 

Consequently, the momentum map formula (28.2) yields 

Js(Q,P) = P-dQ, (29.4) 

with the indicated pairing of the one-form density P • dQ with the vector field X. We have set things 
up so that the following is true. 

Proposition 29.1 The momentum map J5 is preserved by the evolution equations (27.10-27.11) for 
Q and P . 

Proof. It is enough to notice that the Hamiltonian in equation (27.12) is invariant under the 
cotangent lift of the action of Diff (5); it merely amounts to the invariance of the integral over S under 
reparametrization; that is, the change of variables formula; keep in mind that P includes a density 
factor. ■ 



Remark 29.2 
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• This result is similar to the Kelvin- Noether theorem for circulation T of an ideal fluid, which 
may be written as T = ^^^^^ D{s)~^P{s) ■ dQ{s) for each Lagrangian circuit c{s), where D is 
the mass density and P is again the canonical momentum density. This similarity should come 
as no surprise, because the Kelvin- Noether theorem for ideal fluids arises from invariance of 
Hamilton's principle under fluid parcel relabeling by the same right action of the diffeomorphism 
group, as in (29.2). 

• Note that, being an equivariant momentum map, the map J5, as with Jsing; is also a Poisson 
map. That is, substituting the canonical Poisson bracket into relation (29.4); that is, the relation 
M(x) = Pj(x)VQ'(x) yields the Lie-Poisson bracket on the space o/M's. We use the 
different notations m and M because these quantities are analogous to the body and spatial 
angular momentum for rigid body mechanics. In fact, the quantity m given by the solution 
Ansatz; specifically, m = JsinglQ^P) gives the singular solutions of the EPDiff equations, while 
M(x) = J5(Q,P) = -Pi(x)V(5*(x) is a conserved quantity. 

• In the language of fluid mechanics, the expression of m in terms of (Q,P) is an example of 
a Clebsch representation, which expresses the solution of the EPDiff equations in terms of 
canonical variables that evolve by standard canonical Hamilton equations. This has been known 
in the case of fluid mechanics for more than 100 years. For modern discussions of the Clebsch 
representation for ideal fluids, see, for example, [HoKul983, MaWel983]. 

• One more remark is in order; namely the special case in which S = M is of course allowed. In this 
case, Q corresponds to the map rj itself and P just corresponds to its conjugate momentum. The 
quantity m corresponds to the spatial (dynamic) momentum density (that is, right translation 
o/P to the identity), while M corresponds to the conserved "body" momentum density (that is, 
left translation of P to the identity). 

29.3 Brief summary 

Emb(5', M") admits two group actions. These are: the group Diff(S') of diffeomorphisms of S, which 
acts by composition on the right; and the group Diff(R") which acts by composition on the left. The 
group Diff(M") acting from the left produces the singular solution momentum map, Jsing- The action 
of Diff(5) from the right produces the conserved momentum map '■ T*Emb(S', M") X{S)* . We 
now assemble both momentum maps into one figure as follows: 



T* Emb(5, M) 




X(M)* X{S)* 



30 The Euler-Poincare framework: fluids a la [HoMaRal998a] 

Almost all fluid models of interest admit the following general assumptions. These assumptions form 
the basis of the Euler-Poincare theorem for Continuua that we shall state later in this section, after 
introducing the notation necessary for dealing geometrically with the reduction of Hamilton's Principle 
from the material (or Lagrangian) picture of fluid dynamics, to the spatial (or Eulerian) picture. This 
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theorem was first stated and proved in [HoMaRal998a], to wliicli we refer for additional details, as 
well as for abstract definitions and proofs. 



Basic Eissumptions underlying the Euler-Poincare theorem for continua 

• There is a right representation of a Lie group G on the vector space V and G acts in the natural 
way on the right on TG x V*: {Ug, a)h = {Ugh, ah). 

• The Lagrangian function L : TG x V* ^ M. is right G-invariant under the isotropy group of 

ao G V*}^ 

• In particular, if oq G V* , define the Lagrangian La,, : TG ^ M by LaQ^Ug) = L(Ug,a()). Then 
Lag is right invariant under the lift to TG of the right action of Gag on G, where Ga^ is the 
isotropy group of 0,9. 

• Right G-invariance of L permits one to define the Lagrangian on the Lie algebra q of the group 
G. Namely ^ : x F* ^ M is defined by 

i{u, a) = L{Ugg~\t),aog-\t)) = L{Ug, ao) , 

where u = Ugg~^{t) and a = aQg~^{t) . Conversely this relation defines for any £ : q x V* ^ R 
a right G-invariant function L : TG x V* ^R. 

• For a curve g{t) G G, let u{t) := g{t)g{t)~^ and define the curve a{t) as the unique solution 
of the linear differential equation with time dependent coefficients a{t) = —a{t)u{t), where the 
action of an element of the Lie algebra it G 5 on an advected quantity a e V* is denoted by 
concatenation from the right. The solution with initial condition a(0) = ao eV* can be written 
as a{t) = aog{t)~^. 



Notation for reduction of Hamilton's Principle by symmetries 

• Let g{V) denote the space of vector fields on V of some fixed differentiability class. These vector 
fields are endowed with the Lie bracket given in components by (summing on repeated indices) 

The notation adu v := [u, v] formally denotes the adjoint action of the right Lie algebra of 
Diff (P) on itself. 

• Identify the Lie algebra of vector fields g with its dual g* by using the L^ pairing 

(u,v) = / u-vdV. (30.2) 
Jv 

• Let q{T>)* denote the geometric dual space of g{T>), that is, q{T>)* := A^{V) (g) Den(D). This is 
the space of one-form densities on V. If m ^dV G A''"(P) (8)Den(P), then the pairing of mi^dV 
with u G 5(1?) is given by the L^ pairing, 

{m^dV,u)= [ m udV (30.3) 

Jv 

where m • u is the standard contraction of a one-form m with a vector field u. 



^^For fluid dynamics, right G-invariance of the Lagrangian function L is traditionally called "particle relabeling 
symmetry." 
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• For u G 0(^5) and m (g) dV G 0(^5)*, the dual of the adjoint representation is defined by 

(ad*(m(g)dy),v) = - /m-aduVdF = - / m-[u,v]dF (30.4) 

Jv Jv 

and its expression is 

ad* (m ® dV) = (£„«! + {diYdv u)m) ^ dV = £u{ni (8) dF) , (30.5) 

where div^i/u is the divergence of u relative to the measure dV, that is, £udV = (divd\/u)(iy. 
Hence, ad* coincides with the Lie-derivative £u for one-form densities. 

• If u = u^d/dx^ , m = midx^, then the one-form factor in the preceding formula for ad* (m® dy) 
has the coordinate expression 

(\ ■ I ■ din ' duP \ 

ad* mj dx^ = iw' -^-j + rrij -^—r + {divdv u)mi 1 dx^ (30.6) 

The last equality assumes that the divergence is taken relative to the standard measure dV = 
in M". (On a Riemannian manifold the metric divergence needs to be used.) 

Conventions and terminology in continuum mechanics 

Throughout the rest of the lecture notes, we shall follow [HoMaRal998a] in using the conventions and 
terminology for the standard quantities in continuum mechanics. 

Definition 30.1 Elem,ents ofT> representing the material particles of the system are denoted by X; 
their coordinates X^,A = l,...,n may thus be regarded as the particle labels. 

• A configuration, which we typically denote by r), or g, is an element of Diff (P) . 

• A motion, denoted as rjt or alternatively as g{t), is a time dependent curve in DifF(P). 

Definition 30.2 The Lagrangian, or material velocity U(X, t) of the continuum along the m,otion 
rjt or g{t) is defined by taking the time derivative of the motion keeping the particle labels X fixed: 



X 



These are convenient shorthand notations for the time derivative at fixed Lagrangian coordinate X. 

Consistent with this definition of material velocity, the tangent space to Diff (I?) at r] & Diff (I?) is 
given by 

r^Diff(P) = {Vr,:V^TV\ U^(X) e r^(x)2^}- 

Elements o/T^ Diff (P) are usually thought of as vector fields on V covering 77. The tangent lift of 
right translations on rDiff(P) by cp E Diff(I') is given by 
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Definition 30.3 During a motion rjt or g{t), the particle labeled by X describes a path in V, whose 
points 

x{X,t) -VtiX) := 9{t)-X, 

are called the Eulerian or spatial points of this path, which is also called the Lagrangian trajec- 
tory, because a Lagrangian fluid parcel follows this pa,th in space. The derivative u(x, t) of this path, 
evaluated at fixed Eulerian point x, is called the Eulerian or spatial velocity of the system: 



u{x,t) := u{rit{X),t) ■.= \J{X,t) := ^ 



Vt{X) := g{t)-X:=m9'\t)-x. 

X 



Thus the Eulerian velocity u is a time dependent vector field on T>, denoted as Uj G 0(2?), where 
Ut(x) := n{x,t). We also have the fundamental relationships 

Ut = uj o r]t and Ut = g{t)g~^{t) , 

where we denote \Jt{X) := \J{X,t). 

Definition 30.4 The representation space V* o/DifF(P) in continuum mechanics is often some 

subspace of the tensor field densities on V, denoted as T(P) (X)Dcn(P), and the representation is given 
by pull back. It is thus a right representation o/DifF(P) on T(X') 0Den(T>). The right action of the 
Lie algebra q{T>) on V* is denoted as concatenation from the right. That is, we denote 

au := £^0. , 

which is the Lie derivative of the tensor field density a along the vector field u. 
Definition 30.5 The Lagrangian of a continuum mechanical system is a function 

L : rDiff(P) X l^* ^ M, 

which is right invariant relative to the tangent lift of right translation o/Diff('D) on itself and pull back 
on the tensor field densities. Invariance of the Lagrangian L induces a function i : 5{T>) x V* ^ R. 
given by 

£{u, a) = L{u o r), rfa) = L(U, oq) , 

where u G 0(1') and a & V* C ^■{T>) ® Den(X'), and where rfa denotes the pull back of a by the 
diffeomorphism rj and u is the Eulerian velocity. That is, 

U = u o 77 and oq = ri*a . (30.8) 

The evolution of a is by right action, given by the equation 

a=—£ua = —au. (30.9) 

The solution of this equation, for the initial condition ao, is 

ait) = rit*ao = aog~^{t) , (30.10) 

where the lower star denotes the push forward operation and rjt is the flow of u = gg~^{t). 

Definition 30.6 Advected Eulerian quantities are defined in continuum mechanics to be those 
variables which are Lie transported by the flow of the Eulerian velocity field. Using this standard 
terminology, equation (30.9), or its solution (30.10) states that the tensor field density a{t) (which 
may include mass density and other Eulerian quantities) is advected. 
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Remark 30.7 (Dual tensors) As we mentioned, typically V* C 1(1?) Den(P) for continuum 
mechanics. On a general manifold, tensors of a given type have natural duals. For example, symmetric 

covariant tensors are dual to symmetric contravariant tensor densities, the pairing being given by the 
integration of the natural contraction of these tensors. Likewise, k-forms are naturally dual to (n — k)- 
forms, the pairing being given by taking the integral of their wedge product. 

Definition 30.8 The diamond operation o between elements of V and V* produces an element of 
the dual Lie algebra q{V)* and is defined as 

{boa,w) = - [ b-£^a, (30.11) 
Jv 

where b ■ £^ a denotes the contraction, as described above, of elements of V and elements of V* and 
w G g(I'). (These operations do not depend on a Riemannian structure.) 

For a path r]t G Diff (P), let u(x, t) be its Eulerian velocity and consider the curve a{t) with initial 
condition oq given by the equation 

a + £ua = 0. (30.12) 

Let the Lagrangian La„i\J) := L(U,ao) be right -invariant under Diff(P). We can now state the 
Euler-Poincarc Theorem for Continua of [HoMaRal998a]. 

Theorem 30.9 (Euler Poincare Theorem for Continua.) Consider a path rjt in Diff(P) with 
Lagrangian velocity U and Eulerian velocity u. The following are equivalent: 

i Hamilton's variational principle 

6 f^L{X, Vt{X),ao{X)) dt = (30.13) 

Jti 

holds, for variations Srjt vanishing at the endpoints. 

ii r)t satisfies the Euler-Lagrange equations for LaQ on Diff (P) . 

iii The constrained variational principle in Eulerian coordinates 

'■t-2 



'tl 

holds on 0(P) X V* , using variations of the form 



6 j \{u,a) dt = (30.14) 



(9w (9w 

(^u = ^r- + [u, w] = — - + aduW , 6a=-£y/a, (30.15) 
at ot 

where Wt = Srjt o rj^^ vanishes at the endpoints. 

iv The Euler-Poincare equations for continua 

d Si „ Si Si ^ Si Si 

hold, with auxiliary equations [dt + £u)o, = for each advected quantity a{t). The o operation 
defined in (30.11) needs to be determined on a case by case basis, depending on the nature of the 
tensor a{t). The variation m = S£/6u is a one-form density and we have used relation (30.5) 
in the last step of equation (30.16). 

We refer to [HoMaRal998a] for the proof of this theorem in the abstract setting. We shall see 
some of the features of this result in the concrete setting of continuum mechanics shortly. 
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Discussion of the Euler-Poincare equations 

The following string of equalities shows directly that iii is equivalent to iv: 
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dtSu "(5u Sa 



dt. 



(30.17) 



The rest of the proof follows essentially the same track as the proof of the pure Euler-Poincare theorem, 
modulo slight changes to accomodate the advected quantities. 

In the absence of dissipation, most Eulerian fluid equations''^^ can be written in the EP form in 
equation (30.16), 

d 6£ d£ 51 

— — + ad*— = — oa, with idt + £u)a = Q. (30.18) 
ot5n 5vl 5a ^ ' 

Equation (30.18) is Newton's Law. The Eulerian time derivative of the momentum density m = 
51/ 5u (a one-form density dual to the velocity u) is equal to the force density {51 /5a) o a, with the o 
operation defined in (30.11). Thus, Newton's Law is written in the Eulerian fiuid representation as,^^ 



d^ 
dt 



S£ 



Lag 



m := [dt + =Cu)m = — oa 



with — 
dt 



Lag 



:= {dt + £u)a = 0. 



(30.19) 



The left side of the EP equation in (30.19) describes the fluid's dynamics due to its kinetic energy. 
A fluid's kinetic energy typically deflnes a norm for the Eulerian fluid velocity, KE = |||u|p. 
The left side of the EP equation is the geodesic part of its evolution, with respect to this norm. 
See [ArKhl998] for discussions of this interpretation of ideal incompressible flow and references 
to the literature. However, in a gravitational fleld, for example, there will also be dynamics due 
to potential energy. And this dynamics will by governed by the right side of the EP equation. 

The right side of the EP equation in (30.19) modifles the geodesic motion. Naturally, the right 
side of the EP equation is also a geometrical quantity. The diamond operation o represents the 
dual of the Lie algebra action of vectors fields on the tensor a. Here 51/ 5a is the dual tensor, 
under the natural pairing (usually, pairing) ( • , • ) that is induced by the variational derivative 
of the Lagrangian £(u, a). The diamond operation o is defined in terms of this pairing in (30.11). 
For the L^ pairing, this is integration by parts of (minus) the Lie derivative in (30.11). 

The quantity a is typically a tensor (e.g., a density, a scalar, or a differential form) and we shall 
sum over the various types of tensors a that are involved in the fiuid description. The second 
equation in (30.19) states that each tensor a is carried along by the Eulerian fluid velocity u. 



^''Exceptions to this statement arc certain multiphase fluids, and complex fluids with active internal degrees of freedom 
such as liquid crystals. These require a further extension, not discussed here. 

^^In coordinates, a one-form density takes the form m ■ dx g) dV and the EP equation (30. 16) is given neumonically by 



dt I Lag 



(m • dx (gi dV) 



dt 



Lag 



dx (g) -I- m-du^dV + m • dx ® (V • u)dV 



Sa 



o a 



Advection 



Stretching 



Expansion 



with 



dx 

I -Lag 



[dt + i^u)dx = du = ujdx^ , upon using commutation of Lie derivative and exterior derivative. Compare 



this formula with the definition of adu(m (g> dV) in equation (30.6). 



D. D. Holm 



Imperial College London 



Applications of Poisson Geometry 



122 



Thus, a is for fluid "attribute," and its Eulerian evolution is given by minus its Lie derivative, 
— XuCi. That is, a stands for the set of fluid attributes that each Lagrangian fluid parcel carries 
around (advects), such as its buoyancy, which is determined by its individual salt, or heat 
content, in ocean circulation. 

• Many examples of how equation (30.19) arises in the dynamics of continuous media are given in 
[HoMaRal998a]. The EP form of the Eulerian fluid description in (30.19) is analogous to the 
classical dynamics of rigid bodies (and tops, under gravity) in body coordinates. Rigid bodies 
and tops are also governed by Euler-Poincare equations, as Poincare showed in a two-page paper 
with no references, over a century ago [Pol 901]. For modern discussions of the EP theory, see, 
e.g., [MaRal994], or [HoMaRal998a]. 

Exercise 30.10 For what types of tensors ao can one recast the EP equations for continua (30.16) 
as geodesic motion, by using a version of the Kaluza-Klein construction ? 



30.1 Corollary of the EP theorem: the Kelvin-Noether circulation theorem 

Corollary 30.11 (Kelvin-Noether Circulation Theorem.) Assume u{x,t) satisfies the Euler- 
Poincare equations for continua: 

and the quantity a satisfies the advection relation 

da 

— + £uO = 0. (30.20) 

Let Tjt be the flow of the Eulerian velocity field u, that is, u = {drjt / dt) orj^^ . Define the advected fluid 
loop 7t := rjt o 70 and the circulation map I{t) by 

m^iyi- (30.21) 

In the circulation map I{t) the advected mass density Df satisfies the push forward relation Df = ti^Dq. 
This implies the advection relation (30.20) with a = D, namely, the continuity equation, 

dtD + dwDu = 0. 

Then the map I{t) satisfies the Kelvin circulation relation, 

d , , f 1 5i , , 

I(t) = d) 00 . (30.22 

Both an abstract proof of the Kelvin-Noether Circulation Theorem and a proof tailored for the case of 
continuum mechanical systems are given in [HoMaRal998a]. We provide a version of the latter below. 



Proof. First we change variables in the expression for I{t): 



m 



It 



1 61 



Vt 



70 



1 61 



70 



-Vt 



'61' 
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Next, we use the Lie derivative formula, namely 



applied to a one-form density at- This formula gives 



dt 



m 



dt 



1 

70 ^0 

l^df. 



51 
5x1 



Vt 



51 



To 



Do 



Tit 



(5u 
dt [ja 



5[_ 
Ja 



'51 




'51 ' 


— oa 




^ oa 


da 




da 



By the Euler-Poincare equations (30.16), this becomes 
again by the change of variables formula. ■ 

Corollary 30.12 Since the last expression holds for every loop 74, we may write it as 

d ^ \ I 51 I 51 



dt 



£ 



D5u D5a 



o a . 



(30.23) 



Remark 30.13 The Kelvin-Noether theorem is called so here because its derivation relies on the 
invariance of the Lagrangian L under the particle relabeling symmetry, and Noether's theorem is 
associated with this symmetry. However, the result (30.22) is the Kelvin circulation theorem: the 
circulation integral I{t) around any fluid loop (jt, moving with the velocity of the fluid parcels u) is 
invariant under the fluid motion. These two statements are equivalent. We note that two velocities 
appear in the integrand I{t): the fluid velocity u and D^^5£/5u. The latter velocity is the momentum 
density m = 51/ 5u divided by the mass density D. These two velocities are the basic ingredients for 
performing modeling and analysis in any ideal fluid problem. One simply needs to put these ingredients 
together in the Euler-Poincare theorem and its corollary, the Kelvin-Noether theorem. 

31 Euler-Poincare theorem & GFD (geophysical fluid dynamics) 
31.1 Vciriational Formulae in Three Dimensions 

We compute explicit formulae for the variations 5a in the cases that the set of tensors a is drawn from 
a set of scalar fields and densities on M^. We shall denote this symbolically by writing 



aG{b,D d^x} 



(31.1) 



We have seen that invariance of the set a in the Lagrangian picture under the dynamics of u implies 
in the Eulerian picture that 

d \ 
Q^ + £u a = 0, 
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where £u denotes Lie derivative with respect to the velocity vector field u. Hence, for a fluid dynamical 
Eulerian action & = J dt i{u; b, D), the advected variables b and D satisfy the following Lie-derivative 
relations, 



+ £ii] b = 0, or 



d_ 
dt 



dt 



+ £^\Dd^x = Q, or 



db 
dt " 
dD 

~dt 



u- V6, 
- V • (Du) . 



(31.2) 
(31.3) 



In fluid dynamical applications, the advected Eulerian variables b and D d^x represent the buoyancy 
b (or specific entropy, for the compressible case) and volume element (or mass density) D d^x, respec- 
tively. According to Theorem 30.9, equation (30.14), the variations of the tensor functions a at fixed 
X and t are also given by Lie derivatives, namely 5a = — £w or 

5b = ~ £^ff b = — -w • V b , 
5D d^x = -£y,{D d^x) = - V • (L>w) d^x . (31.4) 

Hence, Hamilton's principle (30.14) with this dependence yields 







5 J dt 

/* 
/* 

J dtw 
- J dtw 



n;b,D) 



M 
M 

5u \dt 

d_5i_ 
~ diJa 



/dw \ 6£ 



w-Vb- 



Si_ 
SD 



(V • {Dw)) 



5£ 51 ^^51 



5£ ^, ^^5£ 

1 1 Vb- D V — 

)5vi 5b 5D 



(31.5) 



where we have consistently dropped boundary terms arising from integrations by parts, by invoking 
natural boundary conditions. Specifically, we may impose n • w = on the boundary, where fi is the 
boundary's outward unit normal vector and w = 5rito rj^^ vanishes at the endpoints. 



31.2 Euler Poincare framework for GFD 

The Euler-Poincare equations for continua (30.16) may now be summarized in vector form for advected 
Eulerian variables a in the set (31.1). We adopt the notational convention of the circulation map / 
in equations (30.21) and (30.22) that a one form density can be made into a one form (no longer a 
density) by dividing it by the mass density D and we use the Lie-derivative relation for the continuity 
equation {d/dt + £u)Dd^x = 0. Then, the Euclidean components of the Euler-Poincare equations 
for continua in equation (31.5) arc expressed in Kelvin theorem form (30.23) with a slight abuse of 
notation as 

d „\/15i .\ 1 5£_, , _/6£\ , ^ 

0, (31.6) 



+ 



-V6 • dx 



^D5u J D5b 

in which the variational derivatives of the Lagrangian £ arc to be computed according to the usual 

physical conventions, i.e., as Frechet derivatives. Formula (31.6) is the Kelvin-Noether form of the 
equation of motion for ideal continua. Hence, we have the explicit Kelvin theorem expression, cf. 
equations (30.21) and (30.22), 

d I 1 5e , r 1 5£^, , , 

i -^■d-K = -j> --Vb-d^, 31.7 
J^^i^)D5u .UhAD5b 



dt 



7t{u) 
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where the curve 7t(u) moves with the fluid velocity u. Then, by Stokes' theorem, the Euler equations 
generate circulation of v := {D~^Sl/du) whenever the gradients V6 and V{D~^dl/db) are not collinear. 
The corresponding conservation of potential vorticity q on fluid parcels is given by 

da 1 / 1 5i\ 

-^ + u-Vg = 0, where g = — V6 • curl — — . (31.8) 
ot D \D on J 

This is also called PV convection. Equations (31.6-31.8) embody most of the panoply of equations 
for GFD. The vector form of equation (31.6) is, 

fd \fl5l\ I 51 ^ ^61 ISl, 

.. ■' ■' 



Geodesic Nonlinearity: Kinetic energy Potential energy 

In geophysical applications, the Eulerian variable D represents the frozen- in volume element and b is 
the buoyancy. In this case, Kelvin's theorem is 

with circulation integral 

31.3 Euler's Equations for a Rotating Stratified Ideal Incompressible Fluid 

The Lagrangian. In the Eulerian velocity representation, we consider Hamilton's principle for fluid 
motion in a three dimensional domain with action functional S = J Idt and Lagrangian Z(u, b, D) given 

by 

l{u, b,D) = J pqD{1 + b) Q|up + u • R(x) - - p{D - 1) d^x , (31.10) 

where ptot = PoD{l + b) is the total mass density, po is a dimensional constant and R is a given 
function of x. This variations at flxed x and t of this Lagrangian are the following. 



51 1 

po(l + b){n + R) , - = /Jo^(-|up + u • R - ffz) , 



^ = po(l + 6)(2|up + u-R-5^)-p, ^1 = -{D-1). (31.11) 

Hence, from the Euclidean component formula (31.9) for Hamilton principles of this type and the 
fundamental vector identity, 

(b • V)a + ajVb> = - b x (V x a) + V(b • a) , (31.12) 

we find the motion equation for an Euler fiuid in three dimensions, 

^-uxcurlR + 5Z+— ^Vp = 0, (31.13) 

at po(l + o) 

where curlR = 20 (x) is the Coriolis parameter (i.e., twice the local angular rotation frequency). In 
writing this equation, we have used advection of buoyancy, 

- + u.V6 = 0, 
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from equation (31.2). The pressure p is determined by requiring preservation of the constraint D = 1, 
for which the continuity equation (31.3) imphes divu = 0. The Euler motion equation (31.13) is 
Newton's Law for the acceleration of a fluid due to three forces: CorioHs, gravity and pressure gradient. 
The dynamic balances among these three forces produce the many circulatory flows of geophysical 
fluid dynamics. The conservation of potential vorticity q on fluid parcels for these Euler GFD 
flows is given by 

dq 

— + u • Vq = , where, on using D = 1 , g = V6 • curl(u + R) . (31.14) 

Exercise 31.1 (Semidirect-product Lie-Poisson bracket for compressible ideal fluids) 

1. Compute the Legendre transform for the Lagrangian, 

l{u,b,D) : Xx A°x A^H^M 
whose advected variables satisfy the auxiliary equations, 

2. Compute the Hamiltonian, assuming the Legendre transform is a linear invertihle operator on 
the velocity u. For definiteness in computing the Hamiltonian, assume the Lagrangian is given 
by 

l{u,b,D) = J D (^^{ul"^ + u-'R{x) - e{D,b)^ d^x, (31.15) 

with prescribed function R(x) and specific internal energy e{D, b) satisfying the First Law of 
Thermodynamics, 



de = ^dD + Tdb , 



where p is pressure, T temperature. 



3. Find the semidirect-product Lie-Poisson bracket for the Hamiltonian formulation of these equa- 
tions. 

4. Does this Lie-Poisson bracket have Casimirs? If so, what are the corresponding symmetries and 
momentum maps? 

32 Hamilton-Poincare reduction and Lie-Poisson equations 

In the Euler-Poincare framework one starts with a Lagrangian defined on the tangent bundle of a Lie 
group G 

L:TG^R 

and the dynamics is given by Euler-Lagrangc equations arising from the variational principle 

r-ti 



5 [' L{g,g)dt = 



The Lagrangian L is taken left/right invariant and because of this property one can reduce the problem 
obtaining a new system which is defined on the Lie algebra g of G, obtaining a new set of equations, 
the Euler-Poincare equations, arising from a reduced variational principle 



'to 

where Z(^) is the reduced lagrangian and ^ G fl. 



6 [ ' l{^)dt = 

Jtn 
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Problem 32.1 Is there a similar procedure for Hamiltonian systems? More precisely: given a Hamil- 
tonian function defined on the cotangent bundle T*G 

H : T*G R 

one wants to perform a similar procedure of reduction and derive the equations of motion on the dual 
of the Lie algebra g* , provided the Hamiltonian is again left/right invariant. 

Hamilton-Poincare reduction gives a positive answer to this problem, in the context of variational 
principles as it is done in the Euler-Poincare framework: we are going to explain how this procedure 
is performed. 

More in general, we will also consider advected quantities belonging to a vector space V on which G 
acts, so that the Hamiltonian is written in this case as [?] [HoMaRa98] 

H : T*G X ^* ^ M 



The space V is regarded here exactly the same as in the Euler-Poincare theory. 

The equations of motion, i.e. Hamilton's equations, may be derived from the following variational 
principle 

6 [ \{pit),git)) - HaMt),Pm dt = 
Jto 

as it is well know from ordinary classical mechanics {g{t) has to be considered as the tangent vector 
to the curve g{t), so that g{t) G Tg^^-^G). 

Problem 32.2 What happens if Ha^ is left/right invariant? 

It turns out that in this case the whole function 

Fi9,9,P) = {P,9) - Haoig,p) 

is also invariant. The proof is straightforward once the action is specified (from now on we consider 
only left invariance): 

h{g,g,p) = {hg, TgLhg, TlgL^-ip) 

where TgL^ '■ TgG ThgG is the tangent of the left translation map L^g = hg E G sit the point g 
and T*gLf,-i ■.T*G^ T*gG is the dual of the map ThgL^-i : ThgG TgG. 
We now check that 

{hp, kg) = {TlgLy,-ip, TgLhg) 
= {p, ThgLh-i o TgLhg) 
= {P, Tg{Lh-i o Lh) g) = {p,g) 

where the chain rule for the tangent map has been used. The same result holds for the right action. 
Due to this invariance property, one can write the variational principle as 



S f\{li,i)-h{ii,a)}dt = {) 
with 

^i{t) = g-\t)p{t) e g*, m = 9~\t)g{t) e g, a{t) = g-\t)ao e V* 
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In particular a{t) is the solution of 

a{t) = -^{t) ao. 

where a Lie algebra action of q on V* is implicitly defined. In order to find the equations of motion 
one calculates the variations 

^[''iil^^O - h{fi,a)}dt = S^{dfi,0 + {fi,SO - (Sfi, - (da, dt 
As in the Euler-Poincare theorem, we use the following expressions for the variations 

SC = V+[C,V], Sa = -r]a 
and using the definition of the diamond operator we find 



to 



Sfi / \ Sa 
Sfi, - + {n,rj + ad^ri) + (rja, dt 



so that 



and the equations of motion are 
together with 



Sfj, 



ii = adt u — — o a 
^ oa 



5h 

a = —— a. 



This equations of motion written on the dual Lie algebra g are called Lie-Poisson equations. We have 
now proven the following 

Theorem 32.3 [Hamilton- Poincare reduction theorem] 

With the preceding notation, the following statements are equivalent: 



1. With ao held fixed, the variational principle 

I to 

holds, for variations Sg{t) of g{t) vanishing at the endpoints. 
^- {9{t)jP{t)) satisfies Hamilton's equations for Ha^^ on G. 



5 f\{p{t),g{t)) - HaMt),Pmdt = 

Jto 
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3. The constrained variational principle 

rti 



S [\{fi{t),m-h{l^{t),amdt = 

Jtn 



'to 

holds for g X V* , using variations of ^ and a of the form 

S^ = r)+[^,V], 6a=-r]a 
where rj{t) G Q vanishes at the endpoints 
4- The Lie-Poisson equations hold on Q x V* 



{fi, a) = ^ad|/x 



Sh 



o a, 



Sh 



Sa ' Sfi 

Remark 32.4 More exactly one should start with an invariant Hamiltonian defined on 

T*{G xV) = T*G xV xV* 

However, as mentioned in [][HoMaRe98], such an approach turns out to he equivalent to the treatment 
presented here. 

Remark 32.5 [Legendre transform] 

Lie-Poisson equations may arise from the Euler-Poincare setting by Legendre transform 

SI 



If this is a diffeomorphism, then the Hamilton-Poincare theorem is equivalent to the Euler-Poincare 
theorem. 

Remark 32.6 [Lie-Poisson structure] 

One shows that q* xV* is a Poisson manifold: 



/ ^ 



SH SF 

Sn ' Sfj, 
'SF SH 

Sfi ' Sh 

'SF SH' 
Sfj, ' Sfx 

In fact it can be easily shown that this structure 

'SF SH 



SF\_/SH SF\_ 
Sfi / \ SjjL ^ Sa 
SH SF 



Sa ' S/j, 



SH 
Sa 



o a, 



/SH SF 
\ Sfi ^ Sa 
SH SF 
Sjj, ' Sa 



a, 



{F,H}{i^,a) = -(i^, 



Sfj, ' Sfi 



SF SH SH SF 
Sfi Sa Sfj, Sa 

SF SH SH SF 
' Sfj Sa Sfi Sa 



- { a 



satisfies the definition of a Poisson structure. In particular one finds that any dual Lie algebra Q 
is a Poisson manifold. 



Please note: this structure has been found during lectures for the simpler case without advected 
quantities. 
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Remark 32.7 [right invariance] 

It can be shown that for a right invariant Hamiltonian one has 



+ a 



{fi, a) = - |^ad|/x - 
with all signs changed respect to the case of left invariance. 



6F 5H 

6h 
da 



6F 6H 6H 5F 
S/j, 6a Sji 5a 



6h 

o a, -— a 

OjJL 



33 Two applications 
33.1 The Vlasov equation 

In plasma physics a main topic is colHsionless particle dynamics, whose main equation, the Vlasov 
equation, will be heuristically derived here. In this context a central role is held by the distribution 
function on phase space /(q, p,t), basically expressing the particle density on phase space. Intended 
as a density one defines F := f{c[,p,t)dc[dp: because of the conservation of particles, one writes the 
continuity equation just as one does as in the context of fluid dynamics 

F + V • (uF) = 

where u is a "velocity" vector field on phase space, which is given by the single particle motion 

u = (q,p) G X(T*M^) 

if we now assume that the generic single particle undergoes a Hamiltonian motion, the Hamiltonian 
function /i(q, p,i) can be introduced directly by means of the single particle Hamilton's equations 



(4p) 



dh dh 
5p' 5q 



which shows that u has zero divergence, assuming the Hessian of h is symmetric. Therefore, the 
Vlasov equation written in terms of the distribution function /(q, p,t) is 

/ + u-V/ = 

Expanding now the Hamiltonian h as the total single particle energy 

h{(\,p,t) = ^p'^ + V{(i,p,t) 

one obtains the more common form 

^+ P .^_^.^ = 
dt m d(\ d(\ (?p 

Problem 33.1 Can Vlasov equation he cast in Lie-Poisson form? 
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We show here why the answer is yes. First we write the Vlasov equation in terms of a generic single 
particle Hamiltonian h as 

f + {f,h} = 
where we recall the canonical Poisson bracket 

' dq dp dp dq 

The main point of this discussion is that the canonical Poisson bracket provides the set J^{T*W^) of 
the functions on the phase space with a Lie algebra structure 

[k,h] = {k,h} 

At this point, in order to look for a Lie- Poisson equation, one calculates the coadjoint operator such 
that 

{f,{h,k}) = ifMhk) = W,k) = {-{hj},k) 

where the last equality is justified by the Leibniz property of the Poisson bracket, with the pairing 
defined as 

if, 9) = J fgdcidp. 

In conclusion, the argument above shows that the Vlasov equation can in fact be written in the Lie 
Poisson form 

/ + ad;;/ = o 



33.2 Ideal barotropic compressible fluids 

The reduced Lagrangian for ideal compressible fluids is written as 



liu,D) = J ^\u\'^ - De{D) dx 



where u G X(McM^) is tangential on the boundary dM and D is the advected density, which satisfies 
the continuity equation 

dtD + jC^D = 0. 

Moreover, the internal energy satisfies the barotropic First Law of Thermodynamics 

de = -p{D)d{D-^) = dD 

for the pressure p{D). The "reduced" Legendre transform on this Lie algebra j£(M^) is given by 

m = Du 

and the Hamiltonian is then written as 

h{m,D) = (m,u) - l{u,D) 

that is 



h{m,D) = J ^\mf + De{D)dy: 
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The Lie Poisson equations in this case are as from the general theory 

dtm = -ad*s^fsm - 1^ ^ ^ 

dtD = —Csh/Sm D 

Earher we found that the coadjoint action is given by the Lie derivative. On the other hand we may 
calculate the expression of the diamond operation from its definition 



to be 



Therefore, we have 



5h ^ \ / Sh ^ \ 



-,-dwD,) = (w^,, 



Sh ^ Sh 
— oD = DV — 
SD SD 



where 

Substituting into the momentum equation and using the First Law to find d{e + p/D) = {1/D)dp 
yields 

dtm. = —jC-u m — 

Upon expanding the Lie derivative for the momentum density m and using the continuity equation 
for the density, this quickly becomes 

dtu = — u- Vu — 
which is Euler's equation for a barotropic fluid. 

33.3 Euler's equations for ideal incompressible fluid motion 

The barotropic equations recover Euler's equations for ideal incompressible fluid motion when the 
internal energy in the reduced Lagrangian for ideal compressible fluids is replaced by the constraint 
D = l,as 

l{n,D) = j ^|u|2-p(D-l)dx 



where again u G X(Mc M ) is tangential on the boundary dM and the advected density D satisfles 
the continuity equation, 

dtD + AwDu = 0. 

This equation enforces incompressibility divu = when evaluated on the constraint D = 1. The 
pressure p is now a Lagrange multiplier, which is determined by the condition that incompressibility 
be preserved by the dynamics. 
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